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Abstract 


This fvixm pmt'iits n survey of the i/uvry of the twtifs ivtaliou tmtl the^eophyshvl 
phQUomma ajfeviim it^ ^eith emphasis on polar motion ami mint ions. The tlteoreth 
ail development in this review be0s with first principles and formulates the problem of 
(Hilar motion and UTI mriations in eonsidemble generality and detail. The traitment 
includes an analysis of the effects of earth deformations and the solid earth tide.s. 
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A Survey of the Theory of the Earth’s Rotation 


I. Introduction 

In recent years the tccliniqiie of long baseline inter* 
feronietry (UBI) lias deinonstratcd its ability to carry out 
incasureinenls of the variations in tlic earth’s rotation rate 
(UTl) and of the variations in the |)osltlon of the earth’s 
rotation axis with respect to the earth’s crust (polar niotlon) 
wlili an accuracy whicli is comparable to tiiat obtained by tire 
classical methods using pliotographic zenltli tubts (PZT’s). 
While the classical methods have probably been extended to 
tlieir maximum capabilities by their present day use, it Is 
expected that coniimied development of the technique of long 
baseline inicrferornetry will ultimately yield measurements of 
these geopliyslcal quantities with accuracies wliicli exceed 
present day capabilities by an order of magnitude. 

Tliis document was prepared in anticipation of these future 
dcvclopmcms, its purpose is to compreliensively review llic 
present day dynamical tlicory of tlic eartlv’s rotation in order 
to provide a coherent liicorelical basis for tlie development of 
future data analysis ))roccdures and software models for tlie 
treatment of l\iiurc high quality long baseline interferometry 
data. Its purpose ihon is to serve as a tutorial lumdbook for 
workers who will be involved in tlic process of extracting 
gcopityslcal informal 'on from tlic interferometry data. In 
order to cnlumce its usefulness in tins respect the work is 
presented in considerable detail. 

II, Coordinate Systems 

The term “axis” will always be associated witli a corre- 
sponding vector. An axis is a stiaiglu line passing Ihroiigli tlie 


origin of coordinates in a direction parallel to the associated 
vector. TJic point of intersection of an axis wltli the surface of 
the eartli or the celestial splicrc is called a ’’pole”} in tlie 
former case a “terrestriai pnie” and in the latter case a 
“celestial pole.” If not explicitly siuied wldch, “terresirtnl 
pole” will be understood. 

In cstablisidng a theoretical and operational framework for 
describing tlie earth’s rotation it is necessary to use two 
coordinate frames: a space* fixed frame spanned by basis 
vectors ^3 and a body-llxed frame spannod by basis 
vectors ?, Coordinate .systems can be grouped into 

"geometrical” and “dynamical" classifications according to 
tlic iiatiire of their fuiulamcnial denning quantities. Hybrid 
coordiiuile systems requiring a combination of gcomctiieal ami 
dynamical quantities for tlicir dcllnition are also po.ssiblc. 

Tile present space*fi,xed Coordinate frame is a dynamical 
coordlnntc frame which uses tlic orbital and eqtiatorial planes 
of the eartli to define /fj ^ 3 . Tliei ’^3 axis is parallel to the 
earlli’s mean orbital angular momentum vector of 1950.0 and 
i?i is contained by the inieiscclion of i|ie mean orbital and 
mean equatorial planes of 1950.0 and points toward the 
ascending node. Tlic I?, axis is ortliogonal to /?, and to 
complete a rlglit-lianded coordinate frame. Tlic origin of tlic 
space-Oxed coordinate frame is placed at tlie center of mass of 
the eartli, including the oceans and atmosplierc. In defining 
the mean orbiial plane of the eartli it Is necessary to reckon 
willi die fact Uiat tile motion of tlie eartli about tlie sun is 
contimially perturbed by tlie gravitational attract ions of the 
ulhcr bodies of the solar system and tliat tlie actual orbit of 


ihtf earf h‘s ceiUer of nuiss is an irregular ami ever varying curve 
in space (VVoolarrl ami Clemence, 

liven in tire case of a rigiil earth certain liifficuUies are 
encountered when operationally defining the celestial equator. 
The instantaneous celestial equator is usually dellned as being 
contained by a plane wldcli is orthogonal to the earth’s 
instantaneous rotation axis. Fora rigid earth the Instantaneous 
rotation axis is defined by the vector sturr of the angular 
rotation rates of polar molion, spin, precession, and nutation 
and is unobservable by conventional astronomical incairs. What 
is observable is a position on the celestial sphere, to be here 
called the celestial reference pole, located at the center of the 
quasi'Cireular diurnal paths of the stars In tire sky. The celestial 
reference pole is in eonliiuml motion across the sky owing to 
lire effects of precession and muailon alone. Polar motion 
changes the latitudes of observatories and not tire declination 
of stars and so docs not contribute to the motion of the 
celestial reference pole on the celestial sphere. Tire motion of 
Ibe celestial reference pole across lire sky implies the existence 
of an additional rotutioir rate (due to precession and nutation) 
which, ulojrg with the rotation rate due to the polar molion, 
must be added vcctorially to tire spin to obtain tire total carlli 
rotation vector. 

The reali/alion of tire above dynaniicnUy dellned space- 
llxcd basis vectors Is provided in air implicit manner by tlic 
coordinates assigned llie stars of llic FK.4 catalogue whoso 
positions luivc been measured whir reference lo the earth’s 
orbital ami equatorial planes. The FK4 catalogue contains 
aboui 1500 stars and lias an Internal consistency of iO'IiS are 
and an overall internal accuracy of 'O'llO arc (Kolac/ek and 
Wciffenbacb, p, 32. 1975), Soon to replace tlic FK4 catalogue 
is llie FK5 catalogue with about 3000 Stars and an internal 
consistency of tO'IlO arc and an overall accuracy of l' 0‘I02 arc. 

long baseline inierleromclry is expected to provide the 
relative positions of louglily 100, compact radio sourcc.s with 
an overall internal accuracy of iO'IOOl are in tlic near future, 
tins celestial coordinate system will be essentially geometrical 
and some effort sluwld be dedicated lo lying ii to the 
dynamically defined FK4 and FKS coordinate systems, at least 
to the level of llie errors inherent In the stellar coordinate 
systems. 

Tlie body.Hxed coordinate system spanned by basis vectors 
?! has its origin at the center of mass of tlic carilr 
including llic oceans and atmosplicre. lienee the space-fixed 
coordinate system and tlic body.fixcd coordinate system share 
a common origin, and the general linear coordinate trans- 
formation rclaUng llic two at any particular lime consists of a 
rotation about some axis and a scaling. A rclaUve scaling of 


unity Is maintained by tulopting the same imit of length In 
both systems. 

The basts vector is defined to be parallel to the mean 
axis of figure of the earth. The axis of figure correspouds to 
the principal eigenvector (the eigenvector of maximum eigen- 
value) of the earth’s Inertia tensor. For a deformable earth the 
body-fixed oriental Ion of the axis of figure Is Urne-dependem 
since the inertia tensor is time-dependent. 

Tlie mean axis of figure Is defined us corresponding to the 
principal eigenvector of the mean inertia tensor. 

Tlie basis vector is orthogonal to thePj basis vectonmd 
contained in the plane of the conventional prime (zero) 
meridian. The Fj basis vector is orihogonal to and 
oriented so as to form a right-handed orthogonal triad. The 
realiwUion of tlic body-fixed basis vectors F, % is provided 
implicitly by the coordinates of a set of fixed observatory sites 
located on tlic earth’s surface. The following Is a brief 
examination of current practices involved in (lie determination 
of tlie tericsiriul coordinate frame. 

In the prcsatcllite era the gcocenter was, to a large extent 
inaccessible geodctically speaking, and geodetic networks were 
essentially "local” coordinate syslcms. The location of Ibe 
gcocenter relative to the earth’s surface could be determined in 
principle by solving the boundary value problem In the theory 
of gravitational potential. However, tlie |ir,actlcal realization of 
this procedure was hampered by the lack of gravity measure- 
ments over the oceanc, and any attempt based on the limited 
data available excluded the carlii’s atmospbcrc since it was 
outside llic surface over wlilcb the measurements were made. 

The artificial satellite senses tlic center of mass of the earlli 
directly since llie osculating orbit plane passes llirougli it. Tlie 
usual procedure is to use satellite tracking data to .solve 
simultaneously for the geocentric position of the trucking 
station and the spbcricat harmonic coefficients in the expan- 
sion for tlie ciirtli’s gravitational potential. In this manner the 
location of the gcocenter relative to points on the earth’s 
surface can be determined to about i0.5m il.Om. 

Fven in the ease of a rigid earth, owing to the presence of 
the fluid porlions, the geocentcr docs not remain fixed relative 
to llie earth’s .solid surface. The seasonal redistvibulion of the 
masses of the oceans and atmosphere and particularly the 
redistribution of global ground water displace the geocentcr 
relative lo the solid earth periodically in a year along a roughly 
elliptical path v/ith u major axis of the order of 0.5cm (Stolz 
1976). 
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Considering the real deformable earilt, the tidal defonnn* 
tions of the solid portion, becouse of their symmetry, do not 
displace the geocenter, However, owing to their asymmetry 
the tidal redistribution of the masses of the oceans does 
displace the gcocenter, fiach tidal constituent will displace the 
gcocenier along a ctirved patli relative to the solid earth with a 
period equal to tliat of the tidal constituent. Brosche and 
Stmdermann (1977) have shown that theAfj tidal constitneiu 
displaces the geocenlcr relative to the solid earth antnod a 
closed curve of the order of 4 cm in extent with a period of 
12.42 hours. The motion of the geoconler due to theAfj tide 
h essentially connned to the plaiie of the equator (the 
displacement parallel to the polar axis being an order of 
magnitude smaller) along the directions X » 45“B, 1 35®W. 

It can be shown (Goldrcich and Toomre, 1969) that for a 
quasi-rigld, evolving, extended body rotating in the absence of 
external torques about its axis of largest moment of inertia, 
the axis of figure is constrained dynamically to coincide with 
the axis of rotation. However, the earth, being subjected to 
external torques, differs In some respects from the body being 
considered by Goldrcich and Toomre. 


Bttlerlan motion, with only a small perturbation of the order 
0(0102 are occurring ns a result of the lunhsolar gravitational 
torques. Although the elastic yielding of the earth greatly 
alters the Kulcrlan motion from what would prevail on a rigid 
earth, the effect of the lunkolnr torques Is simply to displace 
the instantaneous rotation ii.xls In a diurnal circular pntli of 
diameter roughly 0102 arc about the luilerian position. Thus 
the mean axis of rotation of the earth coincides with tl»c 
Kulcrlan (torque-free) axis of rotation and hence, by the 
arguments of Goldrelcli and Toomre, also coincides with the 
mean axis of figure 

In general, at any Instant the rotation axis Is displaced, in 
the body-fixed frame, from its mean position on tlie axis of 
figure, The figure axis of the earth or the basis vector can 
be determined by establishing the mean position of tire 
rotation axis, In practice the procedure of determining the 
location of the ‘'Instantaneous" rotation axis tii the body-llxed 
frame by PZT observations yields the posillon of the "imsian* 
lancous" Kulcrlan or spin axis. However, ilie mean position of 
tlio Kulcrlan axis will also .serve to determine lire carlh’s figure 
axis ortl’ 3 . 


Kor historical reasons the motion of a rigid body in the 
absence of externa! torques Is known us “Kulorian motion," In 
the ca.se of the earth subjected as it is to external torques it 
can bo shown (Woolard, 1953) that tlic motion of tire rotallon 
axis williin the body-fixed frame is almost entirely due to tlic 


Kigurc 1 1-1 (after Smylie and Mansinba, 1971b) lllustrute.s 
the geometrical relationships pertaining lo polar motion. The 
Hgurc is drawn with rcrcrenec to the Bulcrlan or spin axl.s 
rather than the rotation axis, since it is willi reference to this 
axis dial latllndgs on earth are observatlonally delcrmlned. 



Figur* 11-1. The geometry of poler motion. For clerity the figure ii drawn in a body-fixed frame «o 
that the figure of the earth defined by e, appeare fixed in orientation. Whan viewed in 
a epace-fixed frame, the equator and the apin axia appear fixed in orientation while the 
figure of the earth it diaplaced, (Attar Smylie and Manainha, 1971b) 
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In Figure IH the point 0 refers to an observer on the 
earth’s surface or alternately to the body fixed angular 
coordinates of an interferometer baseline. The point C is the 
geoccnter, the origin of tiie terrestrial coordinate system. The 
displacement of the Eulerlan or spin pole from P, the axis of 
figure or reference polej to F’ clianges the latitude of the 
observer or the declination of the interferometer baseline at 0 
from the angle OCM to the angle OCNV and changes local 
sidereal time at 0 or llte sidereal liour angle of the meridian of 
the Interferometer baseline al 0 from tlie angle ECM to the 
angle F.'CM' where Band E' refer to the subecjulnox point on 
the eartli’s equator before and after polar motion respectively. 

There are two systems of reckoning In use today to describe 
polar motion, One is most widely used by astronomers and 
corresponds to the usage of the Bill (Bureau Inlornational dc 
I’Ucure) while the otiicr is most widely used by geopliyslcists 
and corresponds to the usage in this document. The geo* 
physicist orients the surface of the eartli so tliat the positive 
normal points toward the zenith. The use of n rlghtdmnded 
coordinate system tlicn requires lltat the location of tlio 
“instantaneous** Eulerlan axis in the bodyTixed frame be 
specified by , the angular displacement j ^the Eulerian pole 
parallel to the prime (Greenwich) mcri^a.i, and nij, tlio 
angular displacement of tlic Eulerlan pole parallel to tlie 90® E 
meridian. This is illustrated In Figure lj»2, The net anguiar 


Fiflur* li-2. Th« gMHiMrtrical dtfinWon of tho orHiloo, ntf m, uotd by 
goopbyolcitta to IndioHo tho tocotton of tho Eulorion 
polo or ipin oxlo of tho oofth nrioii VO to tho CIO (Convofv 
tiohol Intornotional Orlsin). Tho axlo of tho CIO lo con* 
•ktorod to bo cdncMont with tho figuro mIo of tho 
oorth. 
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displacement of the Eulerlan pole Irom the figure axil is given 
by angle fit where 

d* a //I* f /,!* + 0 , (Ihl) 

The astronomer adopts a rlglitdinnded coordinate system 
on tlie celestial sphere, wllh Ihe “surface” of the celestial 
oriented so that the positive normal points toward the earth. 
As a consequence of this the astronomer Is required for 
consistency to orlem the surface of the earth with the }>osUlve 
normal pointing toward the geocentcr. The use of a right- 
hnnded coordinate system then requires that tlie location of 
(lie “instantaneous” Eulertaii axis in the body-fixed frame be 
specified by .v, the angular displacement of the Ihilerian pole 
panillcl to tlie prime (Grcenwlcli) meridian, and tlie angular 
displacement of tlie Eulerian pole parallel to the 90“W 
meridian, 


Both systems of reckoning use tlie CIO or Conventional 
International Origin as the reference pole defining tlie figure 
axis or Tlie CIO is llie nominal point of 90®H latitude am! 
is defined implieitly by the assigned nominal coordinates of a 
number of observatories around tlie world. The CIO was 
defined initially by the mean latitudes and iongiuidesXj 1 ■ 
1, 2,,.., 5 of 3 observatories as a result of latitude and 
longitude observations made at llicsc sites extending over the 
Interval 1900.0 - 1905.0. The CIO defines the basis 
vector, the axis of wltlcli makes an angle of 90“ - ?, with tlie 
verticals of each of llic S defining observatories. 

It can be sliown that for an observatory (interferometer 
baseline) wltli nominal (referred to tho CIO and the Grcenwlcli 
meridian) geocentric coordinates given by Inlltude and cast 
longitude tlie increments to tlie latitude (OCM' - 
OCM in Figure ll-I) and to the longlUide (E'CM' « ECM 
in Figure ll-I) accompanying a displacement ffif radians of 
the Eulerlan pole arc given by 


A0q ** cos Xq + sin X^ (11-2) 

AXjj * tan 0J, (Ml, sin X„ « mIj cos X„) . (11-3) 

UTO is a “raw** measure of Universal Time based on data 
obtained al a parlleular site where tlio assumptions arc made 
llirougliout the data analysis tliat 

(i) The observatory (interferometer) coordinates arc Its 
nominal coordinates referenced to tlie CIO and tlie 
Grccnwicii meridian. 
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(2) The earth has been spinning about an axis throu^i the 
CIO. 

irn Is a "true” measure or Univer*.»l Time intended to 
Hive the hour angle of the mean sun at Greenwich when the 
effects of polar motion described above are duly taken Into 
account. UTI cannot be measured from a single point on the 
earth’s surface since it Is Impossible to distinguish the effects 
of polar motion which are to be Incorporated into UTI from 
the effects of variations In the earth’s rotation (spin) rate 
which are normally incorporated into UT2. 

UTO and UTI are related by 

UTO * UTI ,j,n Xp Wj cos X^) . 

Ul-4) 

UTC is an atomic lime scale broadcast by national time 
services and maintained continuously by onsite atomic frc« 
quency standards at observatory sites. At each of m observa* 
lory sites it Is possible to observe directly the qmmUtles 


III. Dynamict of Rotating Bodlaa 

before proceeding to a treatment of the rotational 
dynamics of the earth we shall review the dynamical theory of 
rotating bodies generally, arriving finally at the UouvUle and 
Kuler equations. 

A. Rotational Dynamlca of an Aaaamblaga of 
Particiaa 

We consider an assemblage of// particles whose masses and 
positions relative to our origin of coordinates are given by m, 
and r,, / » 1 , 2, 3 • * * //, rcsj^ectively. If Is llie acceleration of 
llie fth particle rclallvc to the origin of coordinates and if S is 
tlte acceleration of tlte origin of coordinates relative to inertial 
space, then Newton’s second law applied to the Uh particle 
gives 


UTO, - UTC 




Of 


/ * 1,2,3 >• *fn 
/ - 1,2.3 


wlierc is the "insiimtaneous" lalltudcofthe/lh observatory 
and (pQi the nominal lalluidc of the /ih observatory. These two 
observfibles are related to m, iiu and UTI UTC by 

UT0| - UTC » UTI +— ^tan <Ao/ \i " 

-UTC (ll-S) 


N 


m 




/o| 
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(lll-l) 


where F, Is the net force on the /th particle due to all 
lnnucncc.s external to the assemblage and R^, is the force on 
the/lh particle due to the/th particle of the assemblage. 

Summing over all particles 




and denoting 


N AT Af /V 


f“l 


lo 1 /a' ) 

l>fil 


1t 


’*>1 * m, cosXq, + /«j sinX^, (Ih6) 

where Xq/ is tlie nominal longitude of the /(h observatory, 

Prom a large set (of the order of 50) of such observations 
llie Bill solves for the quantities ?»( f/fj and UT! ~ UTC by 
least squares adjustment. Tills adjustment procedure produces 
a value of UTI wliich is not reduced to the prime (Greenwich) 
meridian passing through tlic meridian circle of Uic transit 
telescope at Grcenwieii but to a meridian which is displaced 
from Greenwich by several milliseconds of time. The meridian 
lo which UTl is adju.sted Is called tlic meridian of lire “mean 
observatory" or the Greenwich Mean Astronomic Meridian. 
The corrections for tlic displacement of the mean observatory 
from Greenwich (usually of tlie order of 2-5 msec) are 
publisiied in the Bullciin Uorairc of the Bill. 


N 

R-ER, 

mi 

as the total external force on the assemblage and denoting 

w . E 

mi 

as the total mass of the assemblage gives 

N N N 

A(S+E "l,S ’ F+ EE R,/' ("'•2) 

M f 1 / I 
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imrociucing p the position of the center ofntasft 

N 

E 

p .-ill— (IIW) 

E «■, 

(3 1 


allows Btiuatlon (IU*2) to be written as 


N N 

«(S+» - F+2 E *;r 

/Bj ja\ 


Assuming nonrelatlvlsttc mechanics for wlrlclt 
several terms in Equation (IH'5) can be rewritten: 

(1) f; m,r,Xir,.^(2 

/«! \/»t I 

since 0 

N 

(2) 2wi^r^XS-AfpXS, 

/«t 

using Equation (Ul>3)» and if we deHne: 

(3) L the total angular momentum of the assemblage about 
the origin of coordinates 0 a$ 


In Its “weak" form Newton’s third law merely asserts the 
equality of action and reaction and gives 


N 

L *2 


i»i 


K 


7/ 


(4) N the not external torque on the assemblage about the 
origin of coordinates 0 as 


Simply stated this Implies that the force exerted by t'^/ih 
particle on tlic /th particle is equal and opposite to the force 
exerted by the /th particle on thc/th particle. If the weak form 
of Nesvton’s third law is assumed, then In the double sum 


N 


N 

E 

la I 


,,XF, 


then Equation (I II '5) can be written 


N N 

E E «/, 


t-\ t^\ 


N N 

Mp X S + 1 * N rEE'/XK,/- 

/0| /o| 

y<e/ 


(111-6) 


the terms cancel In pairs with the conclusion that the not force 
on the assemblage of particles duo to internal action and 
reaction pairs vanishes. 

Returning now to Equation (IIM) and taking the cross 
product of Equation (III-l) with the vector r and .summing 
over all the particles gives 

£ m, X S + 2 m, r, X if, ■ 2 ^ ^ P/ 

/a I /a I /ai 

N N 

+ EE',xr„. 

/St jm\ 

(III.S) 


The last term, 

N N 

E E'.xRp 

/«) /a I 
M/ 

In Equation (1 1 1-6) represents the net torque on the assemblage 
resulting from Us internal actions and reuelions. In order for 
this term to vanish It is necessary to Invoke Newton’s third law 
In Its “strong” form; namely, that 

R// ** ■' “ «//‘ 


Simply stated this Implies tlial the force exerted by thc/lh 
particle on the /th particle is equal and opposite to the force 
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exerted on the /lit particle by the /th particle and that these 
forces act along the lino joining the centers of the particles. 

If Newton’s third law In Its strong form Is assumed, then 

N JV N N 

2 2 ^ ^ £ £ ®// *■/ ^ *■/ ^***'^^ 

/ 3 | /BJ /») /a I 

/^f }i‘i 

and the terms on the RHS of Equation (lll-8'l cancel In pairs 
since f/X f/ X r, and It follows that the 

dynamical equation (lll‘6) reduces to 

MpXS + L*N. (m*9) 

It should be noted that Equation (1U“7) is not valid for 
electrodynamical Lorentz forces acting between charged parti* 
cle pairs. However, the conclusion of Equation (111*9), which is 
still valid when the mass assemblage Includes charged particles, 
must be obtained by a more extensive argument than that 
presented hero. 

in order for the term A/ p X S to vanish It Is necessary to 
choose an origin of coordinates to coincide v/ith the center of 
mass, in which case p»0 and the dynamical equation (111*9) 
further reduces to 

L = N , (HHO) 

In dealing with the dynaniics of rotating bodies it is tin* 
portant to realize that the dynamical equation governing 
rotation only assumes tliis simple form of (HI-10) wlicn 
expressed in a coordinate system whose origin coincides with 
tlic center of mass of tlic body, The coordinate system need 
not be an inertial frame in that S need not vanisli to acivieve 
this simplification. Tlic simpilfication occurs because p van- 
ishes. However, wliiic the origin of the coordinate system may 
be accelerating arbitrarily with respect to inertial space, the 
coordinate system itself cannot be rotating, for nowhere in our 
dynamics tiave we allowed for tills. 

Tlic anguiar momentum of a particle of mass m reckoned in 
a nonrotating frame wilii origin 0 relative to which tiic 
particle has position r is defined 

L = r X p (HM I) 

where 

p = mr (HI-12) 


Is the linear momentum of the particle relative to the origin of 
coordinates. If wo now considci a second coordinate fratne 
sharing the same origin 0 as the Inertial frame but rotating 
relative to It with an angular velocity w (defined as usual hi 
the right-hand sense), then the velocity of the particle relative 
to an observer at rest in the rotating frame denoted Is 

related to the velocity of the particle relative to an observer at 
rest in the inertial frame by the klnematlcal rclationslvlp 

r— (111-13) 

It follows from these considerations that the angular 
momentum of the particle may be equivalently written 

E»rXm| — ^ + <uXrj. (111-14) 

In considering an assemblage of particles of masses m and 
position vectors r^ 1 » 1 , 2, 3, • • we can express the total 
anguiaf rnomcntiim of the assemblage aS 

N 

la\ 


where 

L, = + (1IM5) 


is the angular momentum of tlie /tli particle. 

Introducing V, where 

V, (HI-16) 

is llie velocity of tlie /tli particle relative to the rotathtg 
coordinate frame, it is possible to write the total angular 
momentum of tlic assemblage of particles as 

N N 

L = 5] [r,x (61 Xr,)j .(111-17) 

1 1 
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Introducing where 

P/ " 

Is the linear momentum of the Ah particle relative to the 
rotating coordinate frame, wo can rewrite Equation (UI*17) as 


L « X) ^ p/ + E 1'/ '■pi * 

r«r l”i 


Introducing h where 

N 

h * 2 r,Xp, (III.19) 

V>*i 

is the "relative angular momentum" of the assemblage of 
particles, that Is, the angular momentum of the assemblage 
relative to the rotating coordinate frame, then Equation 
(I1M8) can be written 


AT 

L = h + 2 1**/ X X fy)] < 0*'-20) 

la I 


Using the triple vector product expansion 

AX(BXC) = (A • C)B-(A • B)C 

the total angular momentum of the aaemblage of mass points 
can be written 

N 

L = h + XI ^1 K**/ ‘ ■ 

i“\ 

(m-21) 

Considering for the moment only the *th component of 
tills vector 

ft 

= ''fc + E K'/ • '/> * 

i=l 

(111-22) 


Now using the Einstein summation and range convention 
r, • <0 ■ (r,)^ 
and 


and so 

ft 

h ” ■** E 

/-I 

which - jn be written as 

N 
la I 

(111-23) 

where 6^y is the Kroneckcr delta defined by 

( 0 k 


The quantity defined by 

= f 1'? «.r W* W 

hi 

Is a second order tensor (although we liavc not proved its 
tensor character) referred to as the inertia tensor, The inertia 
tensor is clearly symmetric in the Indices k, /, and its six 
independent elements consist of the six independent second- 
order moments of the mass distribution of the assemblage of 
particles. 

With tills result we can Write the kt|i component of the 
angular momentum of the assemblage of particles as 

or in coordinate free notation 

L = T* CO + h , (Hl'26) 
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To avoid possible confusion ll should be emphasized that 
the vector w was Introduced In liquation (II M3) as the 
angular velocity of a rotating coordinate ftame and has 
ttothlns to do with the angular rotation rate (if any) of the 
mass assemblage. 

11 is assumed wHhmU loss of generality that at some instant 
during the rotation the basis vectors spanning the nomotating 
frame and tim basts vectors spanning the rotating frame 
coincide, The Instant is called the “nwment of colneldence.“ 

I'iKing attention on the momettl of coincidence irermils us 
to compare components of tlie vector etjuallon (IIM3) and 
the tensor equation (lll‘24) as well as their time derivatives as 
represented in the rotating and nomotating coordinate frames. 
A! the moment of coincidence the basis vectors of the rotating 
and nonroiaiing coordinate frames coincide and such compari- 
sons are malliematicttlly permissible. 

ll should he emphasized that the luertia tensor T'couslsts of 
the secomhorder muss moments of the massdistrihution taken 
about the coordinate axes at tin* mommi of ivinchhu'c. If 
the mass distribution remains fixed relative to the nonroiaiing 
eoordinalD frame* tlietr clearly the body Is not rotating and 
L,csO. However T» co is stiU nonzt'ro but Is cimcollcd c-saetly 
by h. If the nrass distribulion remuins fixed relative to tiro 
rotating frame tlien clearly the body is rotating and 1. / 0. 
llowcvor,T'» CO is exactly the same value as in the nomotating 
case! Tills time, however, li ® 0, 


Tlierc is In general only one circumstance under which the 
vector CO is to he identified with the “rotation rale of a body” 
and that is the case wherein a rotating coordinate frame Is 
found such that h ** 0. In this case 

t^T-co. (Hi-27) 

However, co is still the roiailon rate of the coordinate systenii 
hut the above relationship oecnrs only for a unique value of to, 
wliicli can then be defined as ‘Tlic rotation rate of the body.” 
This unique rotating coordinate system can always bo found 
for the case of rigid bodies by fi.xing tlie rotating coordlnnio 
axes relative to the rigid body itself. 


Tlds discussion Illuminates the essential kinematic nature of 
the lerniT'- co appearing in hcpiatkm (111-26). The magnitude 
of T* CO can be changed at will by a change of coordinates. It 
is tlie sum of T* co plus h which has dynamical slgnificauco 
and yields the quantity L, The termT* coonly has dynamical 
significance if the rotating coordinate system is referenced or 
'‘attached” in some way to tiie rotating body. 


B. Rotational Dynamica of EKtandad paformabla 
BodloB 

Passing from the case of an assemblage of parUcles (o a 
continuous, e.xtcndcd and deformable mass distribution It is 
necessary to inirodtice a mass density distribution function 
p(r,/) which will In general bo Hme-depomlent. The mass 
density distribution function need not bo a dlfreronllable or 
even a continuous Ametlon of position as the body may 
possess Internal density discoiUlnullles. The extended deform- 
able body ocenpies the time-dependent volume T(/) boonded 
by the llmc-det>cndenl enclosing sorlaec AXz). The summation 
over individual parUcles Is replaced by an Integration over the 
volume l'(Z)- 

Once again wo consider two coordinate frames sharing a 
common origin 0, one coordinate frame rotating and one 
coordinate frame nonrolatlng, The angular velocity of the 
rotating frame Is defined by Hie vector w. At tlie moment of 
coincidence of tlie two coordinate frames tlie velocity vector r 
relative to the nonrotaling frame is related to tlie velocity 
vector relative to tlie rotnling frame by Ihe klnomati- 

cal relationship 


L ... .1. , , y „ 

-j-tcoKr 

(m-28) 

It will be convenient to introduce v(r, t) as the velocity 
vector relative to the rotating coordinate system andu(ri z) as 
the velocity vector relative to the nonrotaling coordinate sy.s- 
tern, Thus we have 

nr, ft 

(111-20) 

M(r. f) « r 

(HI-30) 

and so at the moment of coincidence 


p(r, z) » V (r, f) + (0 X r . 

(111-31) 


I. Dynmiiics of roiaUon without the inertia tensor. The 
total angular momentum of the e.xtended deformable body is 


p(r,z)(rX v(r,r)]rfP' (111-32) 
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and if the origin of coordinates is at the center of mass of the 
body then tlic dynamical equation governing the rotation is 
given by Equation (lUdO) as 

L - N (111.33) 

where, as before, the dot indicates a time derivative taken 
with respect to the nonrotating coordinate system, 

In establishing the time derivative of the integral In 
Equation (111*32) the deformable nature of the body must be 
explicitly recognized, To do tliis we consider two instants of 
time t and t + cit and the increment dt to t)ie angular 
momentum wiiich occurs in the Interval dt, 


dL “ 


/. 


p (r, / + (U) (r X pfr, r + dt)] dV 


VU+dt) 

p(r ,0 (rXi^(r,r)]t/F. 


-/ 

Jnt) 


In evaluating the above integrals we have adopted an 
Eulcrlan viewpoint. The vector r refers to a fixed position In 
tlie nonrotaling coordinate frame. 



Figure 111*1, Tht malhcinatical conatruetlona uaad to avaluaia tha 
tima darivativa of Inlagral quanttttaa whoaa valuas da* 
pand on an Intagral lakan throughout tha voluma of a 
body whoaa ahapa ia changing with tlma, 


And so we have 


dL w 



P (r. t +r/r) [r X v (r, t + </r)] dV 


[rXP(r, 0] dV 

3 


The volumes F(r + dt) and F(r) arc related by 


Vit + dt) = F, + Fj 
F(r) « F, . F, 



{p(r, / +(//) (rX P(r, t +r//)] 


where, as shown in Figure 111*1 (after Prager, 1973), 


(1) F, is the volume common to both F(t + d/) and F(0* 

( 2 ) Fj is the volume swept out in the interval dt by tiiosc 

portions ^hc bounding surface 5(f) whoso 

velocity V (relative to the nonrotating frame) has a 
component parallel to the positive direction of the unit 
outward normal fi, S 2 refers to that portion of 5 for 
which Ti • v>0, 

(3) F 3 is the volume swept out in the interval dt by those 
portions 53 of the bounding surface 5(f) whose 
velocity V (relative to the nonrotating frame) has a 
component parallel to the negative direction of the unit 
outv/ard normal «. refers to that portion of 5 for 
wWchfj*n< 0 . 


-p(r.0(rXi>(r,f)]}t/F, 



p(r, / +r/f) [r X d(r, t +dt)] dKy 


p(r, f) (rX p(r, /)] c/Fj , 


It is clear from Figure lll-l that 

c/Fj = ndS^ • Pdt 
(/F3 = -fidSj-i'dt 
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iiml that volume Inicgrals over volume Kj can be 
replaced, to first order In infinitesimals, by surface Integrals 
over the regions respectively 


luivc fids • udt*0. It follows that Including these regions In 
the surface Integrals will not alter the value of the integrals. 
Consequently we may write 


dL » 



|rXp(r,r + r/r)] 



[p{r,t+dt) [rXKr,/ + t/r)l 


-P(r,/) trXu(r.O]}rfK, 


-p(r,0 \rXv(r,t)]}dVdt 



p (r, / + d() [r X M (r, t + <//)} fidS^ » vdt 



p(r, t) [rX u(r, 0] ndS • vdt 



p(r, t) |rX >^(r, t)] iidS^ • i>d( , 


(llh34) 


We now make tire approximation iliat since dt is an infinitesi- 
mal time increment wliose magnitude is going to be allowed to 
slirink in the limit to a vanisiiingly small quantity, we can set 


from wldch it follows that 

L * / (p(r, r + r/r) [rX P(r, tfdt)] 

Jv(t) 

-P(r, 0 [rXp(r,01}</I' 


K, « m 


and write 


dL “ 



X 

dt 


{p{r,l+dl) frX u(r, r+r/0] 


-P(r, r) [rX M(r, ()]}dVdt 



Pit, t) [rXu(r, 01 Pfr. 0 • 


ndS. 

(m-36) 


In liquation (Ul-36) we recognize the Eulcrian time 
derivative taken with respect to the nonrotating frame 



p (r, /) [r X u (r, 0} n dS^ • vdt 



Pit, 0 [rX 0 ] ?tdS^ • t>d( . 

(111-35) 


iP Ch < + [r X P (r. / ^dt)\ 

-p(r,0[rXP(r.011 (HI-37) 

where l!(r, /) given by 

S(r, r) = p (r, 0 [r X jJ(r. /)] (UI-38) 


The entire area of tlie bounding surface 5(0 is divided into 
regions wldclr 

(1) lieiiuVj. 

(2) He In 53 , 

(3) lie iteithcr In 5j nor 5j, 

It follows from the definition of and Uiat for the 
latter classification wliich lie neither in 5j or 5j we necessarily 


is an angular momentum density field defined also relative to 
the nonrotating frame. 

Consequently Equation (1 1 1 -36) can be written 




(IU-39) 
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The quantity Rv appearing In Bqwatlon (111*39) and given by Collecting terms in Equation (111*45) gives 


R M « p(r, 0 Ir X i; (r, n) V (r, t) (111*40) 


Is a sccond*order tensor, It is an angular momentum flux 
density and shall be denoted ^(r, t), 


^(r, t) * p(r, Oir X p (r, /)] P(r. t) . (111*41) 

The elements of the tensor density are 

« P (r. 0 [r X P (r, 0], Uf (r, t) 

and so Equation (111*39) can be written 

dfi 


A'O) 


ndS . 


(11142) 


We may now use Gauss’s theorem to convert the surface 
Integral in Equation (11142) into a volume integral 


II ndS <= I 


V 


where 


which then gives 


bXi 


-f 

<4 i/r,\ ' / 


(111*43) 


+ V • * 0»-44) 


V{() 


It can be shown directly that the integrand of Equation 
(III44) can be written 

.^«|£(rXp) + p(rXP) 

+ V • (pp) (r X P) + pp • 7 (r X P) . (111*45) 


*1^ + V ’ U' * 1^*1^ + 7 • (pp) 
+ p(r XP) + pp* 7 (rXp) 

and by Invoking the equation of continuity 


(rXP) 


dt 


+ 7 • (pp) 0 


(111*46) 


This reduces to 


ii 

9/ 


+ 7 * y' p(f X p) + pp • V (r XP) (111*47) 


and Equation (111*44) assumes a final form 


L - I |p (r X P) + pp • 7 (r X p) I r/F . 

A'co ^ 

(11148) 

The lUIS of the governing equation (111*32) represents the 
torque acting on the extended deformable body. The net 
torque on the body arises as a result of the contbined actions 
of a system of body forces f and surface stresses S', N can be 
expressed quite generally as 


N(/) 


/ rXfr/F+ II S^ndf, 
V(t) JJ S(t} 


dS (11149) 


where it Is understood that in general both f and S are 
timc'depeiulcnt fields. 

The complication which must be borne in mind when 
applying Equation (111*49) is that usually the vector Held f 
representing the system of body forces and the tensor field S 
representing the system of surface stresses are referenced to 
the material of the rotating body. When expressed as vectors 
and tensors in the nonrotating coordinate system the body 
forces will in general appear as the time-dependent vector field 
f(r, t) and the surface stres^s will in general appear as the 
time-dependent tensor field S(r, /), where r refers to a fixed 
position in the nonrotating frame, even if, when viewed by an 
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observer at vest relulive to the body, these fiekls have no time 
ilepeiuleace. Thus in the noniotating fmine 


m 



rXf(r./)(/F+ 11 1 it, t)* fids 

JJm 


(Hi-50) 


The boumling surface S (lefining the volume I' of “tlie 
luuly” can be drawn arbitrarily. Material density p(r, t ) 
excluded from the volume is regarded as an external medium 
nol belonging to 'The body.'* In such cases the external 
medium can In general act on "the body" by a system of 
induced body forces and surface stresses. 

2. Dynumles of rolatlnn with (he inertia tensor. The inertia 
tensor can be Introduced into the dynamics by the klnematical 
expedient of introducing an arbitrary vector w sucli that tite 
velocity field Kr, t) is kinematically decomposed into 

Kr, /) v(r, r) + to X r 


Here the inertia tensor/’ is given by 


■/ 

Jv{ 


T**/ pCr.Ol^^^rrlf/P (111.55) 

’it) 


where / is tlie imlt tensor 


/ » 


I 0 0 
0 1 0 
0 0 1 


and where 


r* » r • f 


rr « 


(111.56) 


(IU‘57) 


Substituting l•^lUulion (111.54) into the dynamical equation 
governing rotation, Equation (111.33) gives 


wlticli when substituted into Equation (111*32) gives the total 
angular momentum as 


7*' to + 7*’ tb + h*N (HESS) 


L 



p(r, /) [r X V (r, t) H' r X fw> X r)J dV 


(111.51) 


where v(r, /) is a "rcsUUial" velocity field (which may or may 
not be small depending on the clioice of to) defined as a 
function of position r In the nonroiating frame. 

Introducing h the "residual" angular momentum vector 


where the indicated lime derivatives are reekoned relailve to 
the nonroiating frame of reference. 

Once agaiir we are faced witlv the problem of establishing 
the time derivative with respect to the notirototing frame of 
quantities, namely rand li, wh.cii depend on volume integrals 
carried out in the nonroiating frame. 

We proceed in the same manner as before by considering 
two instants of time t and t + dt and the increments dli and dT 
to the relative angujar momentum and inertia tensor rospec- 
lively in tlie Interval dt. 




h « / p(r, f) (r X V (r, /)| dV (111.52) 
t'Cr) 


we may write Equation (lll-Sl) as 


■/ 


L « I p(r, /) [r X (w X r)1 t/K+ h (111-53) 


By a set of manipulations identical to those carried out in 
liquations (111*20) through (111-24) it is possible to sliow that 
Equation (111-53) may be written as 

Lsr^co + h. (111.54) 


=/ 

Jift 


dh s / p (r. r + dl) [r X V (r, r + r//)l d\^ 


' / P (r , 

I 


0 (rX V (r./)] (/r (111-59) 


d! ” / p (r, / + r//) [r * / - r r] f/K 
4‘(f+t/<) 


/ 


p(r.r) [rV“rr] (111-60) 
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can be written as 


where i^(r, 0 is an inertia flux density given by 


■/ 

•'rt 


1 


r/h “I . 7 ; + (rX v(r,r + rf/)] 

'■'(r) 

P (r,0 [rX v(r, /))}</ Fr/r 

p(r ,0 [rX v(r, 0 ] nt/f 


7/ 


(r, 0 " P (ri 0 ir^ 7” r r) n(r, t ) . (UI- 68 ) 

,^(r, /) is a third-order tensor density with elements given 
by 

, » P (r. /) (»•’ 6 , 1 ”/*, rp njt (r, () , (III-69) 

We may now use Gauss's theorem to convert tlio surface 
integrals in Equations (111-63) (111-67) into volume integrals to 

(111.61) 


(tr 


f 

•I 


{ft (r, / + r/0 “ P (r, 0 }[r / - r r] dVdt 


+ // P (r. t) [r^ 1-t r] h dS - v dt, (111-62) 
h(t) 

From lliis It follows tiiat 


i-f 


M. 

d( 


+ V3C dV 


(111-70) 


f.f }|£('.')|r«r-r,ln.v.,?j dK 


(UI-71) 


wlicre 


ii-l 4r<'l'+/ X'fil/S (111-63) 

where fi(r, t) is a residual angular momentum density given by 

/»(r,r)“p(r.r)[rXv(r.OJ (HI-64) 

and wliere 3 C is a residual angular momentum flux density 
given by 


5f(r, t) - p (r, 0 [r X v (r, r)] n(r, r). (111-65) 

The quantity ^ is a second-order tensor density whose ele- 
ments are dCfj given by 

dC,f = p (r, r) [r X V (r, r)) , (r, t). (IU- 66 ) 

It also follows from this that 



|i(r»0(^2 


i - r r) r/ F + 


I 3^'fidi 


dS 

(111-67) 


33C: 

V • 3C * — 
bxf 

(111-72) 

K) P m 

ax'" 

(111-73) 


3, Discussion, A theoretical description of tiic rotational 
dynamics of an extended, generally deformable body lias been 
presented boHi with and wUliout the Introduction of the 
inertia tensor. The velocity of the material of the body and the 
deformation of the body with respect to time relative to the 
system of coordinates is accommodated, in tlie first place by 
the introduction of .7^(r, r), the tensor density flux of absolute 
angular momentum* and in the second place by the introduc- 
tion of 3 ? (r, r) and ,/(r, t), the tensor density fluxes of 
residual angular momentum and inertia respectively. 

Althougli this development illustrates the theoretical tools 
necessary to handle problems of tltis sort, the choice of a 
nonrotating frame of reference in which to describe the 
dynamics of rotating bodies is generally a poor one. This fact 
can be illustrated by considering the case of the rotation ofj 
rigid body. Even in this simple case the tmtsor fields K, 
and -:^do not in general vanish although 3C can be made to 
vanish by an appropriate choice of cu. 
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The terni iT • u appearing In Eqinulon (lll>SB) halhi general 
only kliieniatical slgnincance, Its value can he changed at will 
by a change in the choice ofm. In particular it can be made to 
^iniali by choosing co as constant. It ia tlie combination ol' 
T' tl) plush which has true dynamical signineance. I* ffects not 
appearitig in one will appear in the other. 

Pur the case of rigid bodies and quasirigid bodies the 
rotational dynamics can be greatly simplified by transforming 
to a rotating frame of reference. Sucli a transformation is not 
accomplished by simply introducing the inertia tensor into the 
dynamics as in Pquation (IIW5), for this tensor is still defined 
by moments of the ntass distribution taken about the nom 
rotating coordinate axes. Tlie transformation to a rotating 
frame of reference is accomplished by tlic Uouvlllc equation 
in which the necessary time derivatives are taken witli 
reference to a rotating .system of coordinates, 


At this point in the analysis the eicments of tlie tensor T 
and tlie components of the vector li are by deflnilion (Brpia* 
tions and UI*H) reckoned relative to the basis vectors of 
the rotating coordinate frame and in parlieular the appeaiitnce 
of the operator “r//f/r*' rather tium the operator “*** on the 
UlS of liquation (lll»77) indicates that we are to differentiate 
tlicso quanlilies with respect to tlie basis vectors of the rotsu* 
ing frame, Tbntjs, wo arc to consider the rate at which tlie 
components of ^ and It are changing with time when projected 
onto the basis vectors of tlie rolaling coordinate frame. I)y llie 
same token ilic components of the torque N appearing on the 
RIIS of liquation (Ul'77) must also be given with respect to 
the basis vectors of tlie rotating coordinate frame. Tills aspect 
uf tlie Uouvilie equation is discussed at some iengih in Munk 
and MacDonald (1060 pp 1^14), 

Written out in full 


C. The Uouvilki and Eular Equationt 

To obtain tlie Liouville equation we begin wiili tlie 
dynamical equation governing earth rotation expressed in a 
nonrotaiing center of mass coordinate system 

L * N . (Iit‘74) 


, (Ai 

lit lit ill 


i-toX / • w+ wX h « N 


(111.78) 


For a generally deformable body we can use the previous 
inguments to show liiut In the rotating frame 


We (lien consider a rotating center of mass coordinate 
system wliose rotation rate relative to the nonroiatingcoordi. 
nates Is given by the angular velocity vector w. At tiie moment 
of coincidence we can relate the time derivative L taken witli 
respect to the nonroluling frame to lime derivative ilLlilt 
taken witiv respect to the rotating frame by the kinematicul 
relatlonsliip 

l*^ + <oXL. (111.75) 

Tlie angular momentum of live rotating body can be 
expressed as 

L^r-w + b (111-76) 

Equations (Ui-74) (111.75) (111.76) togother yield 

“(f. oj + h) + to X (T- + h) * N (111.77) 

which is the Uouvilie equation fust obtained by Liouville in 
1858. 



wlicrc r now refers to a Tixcd position in tlie rotaiing frame, 
and as before llio relative angular momentum density /i(r, r) is 
given by 

Hr 1 1) « P(t> /) [r X V (r, /)] (I I i-u I ) 

The rotating frame tensor density fluxes and arc 
by definition obtained from their nonrolating frame eqiiiva. 
lents 3f and.^by ilic replacement of t) with v(r, t). 

Thus 

i (r. 0 “ P(r, t) [r X v (r, t)] v(r, t) (1 11-82) 

(n t) * p (r, Olr^ 7'“ r r) v (r, t) . (111.82) 

IS 
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Using Gauss’s llieorem 1» Equations (Ul*79), (IU-80) gives 


(it 




In the case of a rigid body w can be chosen so that v(r, /) 
vanishes with the consequence thaU 

^-0 

3/ 


h * 0 


If to lies along a principal axis of Inertia dhen two compo- 
nents of w , W2 tC) must vanish and so to X / » to vanishes. 

If Equation (IU-86) Is expressed In the principal axes 
coordinate system It becomes 

rfw, rf<o, r/co- 

+®"3T *^~dT • ?I 

+ {N^ - W,W3 (/l-O) + (A^a *• WjtOj (U^A)] % 

(Itm) 


This is Euler’s equation for the dynamics of rigidly rotating 
bodies obtained by Euler In 1765, 


3f * 0 
'’'’rot 

^ « 0 
*^rot 

The maintenance of the above four conditions implies that 
the coordinate system remains In corotation with the rigid 
body or that the coordinate system Is "attached” to the rigidly 
rotating body, Wlicn expressed in the corotating coordinate 
system the dynamical equation governing the rotation of a 
rigid body reduces to 

* N“wX 7* CO. (m-86) 

The quantity - co X 7- co is called the gyroscopic torque and 
vanishes If the rotation axis coincides with any of the principal 
axes of inertia. To prove this we consider the inertia tensor 
expressed in the principal axes coordinate system 

[ AGO' 

G B 0 

0 0 cj 

and the rotation vector expressed in the same coordinate 
system Is 

to = 0)j Cj + COj Cj + tOj Cg 

Tlien 

to X 7 - 60 « 

+ t 0 jt 03 (4'-Q Gj 

03^oo^(B-A)ey (in-87) 


A comparison of the Llouville equation, Equation (III-78), 
valid in the rotating frame with Equation (IU-58), its counter- 
part valid in the nonrotutlng frame procedures 

/ * to + i * w + h ■ N (ill-58) 

to+7- * N - toX L (111-78) 

dt “ di at 

where we liave used 

L 7 - 60 + h 

in Equation (IU-78), 

D. Kinetic Energy of a Rotating, Extended, 
Deformable Body 

Tlie total kinetic energy of a rotating, extended, deform- 
able body Is T where 

f 

r*i l/aiS + H^t/m (111-89) 

Jv 

where tlic Integral is taken over the entire volume of the body 
under consideration. In Equation (III-89) S is the velocity of 
the origin of coordinates relative to inertial space and r is the 
velocity of the mass element dm relative to the origin of 
coordinates, 

Now 

|S + rl^ == |S|^ + lHn2S-f (111-90) 
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nnd so the total kinetlw energy can be written 

jp J i» *'1-* 


idm. 

(IU-91) 


The vector S Is common to all the mass elements In I'and can 
pass through the integrals 

r ■ J/2 IS P / dm ftn I Ir P dm + & • I rdtth 

J y Jy d 

( 111 . 92 ) 

The position of the center of mass of the body relative to 
our origin of coordinates Is denoted by p where 


and 

* m f dm. (IH.98) 

•'p 

If (he origin of our coordinate system is placed at the center of 
mass of the body then p*0 and the translational kinetic 
energy reduces to 


(HI-99) 

where S Is now tlie wlocity of the center of mass of the body 
relative to Inertial space. 

We shall not coticern ourselves further with the properties 
of tire translational kinetic energy but shall Invesllgate In some 
detail the properties of the rotational kinetic energy. 


P 




and where M Is the total mass of the body. 
It follows from EquiUion (111.92) that 



f dm ■ Afp 


(111.93) 


(111.94) 


and using this result In liquation (111.92) the expression for tiro 
total kinetic energy of the body becomes 


The rotation vector w for a rotating, deformable extended 
body with fluid regions can be unambiguously defined as the 
angular rotation vector of the mean body axes frame in which 
the relative angular momentum h vanishes. (Note that there is 
no need for to to be parallel to any of the three body-fixed 
basis vectors of the mean body axes frame.) Itdr^^Jdt Is the 
velocity of a mass element relative to nn observer fixed at 
position r in tire rotating mean body axes frame llien dr^^Jdt 
is related to r, the velocity of the mass clement relative to an 
observer fixed In inertial space, by the formula 

r**~+wXr (UI.IOO) 

and so the rotational kinetic energy can be expressed as 


fm l/2^flsP+^fS-p+l/2 / IrPrfm. 

(1U.95) 

Hquaiion (IU*95) iUustrates liow the total kinetic energy of 
a rotating extended deformable body cun be decomposed imo 
a translational kinetic energy associated witli the motion of 
the center of mass relative to inertial space denoted and 
a rotational khtetic energy associated with velocity of rotation 
of the body about Its center of mass denoted Tror Titus we 
may rewrite Equation (111-95) as 

+ ("1-96) 


where 

W * 1/2A/ISP+MS-P (1H.97) 



which becomes 


dm 


T « i/2 
rot ' 


i 
i 


(w X r) dm 


+ 1/2/ IwXrPr/m 


(lU-lOI) 


Using ih? vector identities* 

(1) A -(BXC)*B . (CXA) 

(2) fAX BP*(AXB).(AXB) 

(3) (AX B) • fCX D)«A • (BX (CX D)] 
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Equation (lIMOl) can be written 

» 1/2 r 1 V I* dm + w f (r X V ) dm 

Jy> 

’H/2w» f lrX(u>Xr)i dm (IU-102) 
A/ 


where we have made use of Equation (UI»29), 
Defining the mass element dm as 

dm » p(i,t)dV 


are respectively the space-fixed and body-fixed expressions for 
the rate of change of hinetic energy. 

Since the Inertia tensor Is symmetric It can easily be shown 
that 


d)*7* (t) U 0) ‘T’ d), (ilhlO?) 

. (111-108) 

and so Equations (IIH05) (UI-106) can be written 


and using Equations (111*52), (111*53) and (111*55), Equation 
(111*102) can be written as 




I V 1^ dm + tb • h + w li + w • T* <*> 


^r<3/ * 1/2 r I V P rfm + to • h + 1/2 cj * 7» (W 

'V 

(111*103) 


+ 1/2 CO " / ‘ 


to 


(111*109) 


It is of interest to obtain an expression for the time dcriva* 
tive of the rotational kinetic energy 7 rot of an extended 
deformable body. Since is a scalar quantity we can 
conclude that Us time derivative can be taken relative to a 
rotating or nonrotating frame of reference with Idcntlcrl 
results. If we persist with the convention of using dot “*”to 
indicate a time derivative taken with respect to the nonrotat* 
Ing frame and dfdt to Indicate a time derivative taken with 
respect to a rotating frame then 


£ik 

dt 



du 

dt 


■ h + w * 




+ 0) ■ 7 “'^+ 1/2 u> *^' 


0 


(111*1 10) 



dT. 


rot 


dt 


Adding and subtracting 1/2 w • f ‘ w and 1/2 w * dlldt ’ to 
(111*104) In Equations (11M09) (IIM 10) gives 


where 


^rot - ‘/2 


Jv 


* li + co‘ h + l/2tb* / • ti) 


+ 1/2 to • / • to + 1/2 w / ’ (I) 


(111*105) 


dt 


itti a 1/2 

dt dt 




I V 1* dm + ^ • h + (o • ^ 


+ 1/2-^ •A<o+ 1/2(0 ‘Ip’ to 


'Kot “ ^ I V 1* + tb • li + to • (7"’ to +T’ (b + h) 


i"' 


- 1/2(0 • to 


(lIMll) 


dt 






+ 1/2(0’ / ’ 


dcj 

dt 


(111*106) 


+/ 


d0 , ifh 

m*!! 


)- 1/2 wi 

/ CJ 


dt 


to 


(111-112) 
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Using HquAtions (1U*S8) (111*70) in Equations (Hi*! 11) and 
(IH*i 12) raduoes them to 


» ill I I V 1^ rfm + (h * h + w * N - 1/2 « » w 

Jv 

(111*113) 




+ w(N«wXL)“ 1/2 6)-~‘6) (IIM14) 


chosen so that v «0 and h** 0 and Equations (111*1 13) (111* 
US) reduce to 

» wN- l/2wf» w, (UM17) 


dT 

« wN- l/3w|*(0 (IIM10) 

where the vector N In Equation (I1M17) represents the com* 
ponents of the torque as seen in the space^fixed frame and the 
vector N In Equation (I1I-118) represents the components of 
the torque as seen in the hoay*fixed frame* If one of these 
vectors is constant the other is generally tlme*dependent* 


Now 


to • wX L * 0 

and so Equation (111*1 14) reduces 


For a rigidly rotating body with w chosen tlmt V » 0 and 
ii • 0 it follows that dlliH » 0 and so ilie two expressions 
reduce to 


l/2wf‘w (111*119) 



+ ^*h + w* N« l/2 w*|i* CO (IIMIS) 


The kinematical rclatlonsliip for vector time derivatives 
taken in the rotating and nonrotating frames gives 


dT. 


fOt 


dt 


« to *N 


(IIM20> 


Equation (111*120) is recognised as the familiar work 
tlieoreni of rotating rigid bodies* 


In applying Equations (111*1 13) or (Ili*U5) to the ease of a 
generally deformable rotating body it is necessary to establish 
the time derivative with respect to the spacc*fixed and body* 
fixed coordinate frames respectively of the time dependent 
integral 


CO 


dui 

dt 


+ cjX CO 


and since 



1 V I® dm 



p(r, t) I v(r, r) 1* c/F. 


(111*121) 


we have 


co X to 0 


(b*^. (IU.U6) 

The time derivative of the rotation vector is the same 
whether viewed from the rotating or nonrotating frame, 

If wo consider the case of a rigid body, then co can be 


It can be shown that in tlie case of the space-Oxcd time 
derivative 

1/2 I p(r,()l»(r.APriF- / §drd f E-flrfS 
Jy{i) •'m ''''Sid 

(111*122) 

wlicre efr, t) is a relative kinetic energy density 

t'(r,/) * l/2p(r,r)lv(r,r)l* (UM23) 
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and where E(r, () li a relillve klneiic enetsy flux deiitliy In H'» «'oH* «l>liroxlmnllim the eeilh Is an axially 
nieaiuietl In the syBcc-llxeil rtarne. lyinmeitlc rlulil hiidv routing with unifotm "tigulai velocity U 

nbuut nil flKis coinddenl with th@ axis of figure in s 
E(r»0 • l/2p(r»/) 1^ Kr,/), (Ill*124) ^'flordlmue frame which Is eorotathig with the earth about the 

axis of figure the ?#roil» order Inertia tensorT® is given by 


In the case of the oody.lked time derivative It can be shown 
that 


Pit, t) hit t) * 


/I 0 0 

7® « 0 A 0 

0 0 C 


(fV4) 


and tlic zeroth order rotation vector is given by SI® where 


* fi dS 


Cm-125) 


if » nK . 


(IV.2) 


ficnceforili when we speak of the Inertia tensor of tire earth 
we shall mean Ute inertia tensor of the earth as measured in 
where it, t) Is a relative kinetic energy Hux density H'e body.flxcd coordinate frame unless we explicitly state 
measured in the rotating body*(lxcd frame otherw. ,e. 


E,.,y(r,0 - l/2p(r,r) Iv(f.r) v(r, /) .(UM26) The Inertia tensor of the "real" earth is7 where 


Using Gauss's theorem, Equations (IU*122) (HM25) can be 
written £ 


7 * 7® +7^ +7^ +7^ + » » , 


(IV.3) 


1/21 IvPr/m 
JV(t) 


■/f 


+ E 


(111.127) 


and the Instantaneous rotation vector of the body-fixed 
coordinate frame Is 

CO * 0“ T film’ + + flm® + < . » (lV-4) 


where the perturbation terms numbered 1 , 2, 3," * appearing 
, d f , . I /oc „ „ \ in these expressions are the result of a variety of perturbing 

' * L» * 'Xw "'X seophyateal phenomena, 

(IIM28) jiigii n convenient to use the notation 


This completes the discussion on the general question of 
the rotational dynamics of extended deformable bodies. We 
will now turn our attention to the application of the theoreti- 
cal tools developed here to the question of the rotation of tlic 
earth. 


7-7®+T 


<0 " + Sim * SliTj + Sim 


where the tensorThas elements 


(IV4) 

(lV-6) 


IV. The Llouville Equation and the 

Dynam’ss of Earth Rotation 7 « r. 

With the exception of its fluid portions which include the 
oceans, atmosphere, ground water and the liquid outer core, 
the rc.st of the "solid" earth is so nearly rigid that departures 
of the actual earth from a rigid body may be incorporated into S'Ven by the sum 

the dynamical theory of the earth’s rotation by a perturbation 
scheme. 7 » ' 


'■u 

^13 



'■ei 


^33 

(IV-7) 

731 


'*33, 


7* +7^ +7® + ■ 

1 f 9 

(lV-8) 
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and where i!ie vector m has components 

m * w,?, (|V*9) 

and Is given by the sum 

m * in^ + + * > • . (IV*10) 

In general both?aml m arc time-dependent perturbations. 

The angular monicntuin of the real earth Is given by L 
where 

l«7'Oi + U, (IV-ll) 

In the above equation to refers to the total instantaneous 
angular rotation rate, including the effects of precession, 
nutation, and spin, of the body-nxed coordinate frame. In 
general to Is a time-dependent vector. The physical definition 
of the vector to is implied by the physical definition of the 
body-fixed basis vectors 2*1 ?2 ^3 fheir orientation or 
rather their rate of change of orientation in inertia! space, 

Observationally however the situation is more ccmpllcatcd, 
The measurement of to is accomplished by combining data 
from a set of observers scattered over the cartli’s surface and 
attached to the earth’s crust at various points. The solid earth 
and in particular its crust is continuously deformed by tides 
and other geophysical processes and is also the subject of 
large-scale systematic geotectonic motions, The question of 
the pltysical measurement of a unique vector co, conforming to 
its definition, and derivable from a set of terrestrial observa- 
tions from scattered positions on the earth’s surface, becomes 
somewhat problematic at the level of ultrahigli precision 
measurements. We sliall consider this problem later in this 
work and for now will proceed on the assumption that a 
unique vector co Is an observable quantity and that this 
observable co conforms to the definition offered in the context 
of this theory. 

While the body-fixed coordinate system is corotating “with 
the crust” in some uniquely definable sense the presence of 
the fluid portions of the earth, namely the oceans, atmosphere 
and fluid core, will contribute to a nonzero value for the 
vector h. Also contributing to h will be those portions of the 
solid earth which, as a result of tectonic processes, are in 
motion relative to the body-fixed coordinate frame, 

It follows that wliile h does not vanish in the body-fixed 
frame of reference the quantity |h|/|L| is very small. To first 


order in small quantities the angular momentum vector L lias 
components in the rotating coordinate frame 

A, * f2//, 

A3 * Ca [1 + W3] -f r33 a + /<3 . (IV-12) 

To appreciate the nature of the approximations being made 
by retaining only the first-order terms in our theory it is useful 
lo recall (Munk and MacDonald I960) that the total relative 
angular momentum in the zonal circulation of the earth's 
atmosphere is of the order of 10®^ gm cm^ sec"' and lhatof 
the earth’s oceans is of the order of gm cm^ scc“', 
wlicreas the angular momentum of the rotating earth is 
roughly 6 X gm cm^ sec"'. It follows that 
|lij/lL|~ 10'’8 for the atmosphere and ~10"® for the oceans. 

Furthermore It has been shown (Smylie and Mansinha 
1971a; Mansinha, Smylie, and Chapmen 1979) that the 
clianges in the earth’s products and moments of inertia r/f 
resulting from the Chilean earthquake of I960 and tire 
Alaskan eartliquakc of 1964 are of the order of 10^® gm cm^, 
which is (0 be compared to the earth’s moments of inertia C~ 
A ~ 10'*'* gm cm^. Thus even for the largest of mass 
movements in the solid earth rfjlC'^rif/A 10“®, 

The Liouville equation governing the dynamics of earth 
rotation expressed in a rotating body-fixed I^ramc of reference 
is given by Equation (111-78) as 

ciT , •> dcj , dll , w v w I 

~7r* w + y • + w X / • w + coX h « N 

dt dt dt 

(IV-13) 

Substituting Equations (lV-5), (lV-6) and (IV-l 1) into Equa- 
tion (IV-13) and retaining only terms whicli are first order in 
small quantities, we obtain the perturbed Liouville equation 
governing earth rotation which Irt component form becomes 


dm, dr.~ dh. 

N. = Ail + n -7~ + ~ 

‘ dt dt dt 

+ n2(C-/l)m2 


(IV-l 4) 

dm^ dr,, dh. 

(C-A)m^ 

+ S2^r,3 + 

(IV- 15) 
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(Ifu 

** CSi + S 2 +.»>-. (IV* 16) 

In obliilning Equations (IV*14) ^ (IV* 1 6) wo have made use 
of the fact that /l,C,Si arc consiaius and their time derivatives 
vanisli. Equations (IV*14) - (IV*16) arc expressed in tire 
rotating coordinate frame and so tiro components A' i iV^ 
of the impressed torque on the earth must be expressed in the 
rotating frame as well. 

Equations (IV*14) * (lV-16) arc “separable” in the sense 
Ural tiro quantities //I 3 rj 3 and /13 appear exclusively confined 
to the equation for Ny TIris is a consequence, in part, of the 
restriction to a first-order theory and does not occur in a 
second-order expansion of tiro Liouville equation. This means 
in effect that tiro effects of the torque can be treated 
separately from tlic effects of the torques A^i The set of 
equations (IV*14)‘* (lV-16) decouples into what is usually 
described as a pair of cquatiorrs governing polar motion or 
“wobble’’ and iirvolving only the wi, Wj porturbations to 
lire rotatioti vector to and a single equation governiirg UTl and 
Involving only lire perturbatloir to tire rotation vector < 0 . 

Multiplying Equation (lV-15) by /, where and 

adding it to Equation (IV* 14) gives lire complex wobble 
equation; 

N, + W, = .l!i |r (m, + fm,) + a yj- fr,, + fr,,) 

+ ■“-(/», + (C-/0(m2 -/w,) 

- SI’ (t-23 - /r,a) - SI (//j - //I ,) . (lV-17) 

Recognizing that 

-j (mi, +imj) 

/fj - //I, » -/(/ij +i7ij) 


7 l » /l, 't'/Zlg 


F » /y, 


die complex wobble equation becomes 


N * ASl . { [n^ (C’‘*'/I)W- SHiJ 


whlclt together with 


(IV.aO) 


(7m, dh. 

Kf, - TO -jf- 1 « -Ji (1V.JI) 


consliiuto llic governing equations for earih toialion. Witli 
some simple manipulations these equations can be wiltten as 


;H . I £;,d sim . -;1- [w- a f -f-i (SiV- siff)] 

(IV.22) 


and 


dm. I / dr... dfi.X 

dl m “ dt di I • ^ ’ 

Writlon in tl\is way Equations (lV-22) and (lV-23) appear 
explicitly as equations governing changes on tiic eartli rotation 
vector u) including bolli polar motion (wobble) and UTi.TItc 
RHS of Equations (IV-22) and (IV*23) appear as forcing 
functions in the dynamics of eartli rotation w and are often 
referred to as the gcophysica! excitation functions for polar 
motion and UTl flucluaiions. 


It is possible in principle (Q use our present knowledge of 
geophysical processes to model tlic excitation functions and 
lienee predict polar motion and UTl from these equations. 
However, our ability to do tills successfully at lliis time is 
limited by a general lack of accurate information regarding (lie 
character of the goopliysical excitation functions. 


and introducing the complex quantities 
m « m, +/MI 5 

(lV-19) 


It seems most scienUfically productive at this point to 
reverse tlic above argument and set about to accurately 
measure m, mij mij by long baseline iiitcrforometiy or oUicr 
metliods with the objective of learning more about tlic 
geophysical excitation functions. Since tlicse functions reneet 
the effects of atmospheric and oceanic circulation, external 
gravitational torques, dislocations due to earthquake faulting, 
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fluUl motions In tho coroi el^etromiigncUc coupling boUvcon 
the core imd inanUe, clumges In sea level, ehnnges in ground 
water content and other important geodynamicalphenomenu, 
tills program sliould hold great potential for scientific 
discovery, 

liquations (1V.22) and (IV.23) are rcfcired to tlie rotating 
body*flxed basis vectors and so all quantities 

ttppeating lit them innst also be expressed relative to these 
basis vectors. If are (bo components of a gravita- 

tional torque whoso magnitude and dlrecUon arc fixed in 
Inerlial space, ilien wlien expressed relative to i?, tho 
torque componeuis are varying periodically with 

period 2ir/|iol. 

Allhuugh sve liave referred, and will conllmie lo refer, lo 
the vector to as the "rotation vector of (he earili" it sliould be 
borne in mind (lini to is in fad llie rotalion vector of a 
geophysical coordiiiaio system and strictly sireakiiig Inis only 
kinenuUlcal slgtilficance, Tlic unique “rointion vector of tire 
earlli" is a vector R for wliich tlie earlli’s angular momentum 
It is expressed as 

L » T- R . ({V-2.1) 


liquation (lV-2‘)) expresses tlie variations of UTI as a function 
of lime. 

The cpumiiUes?!’ and f appearing in the etpjaiion for polar 
molion (lV-22) as well as (he quantilles and appearing 
ill the equation for UTI (lV-28) depend on volume tntograls 
defined in (lie body-flxed frame, Tlie rigorous dennlllon of 
llieir lime derivaUves lor (be case of a generally deformable 
body arc given by liqualions (lU-84) and (111*85) respectively. 
However, later In lids work we sliall examine some useful 
approximate mcliiods for calculating (liese time derivaUves 
whicli treat llic earth as a rigid body wlUi lluld portions, 


V. Rotational Dynarffica of an Axially 
Symmetric Rigid i^arth 

Altliougb tire cuiili Is in rctillty n deformable body, its 
approximation lo a rigid body is sufficienlly good that 
considemble inslglit into llie carlli*s roiiilioirul dynamics can 
bo obtained by examining soUilions to ilio dynaintcai equa- 
tions governing cailli rolailon in Uieir ■/'erolli upproximalton 
mmicly llie special case of a rigid axially .symmetric earUi. 


We can obtain a relationship between R and w by decom- 
posing llie relative angular momentum h as 

li 7‘ 6m . (lV-25) 

I'Tom Ibis we can see lliat 

R«wt‘6m. (lV-26) 


A. Eularixn (Fore* Fraa) Motion of an Axiaiiy 
Symmatric Rigid Earth 

In 1765 lluler iiivesltgaled llie dynamics of rigidly rotating 
bodies in (be absence of extcriiul torques. Such motion has 
come to be known as "Huloiian motion." We sliall investigate 
llie Euleriaii motion of tlio eartli from llie point of view of a 
body-fixed coordinate frame and a space fixed coordinate 
friune. 


but .since |6tol/lwl « i we have 
R £>? W. 


1, Eulcriaii niuUoii of (he earth in a body-fixed frame, The 

rotational dynamics of a rigid earth in the absence of any 
(IV-27) geopliysical oxeltallon is governed by 


Uqiialton (IV-23) may be wrilien 


dm 


A B! - it, f l\/ 

Cil dt I 
LJo 




(lV-28) 



r* A 

“ SiJJT * 0 
A 


(V-1) 



tV-2) 


The ialler equation allows us to inlegrale the equation for 
«ta(f) to obtain 


"ijW • 


I dt ' « •* /ij + Wj 

Jo ' J 


( 0 ) 

(lV-29) 


wliieii are obtained from Equations (lV-22) and (IV-23) by 
setting tile cxcitiition function to %ero. 

Tlicso equations can be integrated directly to give 

Wjfr) ^ constant (V-3) 


23 



ORIGINAL PAGE IS 
OF POOR QUALITY 


and 


m(/) * exp/ (— 


(V4) 


whore 7 is a uomplox consIniU of iiUcgration, 

F * i-i 0^ , 

If wo liitrodiieo the angular rotation rale where 
°r — “ 

WO SCO that liquations (V*4) and (V*5) yield 

W(0 = exp [“ o’ + / ia^t o’’)] . 

Sotting 

d, « exp-0^ 


(V.S) 


(V.6) 


(V^7) 


Vo " -0'^ 


(V.8) 


and recalling that 


m(t) - m^it) i m^(t) 


(V-9) 


we have the solution for the rotational dynamics of a rigid 
earth in tlio absence of gcopliysical excitation 

m,(0 = cos o^(t-f^) 

^ Pe " ^o) 


or 


and 


no = Pa 


mJt) - constant. 


(V-10) 


From Equations (IV*2) and (IV-6) we have 


and so one way of viewing this solution is to see that it 
corresponds to a constant angular rotation rate of magnitude 
f2(l + Mjj) about tire axis of figure combined with 
time-dependent angular rotation rates of magnitude 12w,(/) 
f 2 /« 2 (r) about tiles'! and^j axes respectively. 

A more instructive way of viewing this solution is to 
consider the rotation vector co in terms of its magnitude and 
direction in tiie body-fixed coordinate frame rather than in 
terms of its components in the body-fixed coordinate frame. 
The magnitude of tlic rotation vector is given by w» 
(w • w)'/2 where 


w * (w * S2(l +2m») 


(V-12) 


to first order in small quantities. The direction of the rotation 
vector is specified by the angles Wifr) in^{t). 

For a positive value of m see that this solution 
corresponds to the uniform circular motion of the axis of 
rotation about the axis of figure in a prograde or west to cast 
direction. The axis of rotation moves within the earth on a 
body-fixed cone whose axis coincides with the figure axis 
and whose apex angle is 2Pg, The rotation axis completes one 
revolution about the figure axis in a period 2irlOy, This 
geometry is illustrated in Figures V-l and V-2. 


•3 



(y-u) 


Figur* V>1. Hm Eulwian (tofqtM-fiM) polar motion for tho 
an axially aymnw ^ rigMl aarth. Tha Eularlan 
rotation la confinod to tha aurfaoa of a goomolric 
of apax angla 2/3, aUgnad with th# figura axia. 


A 

*2 


of 

Of 
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FlsufvV’a. A polar vl«w of Eultri«n(torqu»>frM) polar motion, Tlw 
Eularlan pola la aaan to mova at a unifonn rata In a 
prograda aanaa around tha figura axia at a conatant 
angular dlatanca^,. 

From Stacey (.1977) we liavc 
C « B.0378 X 

A « S.OilSX I0*’'’gmcm^ 


ntul so 


C-A 

aa »im r aa . <m » 

A 


- 0.0032828 


1 

304.6 


and since a,. = (C - A!A)il and Injil corresponds to an 
interval of one mean sidereal day we sec tiral tire axis of 
rotation completes one revolution about tiro axis of figure in 
304.6 sidereal days, in tlio ease of a n'sldcwilt its angular rate 
around tlic figure axis is 27r/304,6 radians per mean sidereal 
day. 

2. Eulcrian motion of lire earth in a spacc-flxcd frame. 3'lte 
analysis of tlte eartli’s rotational dynamics from the point of 
view of a space*fixcd frame in the special case of force free 
motion does not require knowledge of the coordinate trans- 
formation equaiiotts relating the body-fixed frame to the 
space-fixed frame. This is because in tire special case of force 
free motion the angular momentunt vector L provides us with 
.an invariant direction in inertial space as a consequence of the 
contservation of angular momentum. 

Now in general 

L»T*w+h (V-13) 


and in the special case of a rigid earth we can choose a set of 
rotating body-fixed coordinatc.s such that li »0 attd so we can 
write 


L * T* 0) (V-14) 

for an appropriate choice of to corresponding to “the rotation 
rate of the earth.” 

Using liquations (IV»1) and (IV-1 1) we have 


L « /j Jl///, + /I ilntj -I* CS2(1 4 Mj) i?3 

(V-15) 

which can be written 

+ (C-/i)a(i4w,)?3 . (V-16) 

Using liquations (IV-l 1) and (V-6) in (V-16) gives 

L - /lto + /ly,.(l +"»3)^3 (V-17) 

or finally, neglecting the small term In o/n.,, 

L /Ifto + o^lTg). (V-18) 

From the results of Equation (V-18) it is clear that L, /I w, 
form a closed vector triangle and ate hence coplanar. 
Since |w| > |tr,,l these vectors can be represented as shown in 
Figure V-3. 

Tire angle |3^ is the displacement of the rotation vector cj 
from tire axis of figure Fj and the angle y is the displacement 
of the rotation vector co from the axis of angular momciuum 
L. The angle y represents a motion of tire rotation axis in 
space and appears as a "nutation." Sucli a motion of the 
rotation axis in space in the absence of externally applied 
torques is called “Eulcrian nutation” since it is associated 
with the Eulcrian (force frc-c) motion of a rigid body. TIris 
Eulcrian nutation has been called "sway” by sonic authors 
to distinguish It from "forced nutation” which arises as a 
result of impressed torques, The term "free nutation” lias 
also been used to denote Eulcrian nutation. 

Obscrvatlonally Eulerian nutation and forced nutation arc 
difficult to separate. Physically, however, tiicy arc quite 
distinct as Eulcrian nutation docs not displace the angular 
momentum vector in space and forced nutation docs displace 
the angular momentum vector in space. 
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Using Bqiuition (V.18) as an expression for L we have 


/I 


(V-22) 


and 


* (to » <0 2(j,« • ?3 + 0*^3 • ?3)'/^ . (V.23) 


vSinco 

0^ 1?* • * 0^ 

•r -3 3 r 

-0,03 *2o^W3 

and by Equation (V-19) 

to • to » |to|^ « ^3 + ^3 Pg 
we can rewrite Equation (V*23) as 



(co^ + to^ lan^P^ + 20^0)3 + 


Figure V>3. The geornetrical relationship between the earth's figure 
axis tj, Eularian rotation axle to, and artgular momen- 
tum axle L for a rigid axially symmetric earth with 
moments of Inertia A, A, and C. The angles /3, and y are 
greatly exaggerated and not drawn to relative scale, 
(After M. G, Rochester, unpublished research notes.) 


which reduces to 


JL 

A 


“ ((tOg + a + to^ lan^ /3 J . 


(V-24) 


Tlic following analysis allows us to deduce a relationsliip 
between the Eulcrian polar motion described by the angle fig 
and the Eulerian nutation described by tlie angle 7. 


From the constructions shown in Figure V-3 (M.G. 
Rocliester, unpublished research notes) we cun deduce: 

(1) By Pythagoras’s theorem 


From Equations (V-21) and (V-22) we have 


-- = sec^ (fi ~ 7) 

(CO 3 + 0^)2 

and using tlic Identity s<2c^(fig - 7) = tan^(/3^. " 7) + 1 we 
obtain 


Iwl ^ [w5 + (co 3 tan 13/ (V-iO) 

(2) By the definition of cos fi^, 

|co| = 0 J 3 sec fi^ (V-20) 

(3) By the projection of L onto 

L * ^3 |i I 

-j-=-^cos03^-7). (V-21) 


tan^ (fig - 7) = 


A I 


(^3 + ^)^ 


(V-25) 


Substituilrig Equation (V-24) into Equation (V-25) gives 


tan^(/3- - 7) = 


(g>3 + Ogf + (Jl tan^^^ 

+ o,f 
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which reduces to 

t!m (Pg ” 7) * h>« Pg ■ (V'26) 

Using tlic standard trigononietrlo formula this can bo 
written as 


However, in tlie case of the “real" cartli, clastic yielding of the 
mantle lengthens tiro period of tire wobble to roughly 435 
days. For the actual earth 

7 

and has a maximum value of roughly 


tan Pg •’ tan y 

i + tan Pg tan 7 * {033 + 0 ^) 


or 


tan 7 ^ 
tan pg CO 3 



and finally 


tan 7 
tan pg 




+ sec^ 


(V-27) 


Since Oy <<033 and 7 , pg are small angles this exact 
relationship can be approximated very well by 



We sec that the Eulerian nutation in space is roughly 
times the wobble amplitude on cartli. 

For a rigid earth: 

Stt 

Og = radians per sidereal day, 

CO 3 = 27t radians per sidereal day, 


and so 


'^r - 304.6 


where 7 ^ denotes the amplitude of the Eulerian nutation on a 
rigid earth. For a maximum value of 2Pg »« 0'.'40 arc we have 

27 ^ « l'.'30X 10"^ arc , 


27 0'.'03 X 10“=* arc , 

B. Pointot G«om«trical Dtscriptlon of Eultrian 
Motion of an Axially SymnHitrIc Rigid Earth 

The famous construction of Poinsot Is a general method of 
geometrically describing Eulerian (torquo*frcc) motion of a 
rigid body without having to integrate the governing 
dynamical equations. Since the Poinsot construction provides 
a complete description of the motion and since the integration 
of the dynamical equations generally involves the use of 
elliptic integrals, the Poinsot construction is quite useful as 
well as elegant, 

The general method of application of the Poinsot construc- 
tion is given in Goldstein (1950, pp, 159-161). The approach 
adopted in this work will be that of M.G, Rochester (unpub- 
lished research notes). We shall compute the time derivatives 
of the rotation vector u, both with respect to a space-fixed 
frame in which the invariant angular momentum vector L 
provides the reference direction, and with respect to the body- 
fixed frame with reference directions provided by the body- 
fixed basis vectors e^ S'j S’j. The result we seek can then be 
obtained by appealing to Equation (1II-116) and equating 
these two time derivatives. 

We have seen in Figure V-3 that for the case of the earth 
executing Eulerian motion the rotation axis w and the angular 
momentum vector L are inclined at angular Pg and - 7 
respectively to the axis of figure and that all three vectors 
60 , L, ^3 lie in the same plane. Since the angle 7 between w 
and L is a constant of the motion, co can only be incremented 
by the motion in the direction of the unit vector | which is 
orthogonal to both co and L. Also since L is a constant of the 
motion we can calculate the space-fixed time derivative of to, 
denoted according to our convention by w, by referring to L 
as an invariant space fixed vector. 

As shown in Figure V-4, in an interval of time dt the 
increment rito to co is in the direction of ?, a unit vector 
orthogonal to both L and to 
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construction for the case of Eulerian motion on a 
rigid aKis symmetric earth, 


The increment doi has magnitude given by 
\dio\ = lwlsin7tfX 


The dynamical equations governing the Eulerian motion of 
the rotation axis in the body*nxed frame are given by 
Equations (V-1), (V-2) and (V-6) as 


(V-29) 




where d\ is the increment in the interval dt to the azimuth, 
reckoned In a spacc-llxcd frame, of the plane containing Land 
60 . Combining the direction and magnitude of dui we have 
from Equations (V-28) and (V-29) 


dt 


^ 0 


(V-36) 

(V-37) 


where 


dui = Icfwlf 
. , , L X CO 


m = m, +/W2 


By definition 


sin 7 = 


ILXcol 


|L1 Itol 

Which when substituted into Equatioir (V-30) gives 




(V-30) 


(V-31) 


iurd which when substituted into Equation (V-35) yield 
d<c 


dt 


- “ m + OS'^) . (V-38) 


(V-32) 


Using the relationsiUps X ^2 ^ ^3 “^1 >^3 ^ 

= d' 2 , between the body- fixed basis vectors we can rewrite 
Equation (V-38) as 

o,.S 2?3 X (MJ,?j + ^ 3 ^- 2 ) (V-39) 
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mid since ?3 X P 3 « 0 we can add zero to the RUS of Equation 
(V.3‘)) to ohtiiin 

X (V40) 

or flmij|y> using Equation (V*34). 

^»o/aXw> (V41) 

The general transrorination rclullng the space*nxcd time 
derivative G to the body-fixed time derivative f/G/r/r of an 
arbitrary vector C is 


G »“ + wXC 
(It 

and so for tlic rotation vector w we Imve in general 


tl) * + to X <*) 

(It 


(V-42) 


(V43) 


which, since to X to « 0, reduces to 

tito 

lit 


to * 


(V44) 


as was siiown In Equation (IIMIO). Substituting Equations 
(V-33) and (V4I) into Equation (V-44) gives 

(V-4S) 

Equating the niagnitudes on botii sides of (lie vector 
Equation (V4S) gives 

(V-46) 
(V47) 



Now from Eigure V-3 we liave 

l^'j Xtoi » |<olsin|3^ 


The result of Equation (V49) expresses algebraically the 
result of tlie Eoinsot construction. Eoiiisot(lS52) showed that 
any continuous rotation of a rigid body is geometrically 
equivalent to the rolting of a cone, Averf vvitldn llte body, on a 
cone ftxeii within space. Tlie cone fixed witliin ilie body is 
called the polliodc cone and the cone fixed within space is 
called the herpolluKle cone. Tlie instantaneous rotation axis of 
the body relative to inertial space w lies along the line of 
coniact belween the two cones. 

The geometry of this arrangement is lllusiruted In Figure 
V-5. The riK tlon is a coniinumis rotation w around the line of 
contact between the cones, during whicli the axis of roialiun 
w describes successive circuits around the cone of apex angle 
2 q in space and also around llie cone of apex angle 2<3^ in the 
mill, 



Figur* V-5. Ttw rKuHanl Pointot conilructlon for Eulorton motion 
on a rigid axia aymmatric aarlh. Tha laiga body-fixad 
cona of apax angla 2/9, oanlarad on tha figura axIa t, 
rotia without aUpping on tha amall apaea^lxad cona of 
apax angia 2y cantarad on tN angular momantum 
vactor L 


and using Equations (V-31) and (V-47) in Eqmition (V46) wc 
liave 


|iol siii’ylV “ |w| siri/3^ 

or finally 

i sin Pg 

Oj, sin y 


(V48) 


(V49) 


During each complete circuit around the cone in space, the 
axis of rotation w progresses along the cone in the cartli only 
througli a distance equal to tlie circumrerciice of the small 
space -fixed cone. Hence as tlie roialion axis successively 
returns to the same position in space at the end of each circuit 
around L It lies at successively different posUioiis Wilhin llic 
earth, with the consequence that t|ie eartli lies in different 
positions in space. 
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This aspect of tlie motion can be understood by recalling 
that Of is the angular rate In tfw body^flxed frame of the 
azimuth of tire moving axis of rotation for which % is the 
angular rate in the sfme‘ftxtii frame of tlie azimuth of the 
moving axis of rotation, liquation (V-IO) can be written 


bodies of the solar system can be shown (Woolard, 1953>) to be 
smaller than tiieir indirect effect manifested through the 
irerturbations these h 'dies produce on the positions of the sun 
and moon relative to the earth, 


sin 

sinT 


m 0 


»''cosT 


tiinT’ 


(V.50) 


and using liquation (V-27) In Equation (V.50) 


X •» cr. 


cos 0,^3 /a^ , \ 


Since all angles are small cos 1, cos 7 » U and sec^ 
Pg *0 1. Equation (V.5J) reduces to its approximate form 


* CO3 + Oj. . (V-52) 

Equation (V*52) shows that the angular rate of <0 about L is 
sliglitly larger tlian the diurnal rate, being equal to the diurnal 
rate plus the polar motion or wobble rate. 

C. Hon^Eulerian (For^) Mcitlon of in Axially 
Symmatric Rigid Earth 

In reality the motion of the earth departs from the ideal 
case of Eulerian motion for two reasons: 

(1) The real earth is subjected to rotational excitation of 
both an internal and external origin. 

(2) The real earth is not an inilnltcly rigid body but a 
deformable solid with a strength comparable to that of 
steel ~10^ dynes cm**^ and in addition possesses fluid 
portions in the form of a liquid outer core, an 
atmospliere, and oceans. 

We shall examine here tire consequences for the earth’s 
rotation of the geophysical forcing functions and leave the 
investigation of the consequences of the departure of tlie earth 
from a rigid body for a later portion of this work. 


The geophysical forcing functions can be broadly classified 
into rotational excitation of external origin and rotational 
excitation of internal origin. Rotational excitation of external 
origin would include effects such as the lunisolar gravitational 
torques, the gravitational effects of the other bodies of the 
solar system, coupling to the solar wind by fluid or electro, 
magnetic processes, meteors passing througli the atmosphere, 
meteorites striking the earth and so on, Of the effects listed 
above only the lunisolar gravitational torques can be reliably 
demonstrated to have any observable effect on the earth’s 
rotation. Even the direct gravitational effect of the other 


External rotational excitations are distinguished by the fact 
that they, aitd only they, may alter the total angular 
momentum vector U of the earth. Such processes which cause 
a change of the magnitude and orientation of t in n 
spscc*fixed frame are studied under the general rubric of the 
tiieory of the precesrion and nutation of the earth. 

Rotational excitation of internal origin would include 
effects such as variations in ocean current systents; variations 
in atmospheric wind systems; redistribution of gtound water; 
changes in sea level; oceanic, atmospheric, and solid earth 
tides; fluid motions In tlie earth’s core; electromagnetic effects 
involving the operation of the geodymuno responsible for the 
main magnetic field of the earth; long term geologic processes 
such as post.glacial rebound, erosion and sedimentation, 
geologic uplift, and continental drift; and so on. 

Iiiicrnal rotational excitations are distinguished by the fact 
that they may not alter the total angular momentum vector of 
the earth. Tliis is true even of the electromagnetic processes 
involving the geodymmm and the earth’s magnetic fleld. The 
total angular momentum of the earth L necessarily includes 
the angular momentum of all its associated fields and in 
partieular tlto angular momentum of the geomagnetic neid. 
When this is doitc, angular momenuun Is conserved on the 
earti) for all internal processes, 

To prescribe (Ite carlli’s orientation in space it is necessary 
to .specify the orientation of the body.fixcd ba.sis vectors 
E, Cj relative to the space-fixed basis vectors By 
The external rotational excitations of the earth, prlncipsitly tlie 
lunisolar gravitational torques, arc related by phy.sical theory 
to the lime derivative of the earth’s angular momentum vector 
L The earth’s angular momentum vector and its time 
derivative arc not directly observable and so in order to deduce 
observational consequences from physical theory it Is ncccs* 
sary to invoke some geophysical modoi for the earth to relate 
the angular momentum vector to some observable gcopliysical 
quantUy. If tire angular momentum vector l> can be related to 
some observable body-fixed vector wltliln tire earth, then in 
observational as well as in theoretical practice the earth’s 
orientation in space can be described by determining the 
orientation of this observable body-fixed vector with respect 
to both the sot of basis vectors E, % Eg and the set of basis 
vectors B^ By Jn actual practice the “observable” body- 
fixed vectors chosen for this role have been the earth’s figure 
axis, Eg itself, and the earth’s instantaneous rotation axis, w, 
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Current astronomical theory liolds the fimdainental refer- 
ence direction in space to be the mean celestial pole of the 
ecliptic being defmed by the mean orbital angular inotnenUim 
vector of the earth. The present theory of the precession Is 
based cm Unown astronomical gravitational torques and known 
moment.-r of inertia of nit assumed rigid earth and ns sucit 
describes the ibeoreticai se,?ular motion of tite cartit’s axis of 
figure relative to the mean pole of the ecliptic, Tltis motion 
consists of tite sitm of the secular motion of about the 
instantaneous celestial pole of the ecliptic, known as lunlsolar 
precession and due to tite gravitational torques of the sun and 
tlie moon on the earth, plus tite secular motion of the 
instantaneous cclcstiai pole of the ecliptic about the mean 
celestial pole of the ecliptic, known as planetary precession 
and due to the perturbations imposed on ;he earth’s orbital 
plane by the other planets of the solar system, Togciltcr tliese 
two motions combine to make up general precession. 

The present tlioory of nutation is also based on known 
astronomical gravitational torques and known moments of 
inertia of an assumed rigid earth and is tabulated in such a way 
that it describes flic periodic motion of the earth’s rotation 
axis CO relative to the mean pole of llie ecliptic. 

The complete motion Is ihe sum of the secular and periodic 
components and strictly speaking sliould be obtained by 
adding the secular motion of2“3 to tlio periodic motion of ^3 
or by adding the secular motion of co to the periodic motion 
of CO. However, tire secular motion of ^3 and the secular 
motion of co arc identical (Goldrelch and Toomre 1969) and 
so the complete motion is described in practice by adding tire 
secular motion of ^3 to the periodic motion of co. 

We see tliat the present tireory of precession and nutation 
togetlicr describe the orientation of co relative to the basis 
vectors B^ for a rigid earth. In order to orient tire earth 
in space it is also necessary (but insufficient) to describe tire 
orientation of co relative to e', ^3. This requires knowledge 

of tlie location of tire axis of rotation relative to body of the 
cartiv or the effects of polar motion. However, tliat this is an 
Insufficient condition to fix the orientation of £'|?j,^3 
relative to B^ B%B^ can be seen from the fact that fixing tlio 
orientation of co In the system ?3 and B\ B^ B^ still 
allows both the set <?i%2'3 and BiBjB^ to bo rotated 
arbitrarily about co, The set of basis vectors Bi Bi B^, being 
space fixed, arc assumed to be not rotating about the direction 
CO and the rotation of the set of basis vectors S', ^3 about 

the direction (o Is measured (very nearly) by UTl, 

(It should be mentioned tliat known errors in the present 
theory of precession indicate that the set of basis vectors 


BiB^B^ is rotating at a rate of roughly 1"1 arc per century, 
Tiie new theory of the precession is intended to reduce this 
error to Hie level of roughly 0?1 arc per century, which, while 
quite small by the standards of conventional astronomical 
measurements, is still roughly ) milli arc second per year and 
probably observable by long baseline Interferometry 
techniques,) 

Alihougli the title of this section, ’’NomEulcrian Motion of 
an Axially Symmetric Rigid Earth,'* clearly embraces both the 
changes In the rotation vector relative to the spsee-flxed frame 
(precession and nutation) ami the changes of the rotation 
vector relative to tlie body-fixed frame (polar motion and 
UTI), it Is the latter phenomena witli whieli tills document will 
primarily concern itself and iUs the latter phenomena whicii 
arc described by llic equations (lV-22) and (lV-23), Conse- 
quently we slinll not be concerned witli a general development 
of the theory of precession and nutation but with a general 
development of the theory of polar motion and UTl, 

The superposition of the two processes, precession and 
nutation occurring simultaneously with polar motion and UTl 
variations, can be understood only approximately by referring 
to Figure V-5 and Imagining the external lunlsolar gravitational 
torques displacing the previously space-fixed vector L around 
on iho surface of a space fixed cone of apex angle about 47° 
(twice live obliquity of the ecliptic) with a period of roughly 
26,000 years. The "space fixed" cone in Figure V-5 is now 110 
longer space fixed but follows the vector L. This description Is 
only approximate and is in error for two reasons. 

First, the above description neglects the effects of nutation 
which would be manifested by small amplitude ('^9" are), high 
frequency ("-2 days- 18,6 years) periodic departures of L 
from the surface of this .space fixed cone. 

Second, the external gravitational torques perturb the 
Eulerian (torquo-freo) motion described so elegantly by tlio 
Poinsol construction. This can be seen by llio appearance of 
tlio torque components A*, Aj A3 in the RHS of Equations 
(lV-22) and (IV-23). In particular the juiiisolai gravitational 
torques displace the instantaneous rotation vector in a circuit 
around its Eulerian position in a retrograde sense, The radius of 
this circle is about 0?02 arc and the period of tlie circuit is 
very nearly one sidereal day. Tims the simple "cone-oivcone” 
description breaks down. 

Tlie problem of the body-fixed rotational perturbations 
wliicli occur in a rigid axially symmetric earth In response to a 
prescribed forcing function is of considerable interest in geo- 
physics for its solution allows us to model a variety of geo- 
physical processes and investigate their possible role in exciting 
polar motion and UTl fluctuations in the earth. A general 
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solution to the equations govo/ning polar motion and UTl 
variations in terms of arbitrarv excitation functions and tiieir 
integrals Is eoslly obtained, if the geophysical excitation func* 
tlorts were sufficiently well known this solution could be used 
to predict tlio position of tlie rotation pole and tiie value of 
Ifri In advance. Howeve'r, suclt a program Is not practicoi on 
tlie basis of our present geopitysical knowledge. 

I > Ccneral solution to the dynamical equations governing 
polar motion and UTI fluctuations. The dynamical equations 
governing tlie motion of tlie rotation axis In a body fixed 
frame are given by Equations (IV.22) and (iV»23). Using 
Equation (V.0) these can be written as 


given by 

f (V.59) 

Equation (V»56) can be written 

^ (V*60) 

The general solution to the homogeneous (c' » 0) equation 
(V.60)is 


« r 

ift ^ 




(V.53) 


(it cit 


whore (f Is the complex wobble excitation function given by 


m(t) - nf 

to which must be added a particular integral 


MpQ) » lo^ f e'(r') (it' 
Jo 


(V42) 


io give the general solution to the inhomogeneous (o' ^ 0) 
r “ e, + j (V»55) Equation (V-dO) as 


and where Equation (IV-22) gives 


Jf w 


An 


m 

N- 


n 


dt‘ 


dJi 

dt“ 




a«) . 


(V.56) 


In Equation (V*56) the quantity dc^fdt is the UTl exctiatlon 
function given by 


m) 


* m 0 ® + fo/ 




tf(t)e^'’''' dt' 


(V-63) 


That Equation (V-63) is the most general solution to Equa* 
lion (V.60) can be verified by direct differentiation. In Equa- 
tion (V-d3) is a complex integration constant. 


do. 


dt cn 




NJt')dt'- nr. 
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(V-57) 


The general solution to Equation (V-57) we have seen in 
Equation (lV-30) is given by 

ni^CO » e^(t) + mj (V-64) 


in Equation (IV-23). 

where 


Since 777, arc dimensionless angles the excitation func- 
tions F, de^fdt, have the units of scc“* or “frequency”. 

The solution to Equation (V-56) can be obtained by the 
usual method of variation of parameters. Introducing the 
dimensionless complex excitation function where 

(V-58) 





(V-65) 


and where is an integration constant, 

2. Some idealized example.^ of polar motion excitation in 
an axially symmetric rigid earth. As an illustration of the 
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usefulness of the solutions (V*63) and (V44) we shall eonslder 
the case of a few Idealized examples of the excitation of polar 
motion (Munk and MacDonaldi 1960). 

It can be seen from the governing equations for polar 
motion (V40) tJiat the complex wobble excitation function 

» e', + f C} and the complex coordinate of the pole of the 
rotation axis 191 »• «t| + /«tj have the same “dimensions’* of 
radians. This leads naturally to the concept of a wobble 
exeitatton axis. It will prove mathematically convenient to 
define ^ m the complex coordinate of the pole of the wobble 
excitation axis where 


produces a discontinuous and jiermanent change in the earth’s 
inertia tensor and acts as an abrupt generator of the perturba* 
tions in Equation (V» 68 ) (Smylle and Mansinha, 

1971a; Mansinha, Smylle, and Chapman 1979), 


Substituting the wobble excitation function (V*69) Into the 
general solution, Equation (V*63), gives 

mo » nf " P 7Hij' - t^) d/ 

Jq 

(V.71) 


*r» 

0 » «» {? 

/ 4), (V.66) 

If the wobble excitation Is small In the sense that 15'|<< l 
then the excitation axis defining the pole of wobble excitation 
Is associated with the unit vector^ where 

$ « 0,4 +0j 4 + [l “'|(05 + 05)]4. (V.67) 

The exciiailon pole f can be expressed In terms of tlie 
gcoplryslcal entities sucli ns the components of the external 
l 0 r(|uc A^, A'j, the perturbations to the Inertia tensor r ,3 i'j 3 , 
and the relative angular momentum components A, as 


Since we are assuming an absence of wobble excitation for 
KIq we have 

» 0 

and so Equation (V*71) becomes 

mo » “ /o^ e *■ 7 / e ^ (ii 

which Integrates to give 

.. . iO (i m f ^ 

nT(r) -7[I“C *• ®). (V.72) 


J 

liTcr 


rl 




(^i 


dr.^ dh , 

Tjj + 0 /fj 


.) 

)],(V. 


The geographic coordinates Wj(r) m^(0 of the pole of 
rotation can be obtained from Equation (V-72) by setting 
7* , + i ffi ^ , to obtain 

m , (0 ^ /, “ /, cos 0 ^ (r - rp ) + 4 sin (f - 


a. Step funedon wobble excitation. Step function wobble 
excitation can be represented mathematically as 

^r)-7//(r-g (V.69) 


whore 7is a complex constant given by7»/, +/ and where 
//(/ - tfO is the Heaviside slop function defined by 


fO t<L 


//(r-rj « < 


(V.70) 


[ 1 t>t. 


Physically we might expect step function wobble excitation 
to be an approximate model for the effect of earthquakes on 
polar motion. In such a simple model the fault dislocation 


in^iO * J 2'*'^ t “ ^0) ^ ■'1 ** ' 

(V.73) 


In this solution we sec that at t t^ the rotation pole 
coordinates arc m , » 0, ^ 0. As soon as the step function 

excitation Is imposed and the excitation pole appears at the 
coordinates 


?«7 


the rotation axis begins to describe a steady progradc clmulor 
path of radius |^| about the excitation pole. This Is illustrated 
in Figure V- 6 , The angular rate of the motion of the pole of 
rotation is 0 ^. The pole completes one circuit in an interval 
27 t/o,. 
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Fiflur* V4. tlM potar motion roouHo from • slop function 
oKeHation at timo t m Tfr# rotation oxia axocutfo 
uniform progrwto dfroutor motion atrout tfra axcMation 
axia, 


b. Delia function wobble mitatlon. Impulsive or 6-Ainc* 
lion wobble oxeltation can be represented mathematically ns 


Since vve are assuming the absence of wobble for t Kt^we 
have 

nf ^ Q 

and so Equation (V»76) becomes 
m(t) ■ ° £ Sif - df (V47) 


which Integrates to give 


m(r) 


0 



(V.78) 


Since 7 has the dimensions of time and has the dimensions 
of time'" i the quantity off is dimensionless and will serve as a 
dimensionless wobble amplitude K 


?i:r)«7a(r»g (v-?4) 


K » oj (V.79) 


where T is a complex constant given by + IJ^ and has 
the dimensions of time and v/herc 8(t - /q) is the Dirac 
5-function. The dimensions of 7 follow from the requirement 
that 



6(r- t^ldt « 1 


(V-75) 


hence 6(r Iq) has the dimensions of time'** , and that !^r) be 
dimensionless, 


and the solution, Equation (V-78), can be written 


m) 




t<io 


>t 


(V-80) 


The geographic coordinates /«j(0wa(t) of the pole of rota- 
tion can be obtained from Equation (V-80) by setting S'" AT, 
+ / /fj • W * m I + / Wj , to obtain 


Physically we miglit expect impulsive wobble excitation to 
be an approximate model for the effects of short lived atmo- 
spheric Storms on polar motion largely as a result of the 
changes in fr, and /ij which ml^t accompany such events. 


m , (0 » K^ sin cos of(j - t^) 

, (V-81) 

m^iO m slno^(<- /q)*- AT, coso^(r» t^) 


Substituting this wobble excitation function (V-74) into 
the general solution, Equation (V-63) glv?s 


m(t) 


me ^ 






dt' 


(V-76) 


In this solution we see that at time r « Iq the pole of 
rotation moves discontinuously at the time of the impulse to 
the complex coordinate -/)T or to W|(r) * K^i m^it) , 
For times t>tQ the pole of rotation moves in a steady 
prograde circular path of radius 0?! » o^|7| about the excita- 
tion pole, which for times t > Iq resides at the origin since for 
times t> tQ the excitation is zero. This is illustrated in Figure 
V-7. 
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Figure V-7. The polar motion which roauita from a fi-functlon axcita* 
tion at tima f B Tha rotation axia la diaplacad diacon- 
tinuoualy at tima t s from Ita original position and 
aubaaquantly axacutaa uniform prograda circular 
motion. 

A comparison of step function wobble excitation and 6 
function wobble excitation leads to the following general 
conclusions 

(1) Step function excitation displaces the excitation pole 
but not the rotation pole at ^ = fo' 

(2) 6 function excitation displaces the rotation pole but 
not the excitation pole (except for an interval of mea- 
sure zero) at 

(3) Both excitations lead to prograde polar motion about 
the excitation pole. 

c. Harmonic wobble excitation, Harmonic wobble excita- 
tion can be represented mathematically as 

^ (f) “ <ios a (t - t^) + J'*' sin a(t- t^) (V-82) 

whore 7® 7^ are complex constants given by P = /j + / 

7^ = and where a is the frequency of the harmonic 

excitation which is arbitrary and not necessarily equal to 
the “resonance” wobble frequency of the axially symmetric 
rigid earth. 

Physically harmonic wobble excitation on the earth occurs 
as a result of the external lunisolar gravitational torques N 2 
which, when viewed in the body-fixed rotating frame, are 
harmonically varying, ft is in fact this excitation which pro- 
duces the retrograde motion of the rotation pole about its 
Eulerian position previously mentioned. 


Alternate forms of the harmonic wobble excitation func- 
tion of Equation (V-82) are 


and 


W) * +7’- (V-83) 

m * 1/1 + 1/1 

(V-84) 


where 


7* =Up-iP) J- «i(J«+/P) 


(V-85) 


and where 


tairl 


/+ 

''JMAG 


•'Rk'AL 


= tan*"* 


'IMAG 


REA Li 


(V-86) 


In the formulation of Equation (V- 84 ) the quantities |7^| 
and |7"1 are the modulii of the complex amplitudes 7^ 7“ 
respectively. IT^I represents the amplitude of a prograde rotat- 
ing excitation function T*" whose phase angle at the epoch 
t~tQ is X'*'. |7"| represents the amplitude of a retrograde 
rotating excitation function 7" whose phase angle at the epoch 
t-tQ is)C. 


It is most convenient to proceed with the solution to the 
problem of harmonic wobble excitation by choosing the form 
of the excitation function given in Equation (V-83). If we 
consider the harmonic wobble excitation to have commenced 
at time r = Tq then substituting the wobble excitation function 
(V-83) into the general solution equation (V-63), gives 


TiKO = e 




la(t‘-iQ) -ia/ 


dt' 



(V-87) 
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which integrates to give 



which becomes 


reference frame has been analyzed in detail by Woolard (1953) 
for the case of an axially symmetric rigid earth. The lunlsolar 
wobble motion can be expressed as the vector addition of the 
complex coordinate Tftp to the fiulerian motion 1ft g, where as 
usual 

ntp “ Wp, +/»/,„ 

*= Wg, +/w^2 . (V.92) 

Formulae for the lunisolar wobble motion arc obtained by 
solving Equation (lV-22) written as 


m(0 - 


dt 




An 


(V.93) 


(V.89) 

a + o„ I- 


Equation (V-89) can be rewritten as 


m(0 


= [nf + — ^ T*-] 

I a - 0 ^ 0 + 0 ^ I 

. 7+ - '0> t -4i- 7- e"'" “ ■ 


(V-90) 


where 7} ^ i 1^2 complex linisolar torque on the 
earth expressed In tlie body-fixed frame, The solution to this 
problem is given in Woolard and Clemence (1966) 

- + 0r0087 sin 0 - 0'/0062 sin (<t> - 2 JL^) 

- 0:'0029 sin (0 >• 2 4©) + * • * (V-94) 

^ + 0"0087 cos 0 - 0!'0062 cos (0-2 Z,^) 

- 0:'0029 cos (0 - 2 /.©) + ' ” (V-95) 


where 


We sec that the solution consists of two wobble compo- 
nents: one occurring at the Eulerian frequency Og and one 
occurring at the forcing frequency o. The Eulerian component 
is a prograde rotation c^>0 of the pole of rotation with 
amplitude \M\ where 


\M\ 



0 + 0 . 



(V-91) 


The wobble occurring at the forcing frequency consists of a 
prograde component of amplitude [o^/(o - o^)] |7^| and a 
retrograde component of amplitude [o^/(o + o^)] 17*1. 


The fact that harmonic excitation at frequency a can also 
excite the Eulerian wobble at frequency o^ is intimately con- 
nected with the presence of both the annual and Chandler 
frequencies in the spectrum of the earth’s wobble, 

3. The lunisolar harmonic wobble excitation in an axially 
symmetric rigid earth. The effect of the lunisolar torques on 
the position of the instantaneous rotation axis in a body-fixed 


(1) 0 is Greenwich Mean Sidereal Time (GMST). 

(2) is the mean longitude of the moon. 

(3) Lq is the mean longitude of tire sun. 

The lunisolar torques superimpose on the slow Eulerian 
motion of angular rate Og a retrograde nearly diurnal circular 
motion of the rotation axis with radius ranging from zero to 
0l'02 arc depending on the positions of the sun and the moon. 
This motion, which results from the non-vanislung complex 
torque Win Equation (lV-22), is illustrated in Figure V-8. 

4. The damping of the wobble and wobble Q. It is apparent 
from the form of the solution for harmonically forced wobble 
on the earth, Equation (V-90), that even infinitesimal 
harmonic excitation at the Eulerian frequency Og will produce 
an infinite wobble amplitude. This unphysical prediction 
results from the fact that we have considered the earth to be 
infinitely rigid and hence free of any internal dissipation. 

The real earth is not infinitely rigid. The earth’s finite 
rigidity, in addition to greatly altering the character of the 
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Flgufv V-a. TTw combliMMi tftoct o( ih« Euitrtan polar motion riii, 
and tha Hinlaolar polar motion TIta rotation ol Is 

pro^rads whUs tlw rotation of Is r atroflrada> TIm 
vactor sum of Hf, and IfL conatmitaa tha total polar 
motion m. TIm ampHtuda |m^| latImaKlapandani. 


Eulcrian response in a niaimor (o be investigated in a later 
section of this work, will also result in Internal dissipation of 
wobble energy within the cartin Since the spectrum of 
naturally occurring wobble excitation is quite complex it will 
inevitably contain some power in the inlinitcsimal frequency 
band containing the Eulcrian frequency o^. It is tlie iniernal 
dissipation within the earth which results in tlio observed Unite 
wobble umplUudo even in tlio presence of continual wobble 
excitation. 


Strictly speaking the question of the damping of the earth’s 
wobble docs not belong in a theorctica! treatment of the 
rolnlionnl dynanties of a rigid carih for its answer itceessarily 
lies in an investigatioit of ilto detailed mecinmism of the 
wobble damping and hence an investigation of the general 
rheological nature of tlic earth iircluding its fluid portions. 
However, by introducing the “specific dissipation” or Q H is 
possible to introduce di.ssipatiou into the theory without coi\» 
fronting the question of the {{etailcd mechanism responsible 
for the dissipation. 

Tire Q of an oscillating system with total energy E and 
internal dissipation rate is dofinod to be 



where the integral is taken over one complete cycle of the 
oscillation. I t can be shown (Munk and MacDonald, 1060) that 
the Q is related to the sharpness of the resonance peak of the 
oscillator by 


(V.07) 


where o,. 1 “ Ao arc the frequencies at the half power points of 
tlie resonance curve, 

The Q of the earth's wobble is generally estimated by (Ids 
inoiliod from spectral analysis of polar motion data. Estimates 
for Q based on present day data arc oidy precise enough to 
place it roughly within the bounds 30 <2 bO (Eedersen and 
Rochester, l‘)72). Improved datu which could yield a more 
precise figure for Q would bo of assistance in understanding 
the mechanism of dissipation within the earth. 

VI. Equilibrium Deformation Fields In a Real 
Deformable Earth 

A complete (hcoreticiil understuuding of the rotational 
dynamics of lire “real” earth which incorporates realistic 
modeis for the rheology of the deformable eartlr and its fiuld 
portions as well as all the forces acliirg on them is a distant 
goal for gcodynumical theory. Such an achievemout is greatly 
hampered by: 

(1) Incomplete knowledge of the properties of the carlir’s 
fluid core and Us interaction with the rest of the earth 
including its electromagnetic and viscous cITccts, 

(2) Incomplete knowledge of tire long term rheological 
properties of (he imrer core, mantle, and crust of the 
earth. 

(3) Incomplete knowledge of the current systems In the 
oceans and atmosphere and their Interaction wllli tiro 
sltcll and each ollrer, to name but a few majirr items on 
a very long list. 

Nevertheless a simple technique developed originally by 
A, E, 11. Love (H)09) allows reasonably rigorous treatnieul of 
the effects of eartli deformations for certain classes of 
deforming force fields. Three Important geophysical disturbing 
forces whtcli arc capable of deforming the cartlv and which can 
be treated adequately by Love’s icchnique are: 

(1) Contrifugal forces and their effects on earth rotnUon. 

(2) Tidal forces and theh effects on earth rotation. 

(3) Surface loads and their effects on earth rotation. 
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Sli\ce Love’s theoretical technique plays a major role in the 
incorporation of earth deformation into the theory of earth 
rotation for a wide class of deforming force fields we shall first 
review Love’s theory in some detail before proceeding to 
discuss the effects of the above three phenomena on earth 
rotation. 

A. Lov* Numbtrt ■nd Equilibrium Earth 
Daformatlona 

The use of Love numbers in geophysics is subjected to a set 
of restricting assumptions which are; 

(1) The earth is assumed to be spherically symmetric in its 
elastic parameters and overall structure. 

(2) Thu earth is assumed to be in static equilibrium with 
the system of deforming forces. Strictly speaking this 
restrictive assumption prevents the use of Love num- 
bers in anything but problems of geostatics, However, 
their use in geodynamics is justified in instances where 
the time scale of the change? in the system of 
perturbing forces is large compared to the elastic 
response time of the earth in which case the internal 
displacement field u(r) is at ail times infinitesimally 
close to being in equilibrium with the deforming forces. 
The time scale of the earth’s elastic response is of the 
order of the transit time of a seismic wave across an 
earth diameter or of the period of the gravest mode in 
the earth’s free oscillation spectrum. Both these inter- 
vals are of the order of one hour and so the use of Love 
numbers to describe the geodynamical response of the 
earth to perturbing forces whose characteristics arc 
changing significantly on time scales large compared to 
one hour is possible, 

(3) The perturbing force field is assumed to be weak 
enough that the resulting stresses are small compared to 
the strength of earth materials, wliich is typically 10^ 
dynes cm-^. In tliis case the response of the earth will 
be linearly related to the perturbing stresses. 

We begin by considering the eartli in equilibrium in its 
unperturbed state characterized by a gravitational potential 
F®(r) and a density profile p®(r), both functions only of 
radius r - |r|. We then consider the earth subjected to a 
perturbing force field f(r) which results in an internal 
deformation field u(r). When the perturbing force field f(r) is 
derivable from a potential KP(r) , 

f(r) = - V F" (r) (VM) 


the dilation V • ii(r) are proportional to the perturbing 
potential. In general for any radius r * Irl within the earth we 
can write 


»,(r) - //V) y^(.r) 
y • u(r) - F'(r) K^(r) 


(VI-2) 

(Vl-3) 


where //'(r), F'(t) arc radial functions which depend on the 
earth’s elastic properties. 

Following Love it will be convenient to define functions 
I/(r),F(r) mch that 


m « g\r)H'(r) 
Fir) “ g\r)FXr) 


(V14) 

(Vl-5) 


where g^(r) is the positive scalar magnitude of gravity at radius 
r within the undeformed earth. In other words 


Ar) 

r 


(Vl-6) 


where m°(r) is the mass contained in a sphere of radius r given 
by 


m®(r) = 4nj 


p\r)r^dr. 


It follows that 



r(r) 

(VI-7) 


FP(r) . 

(VM) 


In general the perturbing potential V^(f) will have a solid 
spherical harmonic expansion 

■» n M 

v^(f) - ^ ('“) (cos0)(C^„^ 


i 

ii 




■i 


then it can be shown (Smylie and Mansinha, 1971) that under valid for r<a', where 0 is the geocentric colatitude and X is the 

the restrictions stated above the radial displacement M,(r) and geocentric east longitude and PJJ* (cos d) is the associated 
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Lcgondro polynomial of degree n and order m, Symbolically 
wo may write 


^"’(r) * E E TW. r<a. (VMO) 

n ni 

wlverc V^'”(r) Is the n, mlh clement in the double sum. 


Using Equations (Vl*7KVI'8) In Equation (VI-16) we sec that 


p‘(t) + l^(r) . 

r(r) ^ -* 

(VI.17) 

It follows from Equations (VI-10XVM7) tlrat p*(r) can be 
written 


It follows that 

"r«'-S7-EE "nwew (Vl-ll) 

S VJ n ffi 


-^2 S V)''rw (VM2) 

S (O n m 

where the quantities //„(0 i^m(0 s*'® radially varymg functions 
depending on the earth’s clastic properties and depending on 
the degree n of the spherical harmonic That H^,(r), 

F^(f) depend only on the degree n and not the order m of the 
spherical harmonic is a consequence of the fact that spherical 
harmonics of' the same degree but different order all Irave the 
same radial dependence. 

The deformation of the earth u(r) produces a perturbation 
in the equilibrium density profile of the earth, If p(r) is the 
density profile after deformation and p®/r) is the density 
profile before deformation, then the perturbation to the 
density p‘(r) is defined 

p‘(r) = p(r)“Po(r) (VI-13) 

and is related to the displacement Held and the original density 
field by 

p‘(r) = - u(r) ’ Vp®(r) - p®(r)V • u(r) . (VI-14) 

Now p® depends only on radius r = |r| and so 

(VMS) 

which allows Equation (VI-14) to be written 

^ W = ^ - P (OV • n(r) . (VI-1 6) 


p'(r) = 0(f) £ £ yr<f> 

n m 

where 

0(r) «(r)tp"(f)P'(r)l . 

g®(r) L J 

(VM9) 

The deformation of the earth and subsequent redistribution 
of the mass of the earth produces a perturbation to the 
equilibrium gravitational potential. If K(r) is the gravitational 
potential after deformation and I^(r) is the gravitational 
potential before deformation then the perturbation to the 
gravitational potential K^r) is defined 

F‘(r) = K(r)- F®(r) , (VI-20) 

Both K(r) and i'®(r) necessarily satisfy Poisson’s equation 
for the density distributions p(r) and p°(r) respectively, 

vV(r) = "4rrGp(r) (Vl-21) 


V^I^®(r) = -47rGp®(r) (VI-22) 

and since the Laplacian is a linear differential operator we can 
conclude from Equations (VM3) (Vl-20) (VI-21) (Vl-22) 
that 

V^K'(r) = -47rGp‘(r). (VI-23) 


It can be shown (Raula, 1968, pp, 61-69) tliat the 
gravitational potential K‘(r) resulting from a density distribu- 
tion p‘(r) given by Equation (Vl-18) will have the form 

>''(>) = E E (vi-ip) 

n m 
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where tiro quanlUlcs K„{r) are rttdhilly varying fimctlons 
dopeiullng on radful integrals Involving the density distribution 
p‘(r) and dcirending on the degree ;i of the spherical harmonic 
That depends only on the degree n and not the 
order ni of the spherical harmonic is again a consequence of 
the fact that spherical harmonics of the same degree but 
different order ait have the same radial dependence. 

\Vc now consider a perturbing potential field K^'(r) o( single 
fixed dogm Ht In other words we lake 

|/P(r) t/n(r) « ^ 

Hi 0 

It follows from liquation (VMi) that the radial displacement 
ir,.(r) within the eartli in response to tire perturbing imtential 
h>)is 

flSr) 

t^,(r) **— l^'(r) CV3.20) 

s\r) 

and that the perturbation I'*(r) to the earth’s gravitational 
potential in response to the perturbltig potential Vf^{r) Is 

h'‘(r) “ (Vl-27) 

Tlic total perturbed potential l''^(r) in tiro region 
inside the cartli is tlic sum of 

(1) F°(r) the original unperturbed potential, 

(2) the perturbing potential of the deforming force 
field causing tlic deformation ii(r), 

(3) F’C) the increment to the gravitational potential 
resulting from tlie redistribution of tlic earth’s mass 
accompanying the displacement u(r), 

F'(r) + VPlt) + I" *(r), r <u . 

(Vl-28) 

For a perturbing potential VH^v) of single fixed degree n we 
can use Equation (VI-24) to write 

|//(r) « k 0(,.) + + A;(r)] F,f(r), r <n . 

(Vl-20) 

Tlicrc ate a number of aspects to tins tlvcorctieal develop- 
ment which require empltasis at this point. 

First, the potential l'^(r) refers to a potential field from 
which the system of deforming body forces is to be derived 
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and can only be defined in the region r<a. Also the radial 
function A’„(r) depending on the elastic properties of the earth 
can only bo defined In the region r In general Equation 
(Vl-28) only has meaning in the region r «*»«. 

Second, the Newtonian gravitational potential of the earth 
F(r) (e.\clusive of the disturbing potentiai) depends only on 
the mass distribution. In the undeformod eartli we hud 

F(r) V\r) (VI-30) 

but in the deformed earth F(0 consists of the sum of 

F(r) » FV)+F‘(i-) (VI-31) 

ami so Is given in tlio region r<a by 

F(r) * F^r) + A;,(r)F,f(r). (Vl-32) 

Tliird, the analysis presented thus far is carried out from an 
Eulcrian viewpoint in which the vector r refers to some fixed 
position relative to the origin of Cviordlnates (gcocentcr) and 
involves the comparison of 1^0) before deformation with 
I’V) and F(r) after deformation (Equations (Vl-28) and 
(Vl-32)), For an observer or a particle moving witli the 
deforming earth a Uigrangian description Is more appropriate. 
A L'lgrangian description would Involve a comparison of F®(r) 
and F'[r + u(r)J and F [r + u(r)] . 

A Lagraiiglan analysis begins with the definition 

F(r) « F^OO-l F'(r) (Vl-33) 

from wliich it follows 

l^(r + u) » l'“(r + u) + F ‘ (r + u) . (Vl-34) 

A Taylor series expansion of Equation (Vl-34) gives 

F(r + u) " I'*’ (r) + V F®(r) • u(r) + F ' (r) + V F ‘ (r) • u(r) 

(Vl-35) 

wldch to first order In small quantities reduces to 

F(r + u) - l'«(r) + V F^(r) • u(r) + F‘(r) - (Vl-36) 
Now by definition 

g«(r) » - VF<^(r) 
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and since g®(r) has only a radial component owing to the Hmingmi Viewpoint 
assunwd spherical syminetry of the undeformed earth we can 
write 


II,Xr) 


Vl'O(r) • n(r) » -g«(r) • »(r) » -g®0)ir/r) (Vl.37) 
which when substltnted into Equation (Vl*36) gives 

F(r t* «) « F®(r) » n\r) n^(r) + F‘(r) . (Vh38) 

Using Equations tVH 1) and (Vl-24) In Rquatlun (VMS) gives 

1 


F(r u) * F®(r) -jc®(r) - ^ E E 

,C aH h m 


^EE A‘,,(r) Ff (r), r<«, (V1.30) 


n in 


or 


n m 


r<(t. 


Equation (V140) expresses the Ugrnnglan variation In the 
gravitational potential for an observer movitig with the 
deforming earth. The total efiectivc potential u) sensed 
by a parllelc would Include the perturbing potential F?’(r) 
responsible for the body force deformation field. Using 
■ liquation (V|.60) we can write 


F'(r + u) « F®(r) + E E u -«„w)cw. 


»l HI 

r<tt. (VMt) 


Considering once again a perturbing potential of shrglo 
harmonic degree n we can summuriM the results*. 


liukthm Viewpoint 

F(r) * F®(f) + A„(r) F/;(r) . r < « . (VI42) 


|7f(r) * F»(r) + U +X„(f)] Fj(r) , r <ff , 


(VI43) 






«,(r) » f-’ F^(r) , (VI44) 


F(r + u) * V\) + lA;,(r) - //„(r)l l^J(r) . r «<i . 

(VI45) 


I'Hr + «) “ F«(r) + (1 + A'„(r) - //„(r)] F/;(r) , /' . 

(VI46) 


love numbers of degree n, /i„ arc introduced into the 
theory by sjroclallzlng the formulae(Vl42)-(Vl40)to the case 
of an observer or particle at the surface of the earth r *«. 


Defining 


K * 


(VI47) 

(VMS) 


F(r + ii) « F“(r) + E E I'rw. 

(V140) and selling 


s“ » sV 

ia 

(VI49) 

I'Wlr ,, " 

I'(«) 

(Vl-SQ) 


l^'(«) 

(V1.51) 

W. , = 

F,^(«) 

(Vi.52) 


!',(«) 

(VI-S3) 


where it Is understood the argument n Indicates that the 
quantity is still a function of 0> X, lire geocenlric coiatiUide 
and east longitude respectively, yields the following results: 


Huterian Viewpoint 


The potentials along a spherical surhice of fixed radius 


F(fl) F^(r) + X-„F''(«) 
F^(«) * F» + (l+X*„)FP(fl) 


(Vl.54) 

(VI.5S) 


4t 
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iMgrangian Viewpoint 

Tho iwtontiak along the surface of the deformed earth 

+ 11 . 

u (a) rP(f/) (VI.56) 

S 

V{a + u) - KO((/) + ik„ - //„) K';(a) (VI.57) 

K'(« + II) - V%i) + (I + k„ - h„) V^ia) , (V1.58) 

B. MacCullagh's Formula and Parturbationa to tha 
Earth’a Inartia Tanaor 

The Love numbers allow a calculation of tho surface 
deformations of the real earth hi terms of the potential of the 
perturbing force field and hence arc of great utility In 
gcqdynamics. Another formula of equal utility Is MacCullagh's 
forinula, which when combined with Love numbers allows a 
calculation of the perturbations to the Inertia tensor of the 
real earth In terms of the potential of the perturbing force 
field. 

We begin by considering the gravitational potential l'(r) at a 
fixed point r In the region exterior to an extended spherical 
body of internal density distribution p(R). The geometry of 
the situation is illustrated In Figure VUl. 

Tho contribution <iK(r) to the gravitational potential at r 
due to the mass clement p(R)f/K at tlie position R within the 
body Is 



Figura Vl«1. The dcftnltlona of ttw goonwtrical qtianWiM UMd In 
tho dovotopmont of Iho Ihoory of tho oaith'o gravity 
fMd. 


where; 


dF(r) = (VI.59) 


where 

r * R + S (VI-60) 

The total gravitational potential at r is F(r) given by 


where the integral is to be carried out over the volume of the 
body. Tills integral can be expanded in the usual way as 


Kr) 



p(R) [-j]' P,^ {co^'ii)dV, t>R 

(VI -62) 


( 1 ) r^\tlR^\R\. 

(2) \jj is the angle between r and R . 

(3) /^,(cos i//) is the Legendre polynominal of degree n, 

It can be sliown (Mueller 1969, pp. 3-6) that with the origin 
of coordinates chosen to be at the center of mass of the object 

+ 2 /j 3 (cos 0) sin X + 2 /,j Pj (cos 0) cos X 
+^(^22 ” 

- 7,2 P\ (cos $) sin 2x] + J (VI-63) 
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where: 

(1) M is the total mass of tlie extended body. 

(2) l,j are the elements of the inertia tensor dellned in tiio 
body fixed coordinate system for wliiclr 0 is geoceittrlc 
colatitiidc and X is cast longitude. 

(3) i»JJ'(cos (?) is tire associated Ugondre judynomial of 
degree n and order m, 

liquation (Vl*63) is MacCuiiagh’s fornnda and is valid In 
tiiis form in (lie region exterior to tire body at fixed (Eulcrian) 
positions r. Tlic usefulness of MacCuliaglr’s formula for llto 
ilreory of cartii deformations is Illustrated in lire following 
example. 

Suppose tlic tmdeformed equilibrium figure of tire eartli is 
spherically symmetric witlr radius r * n, in whlclr case the 
undeformed gravitational potential Is given by 


K«(r) 


GM^ 


(VI-64) 


Since, for such an cartlr, lire undeformed inertia lensorT’ Is 
given by 


f 0 « 


/too 
0 /I 0 
0 0 4 


(VI-65) 


tlic undeformed gravitational potential may equally well be 
written 


r r 

se |(/0 +/22 ** 2 133)^3 (eos (?) 


+ 2 1 ^ /'j ^ 2 ^13 ^'2 ^ 

+|(/|j“/®)/>5(COS(?) COS2X 

- (cos (?) sin 2xJ ^r>a, (VI-66) 

since the terms Inside titc brackets vanish. 

Now if the earth is subjected to some deforming force field 
which produces displacements ii(R) internally within the earth 


jtlie inertia tensor is perturbed from the undeformed value 
r® to the deformed value Tgiven by 


4 +r 


ri 


'21 


^31 


'*12 

^13 

A ’»• 

^23 

^32 

^‘+^33 


(V(.67) 


and tiro gravitational potential is perturbed from Its 
undeformed value E®(r) to its deformed value l'(r) wlicre 
according to Equation (Vl-20) 


|/(r) « |^(r)+ l'i(r) 


(VI-68) 


Tlic gravitational potential E(r) of the deformed earth will 
also bo given by MaeCullagli’s formula Equation (Vl-63) and 
so from Equations (Vl-63), (Vl-66), (VI-67), and (VI-68) we 
deduce tliat 


+ 2 Tjj (cos (?) sin X + 2 (“,3 P\ (cos (?) cos X 

•t =|(rjj “ ryy)P\ (cos (?) cos 2X 

-> r „ (cos 0) sin 2xj + 0 . r > n . (VI-69) 


Combining Equation (VI-24) valid for /• <« with Equation 
(Vi-69) valid for r >a m can obtain, for r “ «, 


E E - '*33^ 

Pi (cos (?) + 2 P\ (cos (?) sin X + 2 r ,3 P\ (cos (?) cos X 
1 


+Ty(rj 2 - /•j,)/, (cos (?) cos 2X 
*• P\ (cos 0) sin 2xJ + 0 


(Vl-70) 


Wlicrc E®"’ (r) is the potential from which lire defonning force 
field is obtained. 
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Of particular Interest Is the si>eclal case of a second degree 
|)erturblng potential ^ir) which has a general representation 
valid forr^rtus 

»|sO 

(Vr.71) 


2 

|/n(r) ■ j [ (cos t?) + (cos 0) cos X 
+5l Pi (cos 0) shi X + C2 Pi (cos 0) cos 2X 


+ Pi (cos 0) sin 2X ] . 


(VI-72) 


Substituting Equation (Vl*72) into Equation (Vl-70) and using 
Equations (VI47) gives 

fcj [C^ /»5 (cos 0) + C[ />5 (cos 0) cos X +5{ P\ (cos 0) sin X 
+ Cj /'j (cos d) cos 2X + S\ pI (cos 0) sin 2X] 

. . A [ (^j j + ''22 " ^'‘33) 

+ 2 r23 P\ (cos Q) sin X + 2 P\ (cos 0) cos X 
+4’(''22 “ ''ji^ ^^2 " '‘12 

(V1.73) 

Since spherical harmonic functions arc all linearly independent 
we can equate the coefficients of harmonics of the same 
degree n and order m In Equation (VI.73) to obtain the five 
equations 


2 ^'*21 ” 


a ^2 


^11 + '■22” 253 - 


a 


(VI.77) 


(Vi.78) 


These five equations arc iusufficicnl to solve for (he six 
independent values of/-//. A full solution for the perturbations 
to titc Inertia tensor requires an additional equation to 
supplement the set (Vh74)'-(Vi*78). Tlie usual technique for 
obtaining a “solution” for the is to use as a supplomentury 
equation 


*0* 


(V|,79) 


Tlds equation expresses the (assumed) property that the trace 
of lire inertia tensor T^T Is a dynamical invariant whoso 
magnltudo Is not changed by the deformation of the earth. 

Rociicstcr and Smyllc (1974) liavc criticized the use of 
Eqi^atlon (Vl.79) In geodynumlcs and point out that, while 
T^T is conserved for all earth deformation fields derivable 
from a potential wldch is expandable in solid splicrical 
liarmonics, llicro exists a whole class of cartli deformation 
fields for which this is not true and for which 7'f.T is not a 
dynamical Invariant. 

TJ»c correct supplementary equation to use to provide a 
solution to the tjf is 


''ll '^''22 ■^'•33 “ 


(VI-80) 


where 6(7',. 7) is the variation in the trace of the inertia tensor 
which occurs as a result of the deformation field. 

Tlic correct solution for the r/j Is then obtained from 
Equations (Vh74) “(VI-78) and Equation (VI-80) gives 


-^'2''' c2 
' 12 G 2 


(VI-81) 


^12 ^ G “^2 


(VI-74) 


G *'2 


(VI-82) 


2k^a^ 

r w — 

13 G 2 


(VI-75) 


. cl 
G *^2 


(VI-83) 


2k^ 

r a „---l-.Cl 

G 2 


k 

f„ = (2 c\ -ic? ) (VI.84) 
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r 


n 



“'(F 




(VI.85) 


'*«“3 






‘2“ 

V 


Ct' 


(VI.86) 


In atldlllon (o above, Jhe dyniimtcfll equations (IV*22) 
ami (IV*23) governing polar motion and IJn depend on 
r * r ,3 + frj 3 , From this analysis we see that 




(Ci+iSi). 


(VI.87) 


Tltcse results, Equations (VI*81)-(VI*86), obtained for the 
case of a spljcrlcally syinmotrlc earth vvhose imdelbrmed 
inertia tensorP is given by 


’P » 


A 

0 

0 


0 0 
A 0 
0 A 


CVI-88) 


are not altered at all witen we consider the ilattencd real earth 
whose undeformed inertia tensor/^ Is given by 


ami the gravitational potential for the deformed earth K(r) Is 
given by 


l'(r) « i'®(r)+F*(r). (VM>2) 

It can readily be seen that substituting Equation (VI*9I) into 
MacCullagh’s formula, Equation CV1*63), and using Equations 
(VI.90) and (VI.02) gives 


|/i(r) « 


a 

2t^ 




C'*iF'‘22 




)y»® (eos 0) 


•t‘ 2 »'j 3 ‘ (cos 0) sin X + 2 r ,3 (cos 0) cos X 
+4 (t'jj - 1)/^2 (cos 0) cos 2X 

(cost?) sin 2x] (VI.93) 

which is identical to Equation (Vl»69). It follows that degree 2 
perturbing potentials will, even In the ease of the flattened real 
earth, result in perturbations to tlte earth’s inertia tensor 
given by Equations (VI'81)-(Vh86). 


0 0 
A 0 
0 


(V1.8Q) 


VII. Th« Rotational Dynamica of an Axially 
Symmatrlc Dafonnabla Earth 


in the ease of the flattened earth tlte equilibrium gravita- 
tional potential l'’®(r) is given by 

[(/H-/1 - 2C)e 5 (cost?) 
r 2r E 

+ 2 7 ®3 E{(cos 0) sin X + 2 E* (cos 0) cos X 

+‘^(A-A)fil(cos0) cos2X 
-/;jE5(cosd)sin2xJ (Vl-90) 


A. Th« Efftct of Rotational Doformationa of tha Earth 
on tta Eulorian Motion 

In tlic real earth the centrifugal forces of rotailoit produce 
deformations wblcb greatly alter the character of Us Eulorian 
motion, It is these deformations wlilch are responsible for the 
famous lengthening of the Enlcrian or Chandler wobble period 
from 304.6 days predicted on tire basis of rigid earth dynamics 
to the observed period of 435 days, in addition the yielding of 
the cartli to the centrifugal forces of rotation causes an 
cnbnnccmcnl of the wobble amplitude, 

Tire centrifugal force fleld f^(r) on the rotating earth Is 
given by 


for the undeformed earth, Tire deformed inertia tensor for a 
flattened eartIrTls given by 




'*13 

''21 

^+'’22 

^23 

^31 

'*32 



(VI-91) 


E‘(r) » - wX(coXr) (VIM) 

where w is "the earth’s rotation vector” and the superscript c 
will be used to denote phenomena associated with centrifugal 
effects, 

Olcrc is an instance where, from the standpoint of strict 
mathematlcui rigor, we make an error, it is to be recalled tliat 
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w Is the rotation vector of some rotating geophysical coordl* 
nate system and from the standpoint of formal logic may or 
may not (depending on Its definition) relate to the properties 
of the rotating earth. In practice, of course, because the 
geophysical coordinate system Is denned to be nearly rigidly 
attached to the nearly rigid earth the consequences of the 
“error” are entirely negligible,) 

The centrifugal force field can be obtained from the nega* 
tive gradient of u centrifugal scalar potential field t/^'Cr). 


(VII.2) 

where 

(/^•(r)»i»i.(wV»(w*r)*] (V».3) 

tM 


and where 

w* * 0) ' w , (VIM) 


Jn terms of the body-fixed basis vectors ^3 we have 
r ** X?, Cj 

(VII-5) 

CO « 0), tj + Wj Sj + a>j ^3 

and so 

U%t) » [(io^ + + cop {x^ + z^) - 

- cop^ - cop* + 2 coj cOjXy + 2 co^ oi^xz 
+ 2 ui^co^yz] (VII-6) 


Replacing x, y, z In Equation (VH-6) by their equivalents In 
spherical polar coordinates 

X ® r sin 0 cos X 

y r sin 0 sin X 

z » r cos 0 


the centrifugal potential becomes 

U%r) * - ^ r^ (co* - coj sin* 0 cos* X - co| sin* Q sin* X 

- coJ cos* d - 2 w, coj sin* 0 sin X cos X 

- 2 CO, cOj sin 0 cos d cos X 

« 2 tOj cOj sin 0 cos d sin X) , r<a, t VI 1-7) 

Equation (VII-7) can be rewritten In terms of the associated 
Legendre polynomials jPJJ' (cosd) of second degree 0i»2) 
where 

P\ (cos d) * "I* (3 cos* d " I) 

P]^ (cos d) » - 3 sin d cos d (Vll-8) 

P\ (cos d) * 3 sin* d 
to become, for r<fir, 

v%<) • + 

-•j CO3 ) 7*® (cosd) + -^ cOjCOjjPj (cosd) sin X 
2 

+ Y oiy CO3 P^ (cos d) cos X 
+ (w* - i>)\)P\ (cos d) cos 2X 

- Y w, Wj /** (cos d) sin 2xj (VI 1-9) 

This can be written as 

* -'|*r*<o*+ ^{-j) /’^’(cosd) [Cj' coswX 
mco 

+ sbi mX] , r*ia , (VII-1 0) 
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where the coefflclcnti of the liarinonlcs are given by 
Cj « -£(wj + a5*2w5) S®» 0 


C‘ 


2fl* 





WjWj 


(Vll.ll) 

Now the earth's rotation vector to Is given by 


From Equations (VIM 3) we see that the centrifugal poten* 
tlal U%t) can be tlecomimsed Into a seciihir part U^it) due to 
the steady mean rotation of the earth and independent of/nj , 
Wj, //I3 and a time-varying perturbation £/j^(r) due to periur* 
bations In the earth rotation and depending on 

U%t) * f/;(r) + (/^(r), (VIM4) 

Selecting those portions of the coefficients wtiich are inde- 
pendent of nj| , Mij, mj In Equations (VIM3) and substituting 
them Into Equation (VIMO) to obtain the secular centrifugal 
potential we find 


CO * ntwi,Fj +(1 

and so the coefficients In the centrifugal potential become 



(1“) ^*2 (cosdij 


to* » 0* + JZ* (2wj + /«* +m| +mp 

C® - ^ a* + U* ( 2m^ - j m] - i ml + «r* ) 


which can be written as 


t/f(r) * - 4 /■*«* tl ” (cos 0)] . (VIMS) 


> A 

5’t m « ^ 

C* ■ - S2* (ml - /np 

2 

52 

which to first order In m, m^ m^ reduce to 
to* » S)* i- 2 Q^m^ 

C® “ ^ fl* fl* + -^ fl* fl*«i3 
Cj * - '1‘fl* fl*/H, 

- -I"' 

cj * 0 
si ^0 


The rotation perturbation potential U^(t) is obtained by 
sclccFi's those portions of the coefficients wlrlch depend on 
rvii quantities Wj, nij, m^ In Equations (VIM3) and substl* 

(.vtH 2 ) Equation (Vll-lO) to give 

l^p(r) »-*“•?•* fl*/«3 + ’y ff* U*mj (■”) F® (cos <J) 
-yn"S 2*m, (cost?) cos X 

I 2 

- *” <i* a* Wj (-^) Pi (cos 0) sin X . (VIM 6) 

The action of the secular centrifugal potential over geologic 
lime has given rise to the earth’s equatorial bulge and the 
observed polar flattening /where/” ‘ 298,256. The observed 
value of the flattening together with the value of the earth’s 
mean radius and mean rotation rate when combined with the 
potential U^(t) allows the computation of a "secular” Love 
number of degree 2 It can be shown that Is very 
nearly equal to the fluid Love number k^f which would 
fvn 1 3t describe the yielding of the earth to the secular centrifugal 

^ potential were the earth a perfect fluid. This is regarded as a 

demonstration of the fact that for deforming force fields 
which act over long Intervals the global rheology of the earth 
closely resembles a perfect fluid. We shall not concern our* 
selves further with the secular centrifugal potential U^(r) but 
will consider the effects on the earth’s Eulerian motion of the 
rotation perturbation potential £/®(r). 


47 


ORIGINAL PAGE IS 
OF POOR QUALITY 


We see from (lie form of liquation (Vll»16) that 

w?o 


The lime derivatives of the perturbations to the Inertia 
tensor reckoned In the rotating frame of the earth are 



■ 0 


(VI1.25) 


! 


(CJ" co< m\ + S^‘ lilii mX) (VIM 7) 

where the harmonic coefficients CJ" B^vcn by 

C\ ■ -yo* 

■ -jffUlhfij (VII48) 

C® « 0 

" 0 , 

From Equations (VI-81)-(V1*86) and Equation (V1M8) we 
see that the perturbations to the earth’s rotation described by 
the dimensionless parameters Wj, m 3 cause perturbations 
to the earth's Inertia tensor given, to first order In m,, Wj, 
m 3 by 


* 0 

(VI1.19) 

mC fg 111 

>3 3(7 '"1 

(VH*20) 

ca 3(7 "'a 

(Vir-21) 


(VIF22) 

1 .V. 2/rj(7®n® 

yo Mi.xrrTi"' 

^22 3 f ' ' 9(7 ® 

(Vn*23) 

1 

''m * 9G ^^' 3 ' 

(Vll-24) 



drf^ r/wij 

dt ■ 30’ tit 

(VII.26) 


dr^^ krj^a^Q} dm^ 

dt 30 dt 

(VI1.27) 

r/rf, 

“r// 

1 d -V. 

^ it ir 

(VIU28) 

dt 

!7 * K(t7\ 

1 *7 m-*» X( / / 1 «j=au.s:.Mw--»»«ge8M 

dt ^ 9(7 dt 

(VI 1*29) 

dt 

3 dt 9(7 dt 

(Vn*30) 


Owing to the presence of the term (2/3y^Sl^M^ In 
rotation perturbation potential t/^(r) (Equation (Vlhl7)) we 
sec that Up(r) contains a term which cannot be Incorporated 
into a splierkal harmonic expansion. The presotice of this term 
in the perturbing potential is sufficient to Insure that the 
deformation field associated with the yielding of the earth to 
perturbations % nij m 3 In its mean rotation will not preserve 
the value of T^/CRcichester and Smylle 1 974). There does not 
exist In the literature at this time to my knowledge a solution 
for the quantities SCT;?) and (d/(/f) d(T^T) required In order to 
solve explicitly for the quantities rff and (clr^ldty/. As a result 
the effect of the yielding of tlie earth to centrifugal forces and 
its effect on its Eulcrian motion Is not a completely solved 
problem in geodynamics today. “Sojutions” to the problem 
which have Ignored the terms in 6(7^7') such as lliat presented 
in Munk and MacDonald (I960, pp, 25) arc In error. 

Fortunately the effect of 8(Tp') is confined to the inertia 
tensor perturbations rf,i* 1,2, 3, and so only enters into tlic 
value of UTl [msf/)) through the term In Equation 
(1V*30). So that while the effect of rotational deformations on 
UTl is at present unknown, the effect of rotational deforma* 
tlons on polar motion, Equation (IV*22), can bo calculated 
because these effects are seen to depend only on the quantities 
7 and c/F/dt where 7« r ,3 + ^ 33 , 

We can see from Equations (VII* 19)- (VI 1*24) and Equations 
(VU*25)-(VII*30) that in (he event of polar motion on the 
earth the yielding of the earth to the perturbed centrifugal 


I 

lit 

I 

l! 


I 

n 

i. 

If 


force field can be nccoinmodatcd Into the dynamical equation Substituting Equations (VU»37) (VlI-38) into Equation (Vl|. 
governing polar motion, Equation (IV’22), by Including tlic 36) gives 


effect In the terms r/f/r/r and 7 appearing In the excitation 
function on the RHS, This Is accomplished in the following 
way. 

L,F]dm ,IC-A-F\r.^ ^ 

(VII.39) 

The Eulerlan (force free) polar motion of a rigid earth Is 
governed by the equations 

Where 


« 0, 

at A 

(Vll-31) 

P »» » 

3C 

(VII-40) 

The yielding of the earth, In response to its Eulerlan motion 
modifies the rigid body motion by generating a centrifugal 
deformation excitation function given by 

From Stacey (1977) we have 
k^ = 0,29 



(Vn-32) 

a « 6.3708 X 10® cm 
s 7.2921 X ir® rad sec"' 

G = 6.6732 X 10"® cm® gm"‘ sec"^ 


where 


C = 8.0378 X 10'*'’ gm cm^ 


7*^ = rU 

(VlI-33) 

A = 8,0115 X 10^'’ gm cm^ 

(VII-41) 


from which we can deduce 


df ^ ^^^3 . *^^2^3 

dt “ dt * dt 

(VII-34) 

F = 8.082 X 10'" gm cm^ 

(VII-42) 

are given by Equations (VII-20), (VII-21), (VII-26), (VlI-27), 
above, The equation governing Eulerlan motion in a deform- 
able earth is then 

and 

•— = 1.008 X 10"® . 

(YII-43) 

dt A 

or 

(VII-35) 

Equation (VlI-39) can be written 

l, .F\dm .C~A 1, F 

V^Aldl-' A = 

fk 


U 


fo) . 
(VII-36) 

and since 

(VlI-44) 

We see that 


C-A = 

(VIl-45) 

3C “ 

(VII-37) 

we see that with an accuracy of 1 part in 10® we can write 
Equation (VII44) as 

k^a^Sf^ dm 
dt ' 3G dt 

(VU-38) 

dm .C’^A /. F \ 

* ■' A (‘' 0 -/ 1 )"”'“’ 

(VII-46) 
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The approxhiKiiloti which toads to Bquatiuu (V1M6) depends 
on FfA being mucli smaller than FjC^A or on the condition 
that 



S 


^™<<l 


which Is of course true for the earth. 


(VlI-47) 


With tlie approximation of Equation (VIM6) we see that 
the dynamical equations governing Eulcrian polar motion on 
the deformable earth can be wi'Uten to an accuracy of 1 part 
in 1 0^ ns 


M „ j Qjff „ - . (Vn.48) 

Thus we may Introduce the approximate centrifugal deforma-- 
tion excitation function to be used instead of the exact 
centrifugai deformation excitation function 2^' (Equation 
VII-32) because of the simplifications it brings to the theory, 
and write 


s F (Vlh40) 

for the approximate equation governing Eulcrian polar motion 
on a deformable earth where 



(VIhSO) 


Two equivalent forms for the equation governing Eulcrian 
polar motion on a deformable earth are 

■f-ia„7ir=0 (V1I.50 

and 

(V11.52) 

where Equation (VII-51) was obtained from Equation (Vll-46) 
by setting 


as the angular rate of Eulcrian polar motion on a deformable 
earth and where Equation (VII-52) was obtained by using 


C^A^ 

0^ » 


(Vll-54) 


as the angular rate of Eulcrian polar motion on a rigid earth. 

We see from Equations (Vll-51) and (VU-53) that the 
principle effect of the dcformabillty of the earth is to reduce 
the angular rate of Eulcrian polar motion by an amount 


‘c^«- (VII.5S) 

% 

This effectively lengthens the period of the Eulcrian motion or 
the "Chandler wobble" to 


2it ^ 2nA 


(Vll-56) 


Tins is a period of about 439 days which is an increase by a 
factor of 1 ,44 over the pciiod predicted for an equivalent rigid 
earth. 

In order to understand more fully the effects of centrifugal 
earth deformation on Its Eulcrian motion it is instructive to 
introduce the dimensionless centrifugal defortmtion excitation 
function in its approximate form, (as opposed to its exact 
form defined, in accordance with Equation (V-59), by 

^c> ^11. (VII-57) 

along with the corresponding complex coordinate of the cen- 
trifugal deformation excitation pole defined, in accordance 
with Equations (V-66), by 

55’’" = - . (Vll-58) 

The definition of by Equation (VlhS?) allows the 
equation governing Eulerian motion in a deformable earth, 
Equation (VII-52), to be written in the convenient form 


jf-io^im + r') = 0. (VIl-59) 

It can be shown that the instantaneous axis of figure of a 
body with instantaneous moments and small products of iner- 
tia given by yl, yl, C, r, 2 , ri 3 , Tjj respectively is displaced 




I 


•» 


I 

I 


I 

I 

i 
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from (l)C ,Yg coordinato axis by unglos \k^ 

A’, X2 coordinate axes resjwetively where 

(vn.60) 

and that the complex coordinate Jl^ is then given in 

terms of tire complex quunlityT* r ,3 + by 

P' ' 

II follows from Equallons (Vll-50) (VII-54) (VII-57) lliai 

r' > (VII-62) 

and so we see that the centrifugal deformation excitation pole 
given by 

(V'I-63) 

coincides with tiie instantaneous figure axis of the rolationally 
deformed earth. 

From Equations (VU-37) (VU-40) we see that 

7== m (VU-64) 

and so 

r = (VU.65) 

and using Equation (VU-45) we have 

= 0,3073 m (VU-66) 

The result of Equation (VII-66) Jikdicates that the deform- 
ability of the earth allows the instantaneous figure axis 0'"‘ to 
partially adjust itself to the location of the instantaneous 
Eulerian rotation axis m, the amount of the adjustment being 
about 30% of the total displacement of the rotation axis from 
the mean figure axis. 

B. The Effect of the Rotational Deformationt of the 
Earth on He Non>Eulerian Motion 


accompanies the changing earth rotation vector w profoundly 
alters the Eulerian (force free) motion of tite earth. The period 
of the Eulerian (Chandler) motion is lengthened by approxi- 
mately 44% and the amplitude of the Eulerian motion is 
increased by roughly 30%. 

The deformabllity of the earth will also alter the character 
of its non-Eulcrian motion or the forced motions which result 
from a combination of external torques and internal geo- 
phywv 'll excitation. The effects of earth deformations on its 
non-Euicrlan motion can be conveniently broken down into 
two separate aspects. 

The first aspect considers the effect of earth deformations 
on the geophysically Induced polar motion or equivalently the 
effect of deformations on the changes of w relative to the 
basis vectors S'j, 

The second aspect considers the effect of earth deforma- 
tions on precession and nutation or equivalently the effect of 
deformations on tljo motion of w relative to the basis vectors 


A complete treatment of the effects of earth deformations 
on its rotation would include a third aspect; namely, the 
effects of deformations on the geophysically induced varia- 
tions in UTl. As pom ted out previously in this section (Equa- 
tions VII-19 ft) this requires a detailed solution, for the case of 
the real earth, of the variations in which accompany 
variations in This shortcoming in present geodynamical 
theory has been pointed out by Rochester and Smylie (1974) 
and could be remedied by an extension of the work of 
Manshina and Shen (1974) or of Saito (1974) hut to my 
knowledge has yet to be done. 

I . Non-Eulcriaii polar motion on a deformable earth. The 
results of the analysis of Eulerian polar motion on a deform- 
able earth, sunrmarily presented in Equations (VII-49) — 
(VlI-52), indicate that on a deformable earth Eulerian polar 
displacement m gives rise to an additional polar motion excita- 
tion » - / (Fn//1) Tfi arising solely from the deformation of 
the earth in response to the original polar displacement. Thus 
as far as polar motion is concerned the deformability of the 
earth acts as positive feedback and cnliances the motion, Tliis 
viewpoint leads naturally to a prescription for the treatment of 
non-Eulerian polar motion on a deformable earth. This pre- 
scription for tile general treatment of non-Eulerian (forced) 
polar motion on a deformable earth is described below, 


It is clear from the previous analysis that the elastic yielding First, consider the complete ensemble of identifiable geo- 
of the earth to the changing centrifugal force field which physical phenomena capable of exciting polar motion on the 
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cartit and assign an index 1 , 2, 3 , . , to oncli. Equations 
(IV-22) and (IV-23) slww that according to first-order theory 
eacli mciuHr of the ensemble of geophysical phenomena only 
contributes to polar motion and UTl variations through the 
effect It has on 

(1) Tlie external torques / « 1 , 2, 3. 

(2) The perturbations to the inertia tensor 

(3) The perturbations to the relative angular momentum//, 
/« 1,2,3. 


and given by Equation (V-63) as 


T/lJf (/) » (0) + ia^ J S^\t) e rf/'j 


where 


- n 


(VU.70) 


(VI1.71) 


and 


Furthermore, as seen from Equation (IV-22), for considera- 
tion of polar motion above it is necessary to be concerned 
only with the elements A',, /»j, r,3, 


The second step of the procedure is to calculate for eacli 
member of the ensc/nble of geophysical processes its Individual 
perturbing contributions //^ //j, 1,2, 

3, . , , m to the external torque, relative angular momentum, 
and irierlla tensor respectively, Following tl/is it is necessary to 
form the complex polar motion excitation function e^,k“ 1, 
2, 3, , , . m for each member of the ensemble of geophysical 
processes according to the formula. 


-It' s' 

grV yj. 


ia, ' 


(Vll-72) 


Letting fnjif) denote the total polar motion ncciu'i'ing on a 
rigid cartlr as a result of the ensemble of geophysical processes, 
then 


HI 


(Vll-73) 


It follosvs from Equation (VII-69) that riJ^.(t) is given by'thc 
solution to 




k = 1, 2, 3, , 


(Vll-67) 


cU 

and is given by 


dm c- A ^ i. 

^ /■— ^«777,.= 53 ^ 


k"l 


(Vil-74) 


wire re 


jyk = yvj- + ii^k 


Ir = //J -t- ihl 


4* 

'^13 ^ "^23 


k - 1, 2, 3, . . . m 
k = l, 2, 3, , . , m 
k = 1, 2, 3, . , , m 


(Vll-68) 


On a rigid earth each member of the ensemble of geophysi- 
cal processes would give rise to a component of polar motion 
rn^iO k- I, 2, 3, described by the solution to the 
equation 


lUal i£zA 

dt ~ A 


ft , fc ^ 1, 2, 3, . . . ni 


(Vll-69) 


j^m,(0)+ io,, J ^ • 

(Vn-75) 

On a deformable earth each member of the ensemble of 
gcopliysical processes would give rise to a component of polar 
motion k ~ 1, 2, 3, , >>m, However, on a deformable 
earth must be calculated by recognizing that each mem- 
ber of tl/c ensemble of complex excitation functions 5^ k = 1, 
2, 3, . . . m will be affected by the effects of “positive feed- 
back*’ resulting from the yielding of the earth, The earth 
deformation will produce, for each member /r = 1, 2, 
3 . . . m of the ensemble of complex excitation functions, an 
additional excitation k= I, 2, 3.,.m> arising solely 
from the effects of the deformation itself where, according to 
Equation (Vll-SO) 




kc 


. -fc 
- / — 7~ nf 
A 


(Vn-76) 
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On a deformable earth then, each member of (he ensentble of 
geophysical processes will have an effective complex excitation 
Amcllon given by A’“ 1 , 2, 3 . . , /n, and will give rise 

to a component of polar motion Wi^(f)k “ I » 2i 3, . . , m 
described by the solution to the cciuatlon 


tfr 



ilTif * . 


k ® l , 2, 3 •••«), 
(Vll-77) 


ami given by liquation (V’(i3) as 

iTl'^XO) + /Oj. I “i- 

Jo 

(Vll‘7B) 

where 

p'c a 1,2,3.* ■•//(. (Vlh70) 

An alternative and completely equivalent formulation of 
polar motion on a deformable earth can be obtained by 
substituting liquation (Vlb76) into Equation (Vll-77) to 
obtain 


m^t) 






“f'Vf ^ 


mi''' 

ttnurnrxtM 

tit 


. / C-A-h'\ 


iWi^ ss ® l, 2, 3, • • * m 

(Vll-80) 


which is shown by Equation (V-63) to Irave the solution 


m^it) 



k * 1, 2, 3, • • * m , 


(Vll-81) 


where we have made use of Equation (Vll-53) 


. (Vli-82) 

0 A 

Letting m(t) denote the total polar motion occurring on a 
deformable earth as a result of the ensemble of geophysical 
processes then 


IN 


mO) * 2 

fcal 



(Vll-83) 
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It follows from Equations (Vll-77) and (VH-80) llml )W(0 
Is the solution to either of the two equivalent difforcnllal 
equations 

M.. i £ ® y. (Vll-8d) 

(i{ A 


or 

<a . , (r\-4-£i am . (VIS-85) 

* \ A I f-, 

The general solullons to liquations (VII-84) and (VH-85) arc 
given by 


mf(t) ® h)7(0)*t*fo„ / y (s''* r//'l 

L Jo k-i J 


(Vn-86) 


and 


jn(t) ^ ['r??(0)'t*/aQ f y dA 

L Jo k n J 


(Vll-87) 


respectively, 


2. Non-Eiilcriait precession and nutation on a deformable 
earth, In principle the clastic yielding of the eartli to the 
system of body forces which give rise to the torques responsi- 
ble for precession and nntution alters (he observed precession 
and nutation from that whlcli would prevail wore tire earth a 
rigid body. It has been ,sbown (Lamb 1945. pp. 724 fl) that in 
the case of a disturbing force distribution which Is fixed In 
inertial space tlic precession of a rotating mass of ideal iluld 
proceeds exactly the same ns If the mass were solid Uuougli- 
out, Furtltcrmorc it has been shown (Limb 1945, pp, 724 fO 
(hat when the disturbing force distribiitlon varies slowly rela- 
tive to Inertial space with a period Injih then the precession of 
a mass of ideal fluid rotating with angular velocity w still 
proceeds almost exactly (he same as if (he ma,ss were solid 
throughout, providing tliat (he ratio of co//j is small compared 
to c where e Is the ellipticity of the rotating mass, For the 
carthe = 3X 10“^ 


a, PrecessiOit, Considering tlic case of precession on a 
deformable earth wltcre 27 t/w « I day and lit/n « 26000 
years we see that co/« « 10"'^ , wldcli is very much less than 


S3 



3 X 10"^. Wc conclmlti lhal lb effect of the ylcliling of a 
deformable earth to lire system of body forces driving the 
precession will have a negligible effect on the observed prcces* 
Sion, the earth will essentially prcccss as though it were a rigid 
body. 

At the lAU General Assembly in Grenoble in 1976 It was 
recommended that an iirrproved theory of the precession based 
on the PRS star catalogue and referred to a now standard 
epoch of 2000.0 be adopted. The dclalls of the new theory of 
earth precession have been published by Ucsko ct al, (1977). 
The corrections to the previous theory arose principally fronr 
a failure to distinguish between the “catalogue equinox" of 
tlio FK4 star catalogue (the zero point of riglil ascension on 
the catalogue celestial equator) and the “dynamical equinox" 
(the intersection at the ascending node of tlie ecliptic and tlio 
terrestrial equator) combined with new values for the masses 
of the planets of the solar system and some effects due to the 
galactic rotation of the PK4 catalogue stars. The new theory 
oniic earth’s precession is expected to be accurate to lOl'j arc 
per century or roughly ±1 .0 milliarc second per year. 

h, Nutation, Considering the case of the nutation on a 
deformable earth where 2ir/n ^ 18,6 years we sec that cj/n » 
1.4 X lO"*!, which Is sufficiently near to 3 X 10”^ to 
expect that tlic effect of the yielding of a deformable earth to 
the .system of body forces driving the nutation might produce 
measurable discrcpaiicics when precise observations of earth 
nutation arc compared against that which is tlicorctlcally pre- 
dicted for a rigid earth, Woolard (1953, pp, 1 36) was evidently 
aware of this problem and adopted an “observational” value of 
9','210 arc for the constant of nutation in preference to 
“theoretical” values of the order of 9','224 arc derivable from 
the relationships of the constant of nutation to other astro- 
nomical constants wliich hold in the case of a rigid earth. 
Recent observational determinations of the nutation constant 
place It in the range 9"201 - 9"206 arc, with the discrepancy 
between theory and observation now of the order of 0'.'02 arc 
being attributed to the effects of the deformation of the earth. 

c, Woolard ’s theory of the nutalknt/. The present lAU 
theory of the earth’s nutation is due to Woolard (1953) and 
describes the thcoroticui motion of the instantaneous rotation 
vector CO of an assumed rigid earth with an axially symmetric 
mass distribution relative to the sot of space-fixed basis vectors 
^ description of Woolard’s theory requires the 
introduction of three reference equators: the equator of angu- 
ht momentum, the equator of figure, and the equator of 
rotation. Each equator is defined as passing through the earth’s 
center of mass pcr[)cndicularly to the angular momentum axis, 
0\e figure axis, and the rotation axis respectively. It can be 


seen from Figure VU-I titat tht cquuior of figure and the 
equator of rotation Intersect each other with ufi angle where 

d, IWJ (VU-S8) 


is the amplitude of the Eulerlan polar motion, It can also be 
seen from Figure VIl-l that the equator of angular motnottlum 
aitd the equator of rotation intersect each other with an angle 
7 where, from Equatloits (V-6), (V-l 1), and (V-27), we have 

X (Vll-89) 

F|* ' ■ 

where 

d,, lwr,| 

is the untplitude of the Uulcrian wobble (of Figure V-7), 

It will also be convenient to introduce: 

(1) the ascending node of the fixed ecliptic of the 
fundamental epoch on the equator of angular 
momentum. 

(2) 7 ,. the ascending node of the fixed ecliptic of the 
fundamental epoch on the equator of rotation. 

(3) 7 y. the ascending node of the fixed ecliptic of the 
fundamental epoch on the equator of figure. 


In proceeding with his solution Woolard actually Integrated 
a set of differential equations ktiown as Poisson’s equations, 
which approximately describe the motion in inert iat space of 
the earth’s angular momentum vector L (Kinoshita ct al„ 
1979). The solutions to Poisson’s equations we shall denote as 
\pQ Oq, where Oq arc time-dependent Euler angles defined 
by Woolard referenced to a set of space-fixed basis vectors j?, , 
^ 2 , ^2 defined by the fixed ecliptic and the fixed mean 
equinox of the fundamental epoch of 1900.0. According to 
the conventions adopted by Woolard, the angle 27 T’- 
approximately represents the combined effect of lunisolar 
precession (since it is referred to the fixed ecliptic of the 
epoch) and lunisolar nutation of the angular momentum axis 
since the epoch and the angle Qq approximately represent the 
obliquity of the fixed ecliptic of the epoch on the angular 
momentum equator of date. In Woolard’s work the corre- 
sponding motion of the earth’s axis of figure, described by 
angles and the earth’s axis of rotation, described by 
angles 0^, is obtained to the order of the approxunations 
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Figure VIM, The relative location on the celeatial aphere of the equator of figure, equator of 
angular momentum, and equator of rotation, The anglea /t, and y are greatly 
exaggerated and not drawn to relative scale, 


Involved In llio theory by the nddltions of correction terms 


Ogi' SO^ 


(vn-‘)0) 


cpocli of the !ixis of figure nnd of tlic axis of rotation co 
respectively. The angles Op 0^ represent, to lire order of the 
approximations involved, the obliquity of the fixed ecliptic of 
the fundamental epoch on the equator of figure of dale and 
the equator of rotation of date rospcetively. Woolurd then 
decomposes the angles and into secular tenns 
V/pj>and and.pcrlodlc tenns ^iii//^A0^und with 


and correction terms 5^^ 60^. to give 


0 * + 5d 

0 r 


(Vll-‘)l) 


and 


2it“ i|/^ “ 2rr- 




(Vil.92) 


respectively, The angles 27 t- represent to the 

order of the approximations involved the combined effects of 
lunisolar precession (since they are referred to a fixed ecliptic) 
and lunisolar nutation in longitude since the fundamental 


2tr ' » 27 t - 1//* - 

Q » 0J + AO 

r r r 


(VM.93) 


ss 


The secular terms 2tr ' 0* niui 2n = are Itlentified by 

Woolanl with the luiiisolitr precession of the earth's figure axis 
unci roUillon axis icspectlvely. The periodic terms 
and = are Iclentificd with the lunlsolar nutation of the 

earth's figure axis and rotation axis respectively. 


Note; In his paper Woolard uses and e, to denote the 
angles 0^, and ho uses >lf a.’sd 0 to denote the angles ^^0^. 
Furthermore Woolard does not distinguish in his notation 
betsveen i/'p, and dousing d'. 0 for botlv. 


In Woolard's paper, Folsson's ecpiallons ap|>ear as equations 
(30) on page 34. The solutions of their descendenls, equations 
(44) and (45) on pages 47 and 4B, provide the functions \!»o, 
Oq. Hquations (1‘)) on page 24 govern the motion of the figure 
axis and their solution yields tire functions d'/. liquations 
(24) on page 26 govern the motion of the rotation axis and 
tlrcir solution yields the functions 0^, liquations (53) on 
page 131 govern the correction terms 8\jj^,d0p which yield 
the solutions from the solutions Oq.TIic solutions 
arc given by the expressions (54) on page 132. The 
equations governing tlic correction Icrms arc shown 

by Woolard to lire order of approximations Involved to be 
obtaiirablc directly from equations (53) for the correction 
terms d^^SO/' by multiplying by Ihc factor •* (C “ A)/A, 
Hence the solutions d\|/^6d^ can be obt.iincd directly from 
expressions (54) for by multiplying by - (C » ,4)/,4 

and arc given in expressions (55) on page 1 33. 

(_VU.‘)4) 


6d « - so 

»■ A f 


(Vli-95) 


epoch for tlicse expressions is also January 0 IbOOGrceiiwlcIi 
Meat! Noon (JD 2415020.0). Since tiie conlents of titls table 
arc refercnct'd to a moving ecliptic tlic cxpre.sslotts titerein 
reflect tire combined motion of tire ecliptic and equator since 
tlie epoelr. Tabic 26 represents llie effect of general precession 
and nutation wbiclr is tire smn of lunlsolar precession and 
nutation, wliielt irerturbs only tlic equator, plus planetary 
precession and nutation, wldcli perlnrbs only the ecliptic, 
Table 26 is generated from Table 24 by tlte addition of lire 
contents of Table 25, p, 152, wliicli contains tlte terms witicii 
■account for the planetary perturbations to the ecliptic. 

At the time of this writing the expressions for tlic mitatton 
of the rotation axis of Table 26 of Woolard’s worlv constllulc 
tlic lAU standard scries for tlie nutation and arc reproduced 
on pages 44 and 45 of tlie fuxphituiioty Supplcnmit to tlic 
Astronomical liplimcrls and the Amcrlant Kphmcrls and 
Nautical Almanac (/ 96 J), 

The earlli’s milation is driven by tlic gravitational torques 
of tlic sun and moon, and it is not surprising that the expres- 
sions for the nutation find tlielr most convenient representa- 
tion in terms of tlie angles ciij where 

CTj o C is llic mean anomaly of llie moon 

2%° 06' 16'/59 + 1325'' 198“50’ 56'.'79 7’ + 
33:'09r= +o'.'o5i8r^ 

c- C' is tlie mean anomaly of liie sun 

C' « 358°28' 33'/00 + 99'' 359“02' 59'.'10 7“ - 0'.'547’^ - 

o:’oi2or3 

Oj n p is the mean argument of llie latitude of tlie moon 
F * 11*’ 15' 03:'2- + 1343r82°Ql' 30'.' 547’- ll'.'SCr^ ™ 

0'.'0012r^ 


Presumably motivated by the conviction that llio motion of 
the rotation axis In space is tliat wiiich Is observable by 
astronomic means, Woolard tabulated the results of liis solu- 
tions for In - d*, and in two locations in his 

work. The first appears in Table 24 p. 138 ff, which gives the 
precession and nutation In longitude and obliquity of the 
rotation axis relative to the fixed ecliptic of January 0 1900 
Grcenwicli Mean Noon (JD 2415020.0). The epoch for these 
expressions is January 0 1900 Greenwich Mean Noon 
(JD 2415020.0). Since five contents of this table are refer- 
enced to a fixed ecliptic, the expressions therein reflect only 
the motion of the equator since the epoch, and as such 
represent only tlic lunisolar contribution to precession and 
nutation. 

The second tabulation .appears in Table 26 p. 153, wiilch 
gives Die precession and nutation in longitude and obliquity of 
the rotation axis relative to the moving ecliptic of date. The 


s /) is the mean elongation of the moon from the sun 
D « 350“44' 14'.'95 + 1236'’ 307°06' 5l!'l8r» S'.'nrH 

0'.'0068r^ 

ttj s is the mean longitude oflhe ascending lunar node 

a ^ 259°10'59'.'79 - S'" I34''08'3i:'23f + 7'.’482r* + 
O', '00807'^ 

and where T is measured in Julian centuries of 36525 days 
from the cpoclv January 0 1900 Grcenwicli Mean Noon 
(JD 2415020.0X 

Tlic nutation series is IradiUoiially tabulated in terms of 
amplitudes Ap*it }~ I, 2, 3, and argument 
coefficients Kj), (= t, 2, 3, ...5, in such a way lliai the 
nutation in longitude 6i^(7), reckoned positive to the east 
along the ecliptic, and the mUntion in obliquity, reckoned 
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poslUve if the obilqiilly is increased, are given by expressions 
of the form 


N r s 

/-I 


+ diurnal nutation In longitude (Vll*96) 


HT) - £ J)/m cos 1^2 *„o,(r)l 


+ diurnal nutation In obliquity (Vll*97) 


The diurnal motion in the nutation was shown In section V 
arid illustrated in Figure V*3 of this work to arise from tire fact 
that the axis of llgure and the axis of rotation in general 
depart from tlie axis of angular momentum L by the small 
angles nnd y respectively. Furthermore the Folnsot con- 
struction develo|X5d In section V and Illustrated in Figure V-5 
of the work shows that the figure axis and the rotation axis 
0) make nearly diurnal rotations In a space-fixed frame about 
the angular momentum axis moving on the surface of the two 
cones of a|jcx angles 2^ and 2y respectively. 

(h WoolanVs theory for the nutation of the rotation axis, 
Woolnrd's scries for the nutation of the rotation axis is given in 
Table VIM. 


Tabt* VH-1. Nutirtion mHm for Iht rotation axis (Wootard, 1M3) 


Index 

J 

Period. 

days 

“l 

«2 

“3 

“4 

“s 

Amplitudo 
(O'lOOOl arc) 

Amplitudo 
(0':00 (/l arc) 

^/l 


^/3 

^/4 

*/5 

1 

6798 

0 

0 

0 

0 

I 

-172327 - 173.7 f 

92100 + 9.17’ 

2 

3399 

0 

0 

0 

0 

2 

2088 + 0.2 r 

-904 + 0.47* 

3 

1305 

-2 

0 

2 

0 

1 

45 

-24 

4 

1095 

2 

0 

-2 

0 

0 

10 

0 

5 

6786 

0 

"2 

2 

-2 

1 

-4 

2 

6 

1616 

-2 

0 

2 

0 

2 

-3 

2 

7 

3233 

1 

-1 

0 

-1 

0 

-2 

0 

8 

183 

0 

0 

2 

-2 

2 

-12729 - 1.3 f 

5522 -2.9 r 

9 

365 

0 

1 

0 

0 

0 

1261 -3.1 f 

0 

10 

122 

0 

1 

2 

-2 

2 

-497 + 1.2 r 

216 - 0.6 r 

11 

365 

0 

-1 

2 

-2 

2 

214-0.5 r 

-93 + 0.37’ 

12 

178 

0 

0 

2 

-2 

1 

124 + O.ir 

-66 

13 

206 

2 

0 

0 

-2 

0 

45 

0 

14 

173 

0 

0 

2 

-2 

0 

-21 

0 

15 

183 

0 

2 

0 

0 

0 

16- O.ir 

0 

16 

386 

0 

1 

0 

0 

1 

-15 

8 

17 

91 

0 

2 

2 

-2 

2 

-15 + 0.17’ 

7 

18 

347 

0 

-1 

0 

0 

1 

-10 

5 

19 

200 

-2 

0 

0 

2 

1 

-5 

3 

20 

347 

0 

-1 

2 

-2 

1 

-5 

3 

21 

212 

2 

0 

0 

-2 

1 

4 

-2 

22 

120 

0 

1 

2 

-2 

1 

3 

-2 

23 

412 

1 

0 

0 

-1 

0 

-3 

0 

24 

13,7 

0 

0 

2 

0 

2 

-2037 - 0.2 f 

884-0.5 r 

25 

27.6 

1 

0 

0 

0 

0 

675 + 0.17’ 

0 

26 

13.6 

0 

0 

2 

0 

1 

-342- 0.4 r 

183 

27 

9.1 

1 

0 

2 

0 

2 

-261 

ii 3- o.ir 

28 

31.8 

1 

0 

0 

-2 

0 

-149 

0 

29 

27.1 

~1 

0 

2 

0 

2 

114 

-50 

30 

14.8 

0 

0 

0 

2 

0 

60 

0 

31 

27.7 

1 

0 

0 

0 

1 

58 

-31 

32 

27.4 

-1 

0 

0 

0 

1 

-57 

30 

33 

9.6 

-1 

0 

2 

2 

2 

-52 

22 

34 

9.1 

1 

0 

2 

0 

1 

-44 

23 

35 

7.1 

0 

0 

2 

2 

2 

32 

14 

36 

13.8 

2 

0 

0 

0 

0 

28 

0 

37 

23.9 

1 

0 

2 

-2 

2 

26 

-11 

38 

6.9 

2 

0 

2 

0 

2 

-26 

11 
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TaM*Vlh1.NiiMionMrlMfwl)wratolkmM(«(Wootordi 1WU)(C«f)MMMd) 


Index 

J 

Period, 

days 

Otj 

‘‘3 

Of- 

3 

®4 

“s 

Amplitude 

Af 

(010001 arc) 

Amplitude 

Of 

(010001 arc) 

^/1 

*/2 


*/4 

*/$ 

39 

13,6 

0 

0 

2 

0 

0 

25 

0 

40 

27.0 

-1 

0 

2 

0 

1 

19 

-10 

41 

32.0 

-1 

0 

0 

2 

1 

14 

■*7 

42 

31.7 

1 

0 

0 

-2 

1 

-13 

7 

43 

9.5 


0 

2 

2 

1 

-9 

5 

44 

34.8 

1 

1 

0 

-2 

0 

-7 

0 

45 

1.3.2 

0 

1 

2 

0 

2 

7 

**3 

46 

9.6 

1 

0 

0 

2 

0 

6 

0 

47 

14.8 

0 

0 

0 

2 

1 

-6 

3 

48 

14.2 

0 

-1 

2 

0 

2 

-6 

3 

49 

5.6 

1 

0 

2 

2 

2 

-6 

3 

50 

12.8 

2 

0 

2 

-2 

2 

6 

“2 

51 

14.7 

0 

0 

0 

-2 

1 

*•5 

3 

52 

7.1 

0 

0 

2 

2 

1 


3 

53 

23.9 

1 

0 

2 

-2 

1 

5 

-3 

54 

29.5 

0 

0 

0 

1 

0 

-4 

0 

55 

15.4 

0 

1 

0 

-2 

0 


0 

56 

29.8 

1 


0 

0 

0 

4 

0 

57 

26,9 

1 

0 

-2 

0 

0 

4 

0 

SB 

6.9 

2 

0 

2 

0 

1 


0 

59 

9.1 

1 

0 

2 

0 

0 

3 

0 

60 

25.6 

1 

1 

0 

0 

0 

-3 

0 

61 

9.4 

1 

-1 

2 

0 

2 

*-3 

0 

62 

13.7 

-2 

0 

0 

0 

1 

-2 

0 

63 

32.6 

-1 

0 

2 

-2 

1 

-2 

0 

64 

13.8 

2 

0 

0 

0 

1 

2 

0 

65 

9.8 

-1 

-1 

2 

2 

2 


0 

66 

7.2 

0 

-1 

2 

2 

2 

-2 

0 

67 

27.8 

1 

0 

0 

0 

2 

-2 

0 

68 

8.9 

1 

1 

2 

0 

2 

2 

0 

69 

5.5 

3 

0 

2 

0 

2 

-2 

0 


e, TIte ditifval nutation of the rotation axis. The diurnal 
motion of live rotation axis cu in space can be seen from Figure 
V4 to liavc an amplitude 7 wlrcrc Equation (VU-89) gives 

(VII.98) 

and whore is the amplitude of the Eulcrlan wobble, The 
angular rate of the motion of the rotation axis co about the 
angular momentum axis L in the spacc^fixed frame is given by 
Equation (V-52) as Wj + o^; namely, the sum of vire earth’s 
spin angular rate (diurnal rate) plus the wobble angular rate, 
Since the wobble angular rate is very small (435-day period) 
we see that this motion of the rotation axis in space has very 
nearly a diurnal period. 

Referring to Figure Vll*l and adopting the angular momen- 
tum equator as a slowly moving reference plane relative to 


inertial space we can see that the diurnal motion in space of 
the rotation axis co produces a diurnal nutation In longitude 
5 with amplitude given by 

<y>m 


and a diurnal nutation in obliquity with amplitude given 
by 

56,^ « y (ViMOO) 


where 0 is the obliquity of the ecliptic. To establish the phase 
of the diurnal nutations we refer to Figure V-7 and sec that 
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the body'llxed plane containing the vectors and u has a 
geographic east longitude \ where 

\ (VIMOl) 

and where specifies the body»fixed position 

of the Eulorian axis. The right ascension of this plane In the 
spacc*nxcd frame is \ + 0, where <p is Greenwich Sidereal 
Time, Referring to Figure VIH wo see that when 

\+^>‘0,n (VIM 02) 

the displacement of equivalently achieves its 

maximum negative and positive values respectively and Utal 
when 

(VIM03) 

the displacement y^ y^ or S vanishes. Hence we can set 

+ CVII-104) 

and by a similar argument we can set 

+ (VH-105) 

Now since < O', '2 arc and (C- A)/A 0.00328 we sec that 

taking 0 w 23® implies that the diurnal nutations in longitude 
of the rotation axis 6 are less than 

6t/'^^<0','0018arc (VII-106) 

and that the diurnal nutations in obliquity of the rotation axis 
arc less than 

5e^^<0:'0007arc. (VIM07) 

fr Woolard's theory for the nutation of the figure axis, 

Since even on a rigid earth the rotation axis w docs not remain 

fixed in a body-fixed frame, the motion of the rotation axis in 
space is not precisely shared by body-fixed axes such as 
interferometer vector baselines. In fact if the rotation axis is 
chosen as a reference axis for the nutation, then tlie motion in 
space of a body-fixed axis such as an interferometer vector 


baseline must be obtained from the motion of the reference 
axis by subtracting from the motion of tlie reference axis In 
.space the motion of the reference axis relative to tlie 
body-fixed frame. 

Alternatively on a rigid earth the figure axis is a 
body-fixed axis and Its motion in space is shared directly by all 
otlier body-fixed axes such as interferometer vector baselines. 
If the figure axis is chosen as a reference axis for the nutation, 
then the motion in space of any other body-fixed axis such as 
an interferometer vector baseline is identical to that of tlie 
figure axis, For this reason, In computing the evTcets of 
nutation on the orientation in space of an interferometer 
vector baseline it is useful to have a tlicory for the nutation of 
the earth's figure axis 

In Woolard's theory it is possible to generate a series for the 
nutation of the figure axis from ids series (Tabic 26) for the 
nuta/lon of the rotation axis by 

(1) Subtracting from the entries of Table 26 the terms of 
the equations (55) p. 1 33 of Woolard’s work. 

( 2 ) Adding (noticing the change in tire sign of 5 ^) to tlie 
entries of Table 26 the terms of tlie equations (54) 
p. 132 of Woolard’s work. 

Tlie results of this procedure arc presented in Table VII-z 
below, 

g, The diurnal nutation of the figure axis, The diurnal 
motion of the figure axis ^3 In space can be seen from Figure 
V-5 to have an amplitude 6^-7 wlicrc (3,. is the amplitude of 
the wobble and 

(VIM 08) 

Equation (V-18) and Figure V-3 sliow tliat the vectors co L 
and ^3 arc all coplanar and so (he rate of the motion of the 
figure axis ^3 about tlie angular momentum vector L in the 
space-fixed frame is identical to the rate of tlie rotation axis cj 
about the angular momentum vector L in the space-fixed 
frame and so is also given by Equation (V-52) as Wj + 0 ^. Just 
as in the case of the diurnal motion of the rotation axis in 
space this rate is tlie sum of the earth’s spin angular rate 
(diurnal rate) plus the wobble angular rate, and since the 
wobble angular rate is very smalt (435-day period) we see that 
this motion of the figure axis in space has a nearly diurnal 
period. 

Referring to Figure VIM and adopting the angular momen- 
tum equator as a slowly moving reference plane relative to 
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TaMiVH>2< Nutation mHM for lh*fl9Uir«M(*(Woo4Kd, INS) 


Index 

J 

I’eriQil. 

days 


“2 

«3 

*4 

H 

Amplitude 

Af 

(OlOOOl arc) 

Amptiiuda 
(OIOOOI arc) 


^/2 

^/3 

*/4 

*/5 

1 

6798 

0 

0 

0 

0 

1 

- I 72293. l -173.7 r 

92090+9.1 f 

2 

3399 

0 

0 

0 

0 

2 

2088 + 0.27' 

«904 + 0.4 f 

3 

1305 

-2 

0 

2 

0 

1 

45 

-24 

4 

1095 

2 

0 


0 

0 

10 

0 

5 

6786 

0 

-2 

2 


1 

-4 

2 

6 

1616 

-2 

0 

2 

0 

2 

»3 

2 

7 

3233 

1 

»1 

0 

-1 

0 


0 

8 

183 

0 

0 

2 

-2 

2 

-12804.5-1.37' 

5549.6-2.9 f 

9 

365 

0 

1 

0 

0 

0 

I261-3.ir 

-1.4 

to 

122. 

0 

1 

2 

-2 

2 

-501.4 + 1 . 2r 

217.6-0.67' 

11 

365 

0 

-1 

2 

-2 

2 

2 14 -0.5 r 

-93+0.37* 

12 

178 

0 

0 

2 

-2 

1 

124 + O.ir 

-66 

13 

206 

2 

0 

0 

**2 

0 

45 

0 

14 

173 

0 

0 

2 

..2 

0 

-21 

0 

15 

183 

0 

2 

0 

0 

0 

16 - 0.1 r 

0 

16 

386 

0 

1 

0 

0 

1 

-15 

8 

17 

91 

0 

2 

2 

**2 

2 

- i 5 + o.ir 

7 

18 

347 

0 

-1 

0 

0 

1 

-10 

5 

19 

200 


0 

0 

2 

1 

"5 

3 

20 

347 

0 

-1 

2 

-2 

1 

-5 

3 

21 

212 

2 

0 

0 

-2 

1 

4 

-2 

22 

120 

0 

1 

2 

-2 

1 

3 

-2 

23 

412 

1 

0 

0 

•*1 

0 

-3 

0 

24 

13.7 

0 

0 

2 

0 

2 

2l99-0.2r 

943.2- O.Sr 

25 

27.6 

1 

0 

0 

0 

0 

675 + O.ir 

=9.7 

26 

13.6 

0 

0 

2 

6 

1 

-375.5-0.4 r 

192.9 

27 

9.1 

1 

0 

2 

0 

2 

-291.8 

124.3-0.1 r 

28 

31.8 

1 

0 

0 

-2 

0 

-149 

-1.8 

29 

27.1 

-1 

0 

2 

0 

2 

118.6 

-51.7 

30 

14.8 

0 

0 

0 

2 

0 

60 

—1,6 

31 

27.7 

1 

0 

0 

0 

1 

8.6 

—31 

32 

27,4 

-1 

0 

0 

0 

1 

-54.2 

30 

33 

9.6 

-1 

0 

2 

2 

2 

-57.9 

24.1 

34 

9.1 

1 

0 

2 

0 

1 

-50.4 

24.9 

35 

7.1 

0 

0 

2 

2 

2 

-36.9 

15.8 

36 

13.8 

2 

0 

0 

0 

0 

28 

0 

37 

23.9 

1 

0 

2 

-2 

2 

-27.2 

-11 

38 

6.9 

2 

0 

2 

0 

2 

-30.1 

12.5 

39 

13.6 

0 

0 

2 

0 

0 

25 

0 

40 

27.0 


0 

2 

0 

1 

19 

-10 

41 

32.0 

-1 

0 

0 

2 

1 

14 

^7 

42 

31.7 

1 

0 

0 

-2 

1 

-13 

7 

43 

9.5 

^1 

0 

2 

2 

1 

-10.2 

5 

44 

34.8 

1 

1 

0 

-2 

0 

-7 

0 

45 

13.2 

0 

1 

2 

0 

2 

7 

-3 

46 

9.6 

1 

0 

0 

2 

0 

6 

0 

47 

14.8 

0 

0 

0 

2 

1 

-6 

3 

48 

14.2 

0 

-1 

2 

0 

2 

-6 

3 

49 

5.6 

1 

0 

2 

2 

2 

-7.2 

3 

50 

12.8 

2 

0 

2 

^2 

2 

6 

-2 

51 

14.7 

0 

0 

0 

-2 

1 

-5 

3 

52 

7.1 

0 

0 

2 

2 

1 

-6 

3 

53 

23.9 

1 

0 

2 

~2 

1 

5 

“3 

54 

29.5 

0 

0 

0 

1 

0 

-4 

0 

55 

15.4 

0 

1 

0 

-2 

0 

-4 

0 

56 

29.8 

1 

-1 

0 

0 

0 

4 

0 
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T»bt« VM*2» Nuta(k>n ••Htt for th# tigum omIi (Wooforti, !W)) (ContfmMd) 


imien 

J 

Petlod, 

ilnys 

er^ 

^/» 

*/i 

fta 

^/3 

fljit 

% 

Amplitude 
(OlOflOi arc) 

AmplUude 
(010001 arc) 

57 

26,9 

\ 

0 

“2 

0 

0 

4 

0 

58 

5.9 

2 

0 

2 

0 

1 


0 

59 

9.1 

1 

0 

•» 

m 

0 

0 

3 

0 

50 

25.5 

t 

1 

0 

0 

0 

•*3 

0 

51 

9.4 

1 

"1 

2 

0 

2 

*»3 

0 

62 

J3,7 


0 

0 

0 

1 

-2 

0 

63 

32.6 

“1 

0 

2 

»2 

i 

-2 

0 

54 

13.8 

2 

0 

0 

0 

1 

2 

0 

55 

9.8 

-1 

•"1 

2 

2 

2 

-2 

0 

55 

7,2 

0 


2 

2 

2 

-2 

0 

67 

27.8 

I 

0 

0 

0 

2 

-2 

0 

58 

8.9 

1 

1 

2 

0 

2 

2 

0 

69 

5.5 

3 

0 

2 

0 

2 


0 


Inertial space wo can see that the diurnal motion In space of 
the figure axis S'y produces a diurnal nutation in longitude 
5^^^ with amplitude given by 

(''"•'O’) 

and a diurnal nutation in obliquity 5c,,, with amplitude given 

/ « ■ 


6V,/ 0, ( 2 --j) (VIMIO) 

where 0 is tlio obliquity of the ecliptic. 

By referring to Figure VIM we sec that tlie nutation 
or equivalently is opposite In plrasc to the 
nutation 7, *“ 7/, or 6 as a consequence of tlie fact that the 
figure axis 6*3 and the rotation axis <0 lie on opposite sides of 
the angular momentum vector I,. Hence by arguments similar 
(0 those used in tlie case of the rotation axis we can establisli 
the phase of the diurnal nutations of the figure axis as 

(Vii-iii) 

and 

"/,/■■'’« (2-j) *>"W + V <V"-I12) 

and where C//1 ^ 1, 


/». Coinmems on Wooinrd's ilimy. it is apparent from 
Tables VII-1 and VII*2 tliat the nutations In Woolard's scries 
can be exliaustivoly grouped into nutations wliose period 
exceeds 91 days, which are normally referred to ns the long 
{period nutations, and nutations whose period is less than 35 
days, which are normally referred to as the short-period 
nutations. 

Apparent sidereal time Is referenced to the nutated 
equinox, If the rotation axis is chosen as the reference axis for 
nutation tlicn consistency requires that the rotation equator 
be used to provide tlie reference equinox 7, for apparent 
sidereal time. From Equations (VII-104) (VII-106) we see tliat 
this equinox suffers diurnal nutations in longitude of the order 
of O'.'OOiS arc wlilcli correspond to diurnal in 

equalities of apparent sidereal time of tlie order of 0,12 msec. 
Tills diurnal hi equality In apparent sidereal time Is usually 
ignored in most present-day applications of the theory. 

The motion W * m . t of the rotation axis w relative to 
tiie figure axis ^3 in tlie body-fixed frame is composed of the 
sum of 

(1) A geophysical or Eulcrian (torque-free) component 7n^ 
excited by a variety of internal processes on the earth 
and wliich is unpredictable on the basis of present-day 
gcopliyslcal knowledge, 

(2) A iion-Eulcrian lunisolar component +hMp 2 
excited by the gravitational torques^ of the sun and 
the moon acting on the cartii, often referred to as the 
“dynamical variation in latitude,” and which is pre- 
cisely predictable. Expressions for and arc 
given in Equations (V-94) and (V-95) respectively, 
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The Eulerlan motion has an ampiitudo of the order of Q"2 
arc maximum and the lunisolar motion has an amplitude of 
01*0178 arc maximum. Tlie lunisolar motion or tlie dynamical 
variations in latitude ti a retrograde motion with a nearly 
diurnal period, The center of its diurnal circular path 
corresponds to the Eulerlan imsltlon of tlie rotation axil in the 
hody.flxed frame as shown in Figure V*8, 

Tlie mean amplitude of the lunisolar motion or dynamical 
variations in latitude W, is roughly 0"0089 arc and so 
we see that In a hody.flxed frame the lunisolar motion of the 
rotation axis can he described as motion around a hody.flxed 
cone of mean apex angle 2(/Jp) O'/0 178 arc. In the absence of 
geophysically Induced polar motion or Eulerlan motion the 
precession of the earth's rotation axis w and figure axis in 
inertial space can bo visuallr.ed by a Folnsot construction In a 
manner described by Woolard on p. 31 and Illustrated below in 
Figure VII.2, 



Figure VH*2, The Poinoot conatructlon oppraprlete for Um doacrlp. 

tion of tho motion of tho oorth'o figuro Mie In In^M 
ipoco or lunlaolor proooMion, Tho booio vector it 
normei to tfie ecliptic, 0 io the meen obliquity of the 
ecliptic referenced to the equator of rotatlw, and 
{ p.) la tho meen amplitude of the funiaofir polar 
motion Iff^. 


The lunisolar body-fixed cone of mean apex angle 2<p^) 
rolls witliout slipping on the interior of a space-fixed cone 
centered on the pole of the ecliptic of apex angle 20 where 
0 is tp 1 obliquity of the ecliptic and witli the rotation vector 
CO occupying tlie line of contact between the two cones, Tlie 


retrograde lunisolar motion of tlie bodyfixed cone drives it In 
a retrograde sense around the Interior of the space-flxed cone. 
In the absence of geopliysically induced polar motion or 
Eulerlan motion the Eulerlan jmsltlon of tlie rotation axis In 
tlie hody.flxed frame coincides wltli the figure axlsf^. Hence 
In Figure Vlf.2 83 occupies tlie axis of the body-flxed cone. 
Taking 0 23.*473 (Stacey 1977) the diurnal rotation of the 
body»flxed cone transports and w around the space fixed 
cone In about 26,000 years, 

Figure VU.2 Illustrates the fact that during tills motion the 
angle between ^3 and 83 Is systematically smaller than the 
angle between i?3 and w by the amount (pp) 0*/0089 arc. 
This is tlie klnematlcal reason for tlie fact that the mean 
obliquity of the ecliptic for the figure axis ^3 is less than tlie 
mean obliquity of tlie ecliptic for the rotation axis <0. Tills 
kiiiematicai construction accounts for tlie existence of tlie 
constant term of -O'/OOSOS arc appearing In equations (54) 
p. 132 of Woolard’s work Involving the transformation from 
tho obliquity of the rotation axis to the obliquity of tlie figure 
axis. 

Tlie secular term of “0'.'00043T appearing in equations (54) 
p. 132 of Woolard *s work Involving the transformation from 
tlie node of tlie rotation equator to the node of the flgure 
equation has been iinown by Murray (1977) to be spurious. 
Klnoshlta etui, (1977) mid Kliioslilla (1977) have shown ho / 
the spurious secular term identified by Murray and otiiers like 
it appearing in Woolard's work arise from Incorrec* matlie. 
matical procedures. 

Woolard’s tlicory of tiie nutation is widely regarded as 
being Inadequate for present-day requirements, Fl^?;^f . ihc 
theory treats a rigid earth and there arc indications Ihuc iiiis 
restriction is responsible for errors in the predicted nutation of 
the order of 0','02 arc. This is insufficient accuracy for the 
requirements of modern observational techniques, In particular 
long baseline Interferometry, 

Furthermore It is apparent that the Instantaneous rotation 
axis w is not directly ob>''"i^ablc, Tills conclusion follows from 
the fact Unit any attempt lO observe the vector w necessarily 
requires observations extending over a finite interval of time, 
which for the classical methods Is typically several hours in 
duration and often as long as 12 or 24 hours, during which 
time the vector to will complete a significant portion of Us 
diurnal circuit in the body-fixed frame. The result of tliis Is 
that the axis which is observed in practice is not w but ati axis 
coinciding with some mean position of w averaged over the 
observing interval. 

I The revised theory of the nutation. In 1977 at tlie lAU 
Symposium No. 78 on Nutation and the Earth’s Rotation, a 
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working group was convened whoso task It was to revise tlio 
theory of the cartii’s nutation and to recommend a new series 
for the eartlt’s nutation to be adopted by the lAU at its 17th 
General Assembly in 1979 In Montreal, Canada. At the time of 
this writing this group lias completed its task and its 
recommendations to the lAU General Assembly are two*fold 
(J.C. Williams personal communication): (i)a change of 
reference axis, (ll) the computations of earth nutations for a 
deformable earth, 

/. The chanse of reference axis. Since all attempts to 
observe the present nutation reference axis Invariably yield 
a mean position for w, averaged over the finite observing 
internal in both the body-fixed frame and the space-fixed 
frame, It is recommended that the revised theory of the 
nutation make explicit recognition of tWs fact and replace the 
instantaneous axis o> with a new reference axis defined 
explicitly as the mean position of w when averaged over its 
predictable (l.c„ driven by external gravitational torques//) 
diurnal motion Itl , The new reference axis will be given a 
name to distingulsli It from others with wliich it might 
potentially become confused, Suggested names include “celes- 
tial reference pole,” "mean diurnal axis,” “Eulerian pole of 
rotation” and in Section II of this document it has been 
referred to as the “spin axis,” 


Tho new reference axis is most precisely described mathe- 
matically. Before doing so we require the result of a small 
It^snma: we wish to show that if the matrbcMis an orthogonal 
transformation matrix representing a spatial rotation then the 
matrix A/ M^’, where tiie dot denotes differentiation with 
respect to time, is antisymmetric, The proof of tills is as 
follows: 


The transformation of the components of a body-fixed 
vector denoted to the components denoted r^ of the same 
vector viewed in the space-fixed frame is 

jr, «TNSW£„ (VIM 13) 

* 

whore the matrix: 

P represents the precession of the reference axis of the 
body-fixed coordinate system. 

N represents the nutation of the reference axis of the 
body-fixed coordinate system 

S represents the net rotation (spin) about tho reference 
axis of the body-fixed coordinate system (UTl). 

IV represents the orientation of the reference axis within 
the body-fixed coordinate system (polar motion), 

In reality, observations arc made In a body-fixed frame of the 
body-fixed components of a vector (the direction to a 
quasar or a star) whose components in the space-fixed space 
frame are considered to be constant. 


= W'^ N'^ . (VIM 14) 

The body-fixed motion of the quasar or star is given by 


dt 


= Tibi!; s'^f^rpT ^ j^TpT + ^ pT 

[_ rfr dt dt 

+ ^1-. + A'’’ (VIM15) 


Since M is orthogonal M'^ - AT' and since AfAT* =/ 
where [ is the Identity matrix we haveMA/^ “/, Differentiat- 
ing with respect to time gives 

= 0 

Hence 

Therefore 

and AfA/^ is antisymmetric. Q.E.D. 


Note: Equations (VIM 13) and (VIM 14) deal explicitly with 
the components jTj of the vector and the components of 
the vector and so when differentiating with respect to time 
there is no need to distinguish between time derivatives taken 
m the body-fixed frame and time derivatives taken in the 
space-fixed frame. 

If we consider the case where drjdt ~ 0 which corresponds 
to a constant space-fixed direction and is achieved in practice 
by removing effects of aberration for sources exhibiting no 
proper motion, then Equations (VIM 13) and (VIM IS) give 

^ = \mL s^j^TpT ^^T^^TpT^ gT pT 
dt [_c/f dt dt 

+ IV^ N'^ ^ j PNSWr^^ (VIM 16) 
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which can be written as 



m 



W + w'r + It/'r ^ 5 w 

clt dt 


0^’ « P NS IV (Vll.l 26) 

tincl arc all antisymmetric (Goldstein 1950, pp. 124 fQ. 


+ IV'^’5'^’iV^' 


Ujjch term in the square bracket of Equation (VII *11 7) is 
either an antisymmetric matrix of the form AIM^' or is » 
similarity transformation on such a matrix* Since it is the 
property of similarity transformations that they preserve 
antisymmetry of matrices, we conclude tliat all four terms in 
the square brackets o,f Equation (Vll*117) arc antisymmetric 
matrices and so Equation (VU-1 17) can be written In vector 
form ns 


(It 


PNSW r 


I 


(VIM 17) 


dr, 

X rj (VIM 1 8) 


where 

(1) is a polar motion or wobble rate vector given by 




* 43-1 -^31 

? +d"’ e' 

‘■2 12 *'3 

(VIM 19) 

(2) 

Slg 

is a spin rate vector given by 




- 0^ e +0'^’ 
''23 ‘•1 "31 

'c 0^ ? 

<-2 ‘•3 

(V1M20) 

(3) 

«/v 

is a nutation rate vector given by 




a 0^ “e 
'^23 ‘'1 


(VIM2I) 

(4) 


is a precession rate vector given by 






(VIM22) 


and where the matrices 0^^ 0^ 0^ O'" are given by 

~ IV (VIM 23) 

qS ^ (V1M24) 

dt 

T* 

qN « ,p7’ ^ (VIM25) 
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The instantaneous rotation vector w is given by 

w « + n, + + $ 2 ^, (viM 27 ) 

and consists of the vector sum of the separate rotation rates 
due to polar motion, spin, nutation, and precession, The 
components of these separate rotation rate vectors in the 
body-fixed coordinate frame c, c, e. arc given by Equations 
(VIM19) -(VIM26). 

The magnitudes of the rotation rates in Equation (VIM27) 
are of the order of 


2X rad see*' maximum 

IT2J » 7,29 X 10*® rad see*' 

mArU«,,l « 7.9X I0*'2 rad sec*' 

and the vector rclationsldp of Equation (VIM 13) is illustrated 
in Figure VII-3. 



Figur* Vj|'3. Th« vector rvlationehipi ftortinont to oarth rotation. 

Tho inatantanoous axia u axacutaa a diurnal circuH in 
a body-fixad frama about tha Eularian rotation axia to.. 
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In Wcolfird's theory of the mitallon the vector w Is the 
rercrcncc iixls. The small vector + Up represents the 
diurnal dynamical vjirlallons in latitude. 

In the revised theory of the nutation the vector is the 
reference axis where 

« co»ajv ‘ O,,. (VIM 28 ) 

From their definition it can be seen that the vectors and 
Up arc rotating diurnally in the body-fixed frame and so Wq 
coincides with the mean position of co when averaged over a 
diurital circuit, The angle between <o and is of the order of 
0'.'02 are. 

The small vector represents the small rotation rate 
arisiitg as a result of the slowly varying position of to within 
the body-fixed frame due to geophysical causes other than 
external gravitational torques. The angle between the spin 
vector Ui^ and tOf^ is of the order of 5" X iO""’ arc maximum, 
and so we see that the reference axis of the revised nutntiotr 
theory lies very close to the spin vector. This is a convenient 
choice of reference axis for U coincides with die general 
conception of the earth's motion as consisting of a rapid spin 
about an axis which is in turn changing direction slowly in 
both (lie body-fixed frame and tire space-fixed frame. 

The new reference axis of the revised nutation theory also 
has the advantage of removing the diurnal oscillations of the 
position of the reference equinox on the fixed ecliptic, Since 
Wq coincides with the mean diurnal position of w it follows 
that the equator corresponding to Wq coincides with the mean 
diurnal rotation equator or the mean position of tiro diurnally 
varying rotation equator of the vccto- to (Figure VIM), If Yq 
is the ascending node on the fixed ecliptic of the fundaincrstal 
Of Dch on the mean diurnal rotation equator it is clear that Jq 
coincides with the mean position of and so can be identified 
With Jip Since the angular momentum vector L and the 
equinox are slowly moving in a space-fixed frame, wo see 
that by adopting the new reference axis we remove undesirable 
diurnal nonuntformitics of the order of ±0,12 milliseconds iiv 
the definition of apparent sidereal time, 

This principal advantage of adopting the new reference axis 
is that it is an observable axis. In a body-fixed frame it lies at 
the center of the quasicircuiar diurnal paths of the stars h the 
sky. Apart from the slow Eulerian motion the vector co^ is a 
body-fixed vector and its nutation in space will be nearly 
identical to the figure axis Cj, 

k. The nutation series for a deformable earth. A theory of 
the nutation for a deformable cartli has been developed by 


Kinoshita ct al, (197b) using the Molodcnski II earth model 
(Molodcnski 1961) and based on the accurate rigid cartli 
theory of the .v.v. tion of Kinoshita (1977) with appropriate 
modifications \v- dir ' for the clastic yielding of tlie earth. The 
earth model Molodciuki 11 has a liquid outer core and a solid 
inner core with radially varying clastic constants. 

Tlie modifications to the rigid earth theory are descrll’Ud 
briefly by Kinoshita et al. (1979) as follows, For each circular 
component of nutation with angular frequency 

tie, 

Nj * , (Vll-129) 

wlicro 


6, - £ 

Pi 

is the so-called argument of the nutation (Kinoshita 1977), a 
theoretical ratio, (n/flo)/Vy« is computed for the amplitude of 
the nutation for the defcrmablc earth to the amplitude of the 
nutation for the rigid earth for each nutation frequency A'^. 
This ratio is computed by a two-stage formula, 

p + J.7 (VIM 30) 

0.2159 - 100 

Nj- cOj 

= 1 +0.1124 [d- 4.1] (VIM 31) 

H/tV^ 

where 03 j is the angular rate of earth rotation about the figure 
axis and the numerical constants in those formulae are derived 
from the Molodcnski 11 earth model, 

Kinoshita ct al, (1979) then give the amplitudes of nuta- 
tions in longitude and nutations in obliquity for 
the Molodcnski II “rear earth in terms of the corresponding 
amplitudes of nutations in longitude and nutations in 
obliquity for tlie rigid earth as 



6S 


^ [( 4 / 

Using (his procedure KInoshita et ul. (1979) have gcnorn(cd 
the following lUitatlon scries for (lie figure axis of a deform- 
able earth wliose properties are (Imsc of the model “Molodcn- 
skill” This series lias been recomincndcfd (o the lAU for 
adoption at the 17th General Assembly In Montreal In August 
1979 and is presented In Table VU-3 below. The fundamental 
ojwch for tills scries Is 12000.0 [JED 24S154S.0J. The vari- 
able 7' Is measured in Julian cent uries from the epoch. 

I. Diurnal nuUUiom in the revised theory of the nutation, 
The diurnal motions of the figure axis in the revised theory of 
the nutation are essentially identical to the diurnal motions of 
the figure axis in Woolard’s theory of the nutation. This can be 
seen by combining Figures V-8 and Vll-2 and representing the 
hmisolar idotion of Figure Vll-2 in the presence of geophysi- 
cally induced polar motion or Eulcrian motion. This Is shown 
below in Figure VII-4. 

In this ease the body-fixed hmisolar cone Is centered on an 
axis occupying the Eulerian pole position, as can be seen in 
Figure V- 8 . The hmisolar body-fixed cone ofmean apex angle 
2(fip ) « o', '0178 arc rolls without slipping on the interfor of tlic 
space-fixed cone of apex angle 2 0 , where 0 is the obliquity of 
tiro ecliptic, Tlie rotation axis occupies the lino of contact 
between tlie two cones and the motiou, is retrograde as sliown 
in Figure V- 8 , The result of this motion is the steady retro* 
gr.ade progression of w about with a period of nearly 
26,000 years. 

Figure ViI4 indicates clearly tliat the figure axis F 3 moves 
in a prograde sense around tlie axis tOg as a rcisult of retrograde 
motion of w around the lunlsolar patii in the body-fixed 
frame, The figure axis ^3 is moving in a prograde sense on a 
cone of apex angle 2|3g wiicrc j3^ « 1171^1. 

It follows tiiat the formofee (VI 1-1 1 i and (V1MI2) will 
serve to de, scribe tlie diurnal tuitatlon in ><‘V 5 Hude and obliq- 
uity in the revised theory of the nutation at well, 

It is also of interest to observe that tlie reference axis oig is 
.slowly moving in tlie space-fixed frame as Indicated by Figure 
VIM. In Uie absence of EuJorian motion coincides with Cg 
in Figure Vfi- 2 , 



Fiflura VIM- The ipaoc-flxcd diurrai motion of ttw ooilti’o (Igum 
■xio wMofi occurt 00 a mouN of tfw ktnlaolar polar 
motion iff^ whonovor thorn la 0 nonvanlaMng gooptiya* 
k»Hy Induoad or Eulailan polar motion 


C. ThdEfftet of Solid EirthTiddt on Earth Rotation 

I. Tidal perturbations to the earth’s inertia tensor. Ele- 
mentary considerations are sufficient to sliow (Stacey 1977, 
pp. 90 ft) that the tld.il perturbing force per unit mass f,(r) 
acting throughout tlie body of the eartli due to a celestial 
body of niassAf at a distance /? from (lie earth’s center of mass 
can be expressed as the negative gradient of a tidal potential 
f/,(r)as 

f,(r) - VU,(t) (VI I -134) 


where 

«,W I) (''IM35) 

„S52 
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TaM* VH*3. Ntitailon MflM for «»• ngur* of ttw Motodanokl M daformaMa aarth modal <Klnoal>Ha a( al., ItTd) 


Index 

J 

Period, 

days 

"l 



®4 

®5 

Amiditude 

Aj 

(O'lOOOl are) 

Amplitude 
(O'lOOOl arc) 

^/1 

*/2 

*/3 

*/4 

^/S 

1 

6798.4 

0 

0 

0 

0 

r 

-l720S8-1742.2r 

92044+8.97’ 

2 

3399.2 

0 

0 

0 

0 

2 

2063+0. rr 

-895+0.57’ 

3 

1305.2 

-2 

0 

2 

0 

1 

46+0.07’ 

-24+O.Or 

4 

1095.2 

2 

0 

-2 

0 

0 

11 +o.or 

o+o.or 

5 

1615.7 

-2 

0 

2 

0 

2 

"3 + 0.07' 

1 + o.or 

6 

3232.9 

1 

-1 

0 

"1 

0 

-3 + 0.07' 

0 + o.or 

7 

6786.3 

0 

-2 

2 

-2 

1 

-2 + Q.D?* 

1 + o.or 

8 

943.2 

2 

0 


0 

1 

1 + 0.07’ 

o+o.or 

9 

1 82.6 

0 

0 

2 

”2 

2 

-13152 - 1.57’ 

5719-3.17’ 

10 

365.3 

0 

1 

0 

0 

0 

141|-3.4r 

49-O.ir 

11 

121.7 

0 

1 

2 

-2 

2 

-515 + 1.2f 

224 - 0.67’ 

12 

365.2 

0 

-1 

2 

-2 

2 

217-0,57’ 

-95+0.3r 

13 

177.8 

0 

0 

2 

-2 

1 

129 + O.ir 

-70 + o.or 

14 

205.9 

2 

0 

0 

-2 

0 

48 + 0.07’ 

0 + o.or 

15 

173.3 

0 

0 

2 

-2 

0 

-22 + o.or 

0 + o.or 

1(5 

182.6 

0 

2 

0 

0 

0 

17 -0,17' 

o + o.or 

17 

386.0 

0 

1 

0 

0 

1 

-15 + 0.07' 

8 + o.or 

18 

91.3 

0 

2 

2 

“2 

2 

-15 + 0.17' 

7 + o.or 

19 

346.6 

0 

-1 

0 

0 

1 

-12 + 0.07’ 

6 + o.or 

20 

199.8 

-2 

0 

0 

2 

1 

-5 +0.07’ 

3 + o.or 

21 

346.6 

0 

-1 

2 

-2 

1 

-5 + o.or 

3 + o.or 

22 

212.3 

2 

0 

0 

-2 

1 

4 + 0.07' 

-2 + o.or 

23 

119.6 

0 

1 

2 

-2 

1 

4 + o.or 

-2 + 0.07’ 

24 

411.8 

1 

0 

0 

-1 

0 

-4 + o.or 

0 + o.or 

25 

131.7 

2 

1 

0 

-2 

0 

1 + o.or 

0 + o.or 

26 

169.0 

0 

0 

-a 

a 

1 

1 + Q.or 

0 + o.or 

27 

329.8 

0 

1 

-2 

2 

0 

-i + o.or 

0 + o.or 

28 

409.2 

0 

1 

0 

0 

2 

I + o.or 

0 + o.or 

29 

388.3 

-1 

0 

0 

1 

1 

i+o.or 

0 + o.or 

30 

117.5 

0 

1 

2 

-2 

0 

-1 + o.or 

0 + o.or 

31 

13.7 

0 

0 

2 

0 

2 

-2260 - 0.2f 

972-0.5r 

32 

27.5 

1 

0 

0 

0 

0 

709 + 0.1 r 

-7 + o.or 

33 

13.6 

0 

0 

2 

0 

1 

-384-0.47’ 

199 + o.or 

34 

9.1 

1 

0 

2 

0 

2 

-299 + o.or 

128-O.ir 

35 

31.8 

1 

0 

0 

-2 

0 

-157 + o.or 

-1 +o.or 

36 

27.1 

-1 

0 

2 

0 

2 

123 + o.or 

-53+ o.or 

37 

14.8 

0 

0 

0 

2 

0 

63 + o.or 

-2 + o.or 

33 

27.7 

1 

0 

0 

0 

1 

63 + O.ir 

-33 + o.or 

39 

27.4 

-} 

0 

0 

0 

1 

-58- 0.1 r 

32 + o.or 

40 

9.6 

-1 

0 

2 

2 

2 

-59 + o,or 

25 + o.or 

41 

9.1 

) 

0 

2 

0 

1 

-51 + o.or 

26 + o.or 

42 

7.1 

0 

0 

2 

2 

2 

-38 + O.or 

16 + 0.07’ 

43 

13.8 

2 

0 

0 

0 

0 

29 + o.or 

-1 + o.or 

44 

23.9 

1 

0 

2 

-2 

2 

29 + o.or 

-12 + o.or 

45 

6.9 

2 

0 

2 

0 

2 

-31 + o.or 

13 + o.or 

46 

13.6 

0 

0 

2 

0 

0 

26 + o.or 

-1 +o.or 

47 

27.0 

-1 

0 

2 

0 

1 

21 + o.or 

-10+ o.or 

48 

32.0 

-1 

0 

0 

2 

1 

15 + o.or 

-8 + o.or 

49 

31.7 

1 

0 

0 

-2 

1 

-13 + o.or 

7 + o.or 

50 

9.5 

-1 

Q 

2 

2 

1 

-10 + o.or 

5 + o.or 

51 

34.8 

1 

1 

0 

-2 

1 

-7 + o.or 

o+o.or 

52 

13.2 

0 

1 

2 

0 

2 

7 + o.or 

-3 + o.or 

53 

14-2 

0 

-1 

2 

0 

2 

-7+ o.or 

3 + o.or 

54 

5.6 

1 

0 

2 

2 

2 

-8 + o.or 

3 + o.or 

55 

9.6 

1 

0 

0 

2 

0 

6 + o.or 

o+o.or 

56 

12.8 

2 

0 

2 

-2 

2 

6 + o.or 

-3 + o.or 

57 

14.8 

0 

0 

0 

2 

1 

-6 + o.or 

3 + o.or 
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Tab>»VII-a.NutoMofi— rte tforWwflgWMtoofttwIColodnaklN dafonniM«MrthnMiM(Klno«liHastal.i tt7t)(Con(lfHMd) 


Index 

J 

Period, 

days 

“l 

“2 

*3 

ot 

4 

“s 

AinpUuide 
(O'iOOOl arc) 

Amplitude 

lij 

(O'iOOOl are) 


*/2 

'/3 

*/4 

^/5 

58 

7.1 

0 

0 

2 

2 

1 

-7 + o.or 

3 + o.or 

59 

23.9 

1 

0 

2 

-2 

1 

6 + o.or 

-3 + o.or 

60 

14,7 

0 

0 

0 

-2 

1 

-s + o.or 

3 + o.or 

61 

29.8 

1 

-1 

0 

0 

0 

5 + o.or 

o+o.or 

62 

6.8 

2 

0 

2 

0 

1 

-5 + o.or 

3 + o.or 

63 

15.4 

0 

1 

0 

~2 

0 

»4 + o.or 

0 + o.or 

64 

26.9 

1 

0 

-2 

0 

0 

4+o.or 

0 + o.or 

65 

29.5 

0 

0 

0 

1 

0 

-4 + o.or 

o+o,or 

66 

25.6 

1 

1 

0 

0 

0 

~3 + o.or 

0 + 0.07’ 

67 

9.1 

1 

0 

2 

0 

0 

3 + o.or 

o+o.or 

68 

9.4 

1 

-1 

2 

0 

2 

«3 + o.or 

1 + o.or 

69 

9.8 


-1 

2 

2 

2 

-3 + o.or 

i + o.or 

70 

13.8 

-2 

0 

0 

0 

1 

-2 + o.or 

1 + o.or 

71 

5.5 

3 

0 

2 

0 

2 

-3 + o.or 

i + o.or 

72 

7.2 

0 

-1 

2 

2 

2 

-3 + o.or 

1 + o.or 

73 

8.9 

1 

1 

2 

0 

2 

2 + o.or 

»i+o.or 

74 

32,6 

“I 

0 

2 

*•2 

1 

“2 + o.or 

1 + o.or 

75 

13,8 

2 

0 

0 

0 

1 

2 + O.or 

-1 + o.or 

76 

27.8 

1 

0 

0 

0 

2 

-2 + o.or 

1 + o.or 

77 

9.2 

3 

0 

0 

0 

0 

2 + o.or 

0 + o.or 

78 

9.3 

0 

0 

2 

1 

2 

2 + 0.07' 

“1 + o.or 

79 

27,3 

-1 

0 

0 

0 

2 

1 + o.or 

»i + o.or 

80 

10.1 

1 

0 

0 

-4 

0 

-1 + o.or 

0 + o.or 

81 

14.6 

-2 

0 

2 

2 

2 

1 + o.or 

-1 + 0.02' 

82 

5.8 

-1 

0 

2 

4 

2 

'- 2 + o.or 

1 + o.or 

83 

15.9 

2 

0 

0 

-4 

0 

-1 + o.or 

o + o.or 

84 

22.5 

1 

1 

2 

*•2 

2 

1 + o.or 

«•! + o.or 

85 

5.6 

1 

0 

2 

2 

1 

-1 + o.or 

1 + o.or 

86 

7.3 

-2 

0 

2 

4 

2 

“1 + o.or 

1 + o.or 

87 

9.1 

-1 

0 

4 

0 

2 

1 + o.or 

o+o.or 

88 

29.3 

1 

-1 

0 

-2 

0 

1 +o.or 

o + o.or 

89 

12.8 

2 

0 

2 

-2 

1 

1 + o.or 

»i +o.or 

90 

4.7 

2 

0 

2 

2 

2 

-1 + o.or 

o+o.or 

91 

9.6 

1 

0 

0 

2 

1 

-1 + o.or 

o + o.or 

92 

12.7 

0 

0 

4 

-2 

2 

i+o.or 

o + o.or 

93 

8.8 

3 

0 

2 

-2 

2 

i + o.or 

0 + o.or 

94 

23.8 

1 

0 

2 

-2 

0 

-1 + o,or 

o + o.or 

95 

13.1 

0 

1 

2 

0 

1 

1 + o.or 

o+o.or 

96 

35.0 

-1 

-1 

0 

2 

1 

i+o.or 

o + o.or 

97 

13.6 

0 

0 

^2 

0 

1 

-1 + o.or 

o + o.or 

98 

25.4 

0 

0 

2 

-1 

2 

-1 + o.or 

o + o.or 

99 

14.2 

0 

1 

0 

2 

0 

-1 + o.or 

0 + o.or 

100 

9.5 

1 

0 

-2 

--2 

0 

-i+o.or 

o + o.or 

101 

14.2 

0 

-1 

2 

0 

1 

-1 + o.or 

o + o.or 

102 

34.7 

1 

1 

0 

-2 

1 

-1 + o.or 

o + o.or 

103 

32.8 

1 

0 

^2 

2 

0 

-1 + o.or 

0 + o.or 

104 

7.1 

2 

0 

0 

2 

0 

1 + o.or 

0 + o.or 

105 

4.8 

0 

0 

2 

4 

2 

-1 + o.or 

o + o.or 

106 

27.3 

0 

1 

0 

1 

0 

1 + o.or 

o + o.or 


and where 

(1) r= \r\<a, 

(2) 0 is the zenith angle of the celestial body as seen by 
the observer at r. 


(3) ®) ^1*® Legendre polynomial of degree n. 

(4) G is the gravitational constant. 

Figure Vll-S illustrates the geometry of the situation. The 
geographic coordinates of the observer arc Oq Xq and tire 
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FIgurt Vll*5. Tht daflnHkma of Iho goomofrlcrl quaniHlM uaod In 
tho dtvtiopinoni of Hm thoory of tlM solid sorth tidos. 

geoeraplUc coordinates of a line joining the earth’s center of 
mass and the celestial body of mass Marc 

The addition theorem for spherical harmonics allows us to 
write the expression for the tidal potential t/,(r) in terms of 
the geographic coordinates Qq ^q,0„ X„, , 

»i>„(cos0„,)i>„(cos0o) 



f;;Xme,)cosm%-\) (VIM 36) 

For the sun and the moon (r//?©) <<1 and (r/Ji,^)<< l 
and sufficient accuracy is usually acliicved in bquation 
(VIM35) by truncating the sum at « - 2. 

(/ M ^ /» (003 0) ,t<a, (VIl-137) 


For a celestial body such as the sun or the moon whose 
celestial coordinates are changing relatively slowly we sec tlial 
the argument ’ Xq) is periodic In an Interval of approxi* 
malcly wr • days. The term m * 1 gives rise to the diurnal tidal 
variations and the term m * 2 gives rise to the semidiurnal tidal 
variations. 

Tho tidal deformations of the earth cause perturbations 
to the elements of the earth’s Inertia tensor which can be 
calculated using MacCullagh’s formula. Written out in full 
Equation (VII* 138) becomes 

+ I F*(cos OJ F‘(cos Oq) cos(X,„ - \) 

+ ^ />5(cos Fj(cos Op) cos 2(X,„ » Xp)j , 

r < (I , (VIM 39) 

This expression can be rewritten 

y,w - ■ 

+ I Pl(cos 0,,,)P^(cos Oq) cos X^,, cos Xp 
+ “ 7’‘(cos 0,„) /’’(cos Op) sin X,„ sin Xp 
+ ^ /’=(cos 0„,)/’5(cos Op) cos 2 cos 2 Xp 
+ ^ /’5(cosO,,,)/’5(cosOp)sln2X,„ sin2Xpl , 


which in geographic coordinates Op Xp> 0,„ X,„, becomes 
GMr^ 


UAt) 


R 


1 

m 


-P,(cosO,,,)i’2(cos0a) 


m»i ^ ^ 


cosmfX^^j - X, 


o)]> 


r <a . (V1M38) 


r < a , 

wldcli in turn can be cast in tlie form 


(VII.140) 


* 2 (^) ^2 '(cos Op) fc"' cosmXp 

«i*0 ^ 

+ 5'"sinmXo j, r < fl , (V1M41) 


et 


where 


POOR QUAWlTlf 


CO * . ^ po^cosOJ 
5 ® - 0 

Cl .„^/>»(cosO,„)cosX,„ 


The above five equations In six unknowns arc supplemental 
with the additional equation 

f.ikJe . 5 (VU.148) 

where 5 Is the tldally induced variations In the trace 

of the inertia tensor of the earth* 




S5 - 


S «• 


a m 


GMr^ 

3K® 


GMr^ 

12/? 3 



/»* (cost?,,,) sin X,„ 
Pj(cos <),„) cos 2 X,„ 
/’5(cost),„)sin2X,„ 


(Vli-142) 

it tiien follows from Section VI-2 and specifically from 
Equations (VI-74) - (Vl-78) that 

yMe - — P^(cost? )sin2X... 

' 6/?® ^ 


It lias been shown by Darwin (1910) that for an incompres- 
sible earth in which all possible deformation fields ii(r) have 
the property tliat V • u ® 0, the trace of tlic inertia tensor is 
preserved under earth deformations. In his analysis of the 
effects of earth tides on the earth’s rotation Woolard (1959) 
assumed that the earth was incompressible in order to use tlio 
property that in his solution. l£lcr Rochester 

and Smylie (1974) siiowed that the value of 7,,7was preserved 
under a much wider class of deformation fields than those for 
wlUcii V » u « 0. In particular tlicy show that even on a com- 
pressible earth for wliich V * u ^ 0 tlie deformation field 
arising from tlie effects of tidal perturbations preserves tlie 
value of T^, Following Rocliester and Smylie we can, without 
any restrictive assumptions, take 

ytidc ^yUda ^yMc « 5 (r,./)'"*' = 0 . 

(VII-149) 

Equations (Vll-143) “ (VlI-149) have as their solution 


(VH-143) 


2kJ1a^ 

2^//da ^ — />l ({JOS 0 ) COS X 

3/?® ■ 


(Vll-144) 


')Mde 
"'a 3 


2kMa^ . 

(cos 0,^,) sin \ 


3R^ 


m 


(VII-145) 


yMC « _ p2 QJ sJ,, 2X,„ 



(VII-ISO) 



'■f? = ^l(cosO,„)cosX/M 

(Vll-151) 

kMa^ 


" 3 p 3 (cost?,,) sin X/n 

(VIM 52) 


jide , 
^^2 


.tide 


T) 


k.Ma^ 

3 (cos 0,„) cos 2hn 

6/?3 2 

(Vll-146) 


tide 

II 


k^Ma^ 

6R^ 


[2P® (cos 0,„) 2 (cos 0,„) COS 2XwJ 


(Vll-153) 


+ ^2^ - 2rjf 


2k^Ma^ 

R^ 


P% (cos 0,„) . 


(VIM47) 


.tide 

22 


k^Ma^ 

6R^ 


[2Pj (cos 0 ,„) -Pi (cos 0,„) cos 2Xm] 


(VIl-154) 
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,ritk H 

'^33 


20f«® „ 


(VIM 55) 


Using GqiiuUon (V1I*140) In Equation (V-U460) and Integral* 
Ing over the variables f, X gives 


For the sun and the moon the angles namely , X r , 
arc varying dlurnally with a period near I day. Hence r^, 
are dlurnally varying products of Inertia wliich are 
IHjrlodIc In about 1 day. The moments of Inertia 

have long period components which depend on the angles 
, namely Oq^O^, which arc In fact given by 


Oe 

0* -90“-6.j 


(VIM 56) 


N, 


2n GMa^ 
5 j^3 



iW\ (cosO) 

“rfT^ 


(10\ 


(VIM6I) 


The Integral In the above equation vanishes with the result that 


N, « 0 (VIM 62) 

and so in the absence of dissipation the angular momentum L 
of the earth is conserved under the action of the tidal forces. 


where 6„, 6, are the deeUnatioes of the sun and moon Ifthe unperturbed earth is described by an Inertia tensor /„ 
respectively. ^ 


2. The effect of the solid earth tides on UTt. It is relatively 
easy to show that in the absence of Internal dissipation which 
introduces phase lags in the tidal response of the earth the 
tidal forces exert no net torque on the earth. 


/ * 
'o 


A 

0 

Lo 


0 

A 

0 


0" 

0 

c. 


(V1M63) 


The net torque N, on the earth due to the tidal forces f,(r) 
is 

N, =/ tXi,(x)dV (V1M57) 

V 


and a rotation vector where 

cOq “ nFj (VIM 64) 

then the tidally perturbed earth is described by an Inertia 
tensor Twhcrc 


where the integral is taken throughout the volume of the 
earth. Ifthe vector field f,(r) has components 

f,(r) = /,,(r)rV,o (0? +/,x«'X (V1M58) 

then 

N, = / Ko(r) X“ r/,x (0^1 dV (V1M59) 

V 

Using Equation (VIM 34) in Equation (VIM 59) this becomes 



A + 

^tide 

12 

-tide 

'13 

1 - 

ytUe 

'21 

A 1* r 
/I rfjj 

-tide 

Jr3 


ftidc 

'31 

-tide 

'^32 



(VIM65) 

and a rotation vector <o where 

oi » Ej Cj +(1 +m^‘')c3l 

(VJM66) 


N. 




L' 30 slnO 3X 


(VIM 60) 


Conservation of angular momentum under the influence of 
the tides gives 


7q-co^=7'-w (VIM 67) 
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which gives rise to tiie titreo equations 


S2((/J + (I -t-wf")! » 0 

(Vll.168) 

O [rj^f m + (A + rj^O 0 + "‘f 0] * 0 

(VU»169) 

n[rjf wif" + r 'f /nf " + (C + /-*“") (1 + wifOJ * flC 


(VIW70) 


is a measure of the amount by wlricii tiie earth's rotation lags 
behiiul that of a liypothctlcal uniformly rotating eartir. Tiie 
tidal contribution to r(0 denoted it) Is 


« - f mfHt')clt’ seconds, (VIM 78) 
From liquations (VIM55), (V1M76) (VII478) we have 

T^^^io * ^ r 


—V—/’; (cos 5, „(('))*' 


(VIM 79) 


for the tiuce unknowns ^ ,«i 3 ^‘'.To first order In 

small quantities Equations (VIM68) <=» (VIM 70) reduce to 



(VIM 71) 


(VIM72) 

^(C + On;"*-’ tr;^") = «C. 

(V1M73) 

The first-order solutions for the unknowns 
are then 


it 

(VIM74) 

-tide 

,„tlde _ „ 

'"a A 

(VIM 75) 

yMc 

"a c 

(VIM 76) 


The tidal variations in UTi are dependent on tiie history of 
the quantity . It is customary to express the variations in 
the earth’s rotation rate by parameter r(0 where 


Now 

(cos 0,,,) * (3 cos^ - 1) (VIM80) 

and 


% “90''-6„, (VIM81) 


where is tire declination of the celestial body of mass Af. It 
follows that Equation (VIM 82) can be written 

'o 

(VIM 82) 


Woolard (1959) Integrated the above Equation (V1M82) 
for the combined effects of the sun and the moon. 




-^tl-3sin^S^(r')l 



[1-3 


sln^%(f')]U/' 


(VIM 83) 


t(0 « - J m^if) (It' seconds (VII* 


177) 


to obtain an expression for T"‘*%t) wliicli includes both the 
lunar and the solar tides, 
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Woolard’§ expression for Is given In the form 
20 

m sin milliseconds, 

H 

(VII- 1 84) 

whore the ejwch (q Is 1900.0. Equation (Vil-184) can bo 
evaluated from Table V1I4, taken from Woolard’s work. The 
tabulated amplitudes arc in milliseconds and the vaittg of the 
second degree Love number can be taken to be 

m 0.29 . 


directly from Equations (VIM51) (V1I.152) and (VIM74) 
(VIM75), which taken togethirglve 


m 


tide 


(VII.185) 



MLpi 


(cosd,„)sln\„ (V1I.186) 


where M Is the mass of the perturbing body and R is Its 
geocentric distance and where 0„, X,„ arc the geographic 
coordinates of the position vector R. 


Table VIM. Tba alwnanta UMd In Equation (VIMI7) to oanorata Hw 
tbaoratical contribution to t*** « (UTt-UTC)*** duo to tba affacta of 
tba aolld aarth tidao, (Altar Woolard, itst). 


/ 

Ap 

msec 

0/ 

Period, 

dnvs 

1 

0.32 

2U+S 

9.1 

2 

0.13 

2L +g - ii 

9.1 

3 

2,47 

21 

13.7 

4 

1.02 

2 L-n 

13.7 

5 

0.10 

2L-2n 

13.7 

6 

O.ll 


13.8 

7 

0,23 

2b ^20 

14.8 

8 

2.63 

e 

27.6 

9 

0.17 

g + n 

27.6 

10 

0.17 

g-n 

27.6 

11 

0.14 

2L-^ 

27.1 

12 

0.06 

2L-/r-n 

27.1 

13 

0.58 

2L"X“ 26) 

31.8 

14 

0.60 

20+rr' 

122 

15 

15.29 

20 

183 

16 

0.37 

20 -n 

365 

17 

4.88 

4'* 

365 

18 

0,23 

20 -g' 

6793.7 (18.6 year) 

19 

515.0 

n 

3396.9 (9.3 year) 

20 

2.7 

2n 



The arguments 4>f of the sine function arc defined by taking; 
L mean celestial longitude of the moon 
S mean anomaly of the moon 

S2 mean celestial longitude of the moon’s ascending node 
O mean celestial longitude of the sun 
g mean anomaly of the sun 

3, The effect of the solid earth tides on polar motion. The 

effect of the solid earth tides on polar motion can be deduced 


Due to earth rotation the argument X,„ decreases by 2n In 
slightly more than one sidereal day, allowing for the eastward 
progression of the tidednducing body, be It the sun or the 
moon. 


The complex quantity given by 

« ,ft rWe + / „,rWc (VIM87) 

defines the angular motion of the tidal perturbations to the 
rotation axis In the body.fixcd frame. Equations (VIM85) 
(VIM 86) (VIM 87) together gives 

0idc « . t MIsLpt (cos OJ (VII- 1 88) 
2 2 m 


and since X,„ continuously decreases (moves continuously 
westward) we sec that Is a retrograde motion of the 
rotation axis as shown in Figure Vn*6, 

Using the formula P\ (cos 0„) » cos sin 0^ and 
recognizing that the coordinate is related to the 
declination of the celestial body by 

" 90® » 
m m 


we can rewrite Equation (VIM91) as 

fjjtide cos 6^ sin , (VIM 89) 
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FIguraVM. TTw r«fct>a f dR polaf moWon W** f—u>Hog from tktol 

pvfiUfV^PQViv fO ww vsnn •mwmiOTmf* 

A general formula which Includes the tidal perturbations to 
the rotation axis of both the sun and the moon can be written 
directly from Equation (VIM 89) as 


ffjilde , < 


Mr. 




Q 


i?3 

\ *'0 


cos Sq sin ^0 e ^ 


^aa^i cos 6(j sin 6(^ e 


(VII-190) 


where: 

(1) Mq Mf^ arc the masses of the sun and moon 
respectively. 


D. Tht OMiphytIcal CKcHaUofi FunctlofM for Polor 

Motion ondUTI VaHoliont 

As pointed out previously in this work one of the 
important objectives of obtaining precise measurements of 
variations In UTl and polar motion on the earth is to ’learn 
more about the global geophysical processes wlilch are 
responsible for the variations. In order to carry out such a 
program It will be necessary to use existing geophysical 
knowledge to generate realistic tlieoretical rotation excitation 
functions whose predicted consequences for UTJ and polar 
motion variations can be compared against observed data. 
Suclt a program will require as Inputs, in addition to precise 
polar motion and UTi observations, considerable global 
synoptic data concerning the state of the earth’s oceans and 
atmosphere as well as information concerning the Internal 
state of the earth, 

Our ability to generate models for the oceans and 
atmosphere of the earth has advanced dramatically in the last 
decade With the development of earth satellites capable of 
monitoring the global state of the earth's atmosphere and 
oceans. In addition the ability to theoretically model a number 
of important Internal processes in the earth such as cartirquake 
faulting and to accurately deduce their effect on the earth’s 
inertia tensor for the case of a realistic oartli has advanced 
considerably in the past decade. These trends will no doubt 
continue Into the future and the forthcoming data will provide 
the basis for generating a priori rotational excitation functions 
which could bo rellncd by the precise UTl and polar motion 
measurements. 

Our analysis of the geophysical excitation functions begins 
by decomposing the density field of the cartli p(r,r) Into a 
mean density p®(r) and a geophysical perturbation to the 
mean Ap(r, r). Tlie position vector r refers to a fixed position 
in the rotating geophysical coordinate franie 

p(r,r) « p° (r) + Ap(r, t) (VIM9 1) 


(2) Rq R(i arc the geocentric distances to the sun and 
moon respectively, 

(3) 60 arc the declinations of the sun and moon 
respectively. 

(4) * a© “ G.A.S.T. and X(j; ■ «([ - G.A.S.T. where 
«0 arc the riglit ascensions of the sun and moon 
respectively and CAST refers to Greenwich Apparent 
Sidereal Time. 

The maximum combined amplitude of polar motion due to 
the lunlsolar solid earth tides is of the order of 6" X 10"^ arc 
or 18 cm of motion, The motion is retrograde with a nearly 
diurnal period. 


The inertia tensor of the earth I is similarly decomposed 
according to Equation (lV-5) into a mean inertia tensor r® and 
a geophysical perturbationT 

r- 7® +r (VIM92) 

where 


/ = 



p(r, 0 (r^ 7- r r J r/K 



'’ja 

'■13 

^21 

^+'■22 

'*23 

- '‘31 

'*32 

C+r 


(VIM 93) 
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ai 



p^(r)[r^ 7*“ rr 


(/K • 


’/I 

0 

,0 


0 

A 

0 


0 ' 

0 

c 

$>* 


(VI1.I04) 


f 


^la 

'*13 

T" 1 Ap(r, 0 r** r r ) dV <• 

^71 

1^22 

^23 


/at 

^32 

^33. 


(Vn.195) 


From HquiUioii (IV*22) wo see that tlie geophysical 
excitation function for polar motion 1(0 Is given by 

where 

F*ri3 +ff33 

/?-/», +///^ (V1M97> 

iV^ » iVj + f 


and where 


N • iV,F, +yV 3?3 

h * //jF| 


(VI1.I98) 


arc the external torque and relative angular momentum 
respectively measured in the rotating geophysical coordinate 
frame. The relative angular momentum Ii Is given by Equation 
(111-52) as 


h (r, 0 “ I p(r» 0 [r X V (r, /)] (IV (Vll.199) 

where v(ri/) is the velocity of the material of density p(r,r) 
relative to the rotating geophysical coordinate frame. 

From Equation (lV-30) we sec that the geophysical 
excitation function for UTl is 



(VII.200) 


and S 0 i apart from the question of the externally applied 
torques N, the excitation of polar motion and UTl nuctua- 
tlons depend entirely on the tensor Tand the vector h and 
their time derivatives taken In the rotating frame. 


In order to simplify our dynamical theory (Equaiio.n 
(IIMO)) we have chosen to define “the earth" to include Its 
oceans and atmosphere and according to this formulation of 
Its rotational dynamics the motion of the oceanic currents and 
atmospheric circulation by being part of "the earth" are 
Incapable of exerting an "external" torque on the earth 
through some sort of viscous boundary layer Intcrnctlon with 
the solid surface. The effect of the oceanic and atmospheric 
circulation on the earth's rotational dynamics is entirely 
included In the relative angular momentum term li. 


At the cost of complicating the dynamical description of Its 
rotation we could have defined "the earth" to exclude the 
oceans and atmosphere. In this case the motion of the oceans 
and atmosphere do exert an external torque on "the earth” 
through viscous boundary layer Interaction. In addition to 
being dynamically disadvantageous this formulation of the 
problem of the earth's rotational dynamics Involves the poorly 
understood phenomenon of the boundary layer Interaction of 
the oceans and atmosphere with the solid earth and with each 
other by requiring that we model this process In order to 
express the oceanic and atmospheric torques on "the earth” in 
terms of their respective velocity fields. While these two 
approaches to the problem arc formally cqulvalenii the 
definition of "the earth" (0 Include the oceans and atmo- 
sphere is dynamically simpler and has been shown (Lambcck 
and Cazenave, 1973) to bo capable of a more accurate 
treatment of the effects of oceanic and atmospheric circula- 
tion. Similar remarks could presumably be made for the fluid 
motions of the eartlTs fluid core and its effect on the earth’s 
rotational dynamics Is complicated by the fact that no direct 
measurements of the fluid velocity of the core are presently 
possible. 


St Is our present objective to obtain explicit formulae for 
the contribution to ilie rotational ev, citation arising from 
internal geophysical processes and will not be concerned at 
this point witli developing detailed expressions for the 
lunisolar gravitational torques N. It is our objective to obtain 
expressions for the terms /‘jj ^33 h\ h<x fi^ 
dri^jdt dhildt dh^ldt which appear as the internal geo- 
physical contributions to the rotational excitation functions in 
Equations (VU-196) and (VII-200) In terms of observable 
global flclds such as mass density, velocity, mass displacement 
etc. 
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If Urn geocfiiurlc posUlon vector r In the rotating geophysl* 
cal reference frame Is 

r (VII.201) 

then we have from Equation (VII- 186) 

r.- f (VIi-202) 

Jy 

and from Equation (VII499) 

h, - f pih t) r. Vfe (r, () dV (VII-203) 

Jy 

where 6/y Is the Kroneckcr delta defined by 

1 0 /•jfe/ 

1 i */ 

and where Cjjif Is the alternating tensor dellttcd by 

! +l //k cyclic 123 

0 (/fe not all distinct 

~ 1 ilk noncyclic 123 

The time derivative dhidt, reckoned In the rotating frame, 
Is given by Equation (111-84) as 


ill 

dt 




where Equation (111*81) gives 

>*(riO * P(r./) IrX v(r,0] 
and Equation (111*82) gives 

^rot ^ ^ O (r. 01 

It follows that 


flfh 

d( 




which can be written 
r/h 


dt 


/|[|^ + V‘(pv)j |rXv) 


+ pv* V(rX v)| dV 

which by virtue of the continuity equation 
*1^+ V . (pv) » 0 

reduces to 




rX 


.^1 

DfJ 

(VI 1-204) 


(VII-205) 


^ I p [r x|^] + pv • V (r X V) I rfK . (Vn-206) 
To obtain the components of dUldt we write 




Ar 3V 


dV 


dh 

1ft 




dV 


dh 

dt 




dhf 

IF 


dv, 


r f / -k 


)h 


Since 


this reduces to 


V''/ ’ ° 


- 

dh, f r 


I dV 

dt " J [_ 



•'V L 

' /« / J 
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Recosjnlzlng 


wlilch can be written as 


3 9 0 

IT 


as the Lugranglan lime derivative we conclude that 


h r 

h / >‘Vi-Dr<‘^ 

1/ 


+ V • (pv)] (r’ 7*“ r r ] 


+ pv • 7 T- r r )> dV 


or In coordinate free notation, 




In addition to we also require expressions for dFjdt, 
From Equation (1V*5) we have 


dt dt dt 


(VII-209) 


and since 


wliich by virtue of tlie continuity equation, Equation (HI-46), 
reduces to 


J pv Vlr^r-r r]rfF,(VU.214) 
In component form Equation (Vll-214) gives 

5 '//-O''; I ‘'I' 

!/ 


we have 


dt~dt‘ 

From Equation (HI-85) and (VII-21 1) we have 


(VH-210) 


(VII-21 1) 


which finally gives 




(VH-212) 


where is given by Equation (HI-83) as 


(r» 0 = P (f . 0 V (r, /) (r^ r- r r ) . 


(VIl-213) 


It follows that 


= f P{2/-^ v^5//-r-, V/-r,v^}r/F. (VI1-21S) 

•/ 1/ 


From Equations (VII-195) (VH-I99) (VlI-207) (VlI-215) 
we see that the relevant quantities appearing in the geophysical 
excitation functions are 


= - f Ap(r,t)rtr^dV (VII-216) 

•/i/ 


Apit,()r.r^dV (VH-217) 


= f (VIF218) 

Ju 


= J r- r 9r J f V • (pv) [r* T~ t r] 

+ pv* V [r^i'-r r]| dV 


ft I = / P(J>‘ 

Jy 


0 ('■2V3(f,‘)"''3V2(r,0J dV 
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dr 


* f p(r, /) (rj, V, (r, 0 « r, v^Cr, 01 ci V 

Jy 

(Vll.220) 

/»3 * / P(r, 0 (r» 0 " ''2V,(r, /)] d^ 

Jy 

(VI 1.22 0 

* " / Pift t) [r,V3(r, t) +/J Vj(r, 0] dV 

J U 


SubsUtuUng Equations (VII.2I6).. (VII-225) Into Equa- 
tlons (VIM96) (Vl(.197) (Vll.200) givce expressions for the 
geophysical excitation funcUons for 

(1) Polar Motion 
1 




An 




/ Ap(r,r3-/r,r3)f/K 


J3 

dt 


dr. 


23 

d7 


» - / P(r, 

1/ 


(Vll-222) 

0 [Va^^O + r-gVjfr,/)] dV 

(V11.223) 


fjp [^^Vi+rj ("5^ +WV,) 

'■j (it ■^'"'''' 3)] 


fh 

dt 


^■r, (—" +v vv,) tv 


(VI1.226) 


dj^ 

dt 


«^p(M)[r2(^+vW3) 

- rj + v • VVj)] dV. (VII- 224 ) ( 2 ) UTl 

“ J pit, 0 jr-j + V • V V, j ^ ~ ^ fy ^ 

+ VV 3 )] rfK. (V1I.225) - • (ViI-227) 
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